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SUMMARY OF PROJECT 

The Ontario Ministry of the Environment (OME) awarded a contr- 
act to Concord Scientific Corporation to calculate puff statistics from 
ensembles of individual trajectories for use in the OME Statistical Long 
Range Transport (LRT) Model. 

The first phase of the work program consisted of examination 
of several coordinate systems and methods for computing the trajectory 
statistics such that distortions due to projecting a spherical surface 
(the surface of the Earth) onto a plane would be minimized. A method 
and a coordinate system which would potentially meet these objectives 
were developed, but could not be fully tested due to time and funding 
constraints. It was decided to compute the statistics in the OME 
modified CMC grid coordinates due to practical considerations. The 
first part of the work program was reported in an interim report 
CSC 110.J271.01WP, 

The second phase of the work program consisted of generating 
trajectories using the OME Lagrangian Trajectory Model and the data for 
all 12 months of 1978, 1979, and 1980, and January and July of 1982. 
These trajectories were then used for calculating the statistical para- 
meters to be used in the Statistical Model. Samples of the resulting 
statistics for 1978 - 80 have been reported in the main text of the 
Final Report. 
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Two sets of trajectories and their statistics were generated 
for 1982, namely for surface and 850 mb levels, for the purpose of 
comparing the performance of the Lagrangian Model for these two levels. 
At the same time, a thorough examination of parts of the Lagrangian 
Model Code was made for the purpose of providing an independent review. 
No significant problems were found as a result of this examination and 
comparison. It was found that in some cases 850 mb trajectories move 
farther away than the surface trajectories and in a more intuitively 
expected direction. This part of the work program is reported on in an 
ADDENDUM to the Final Report. 
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1, INTRODUCTION 

The Final Report describes in Section 2 the data used in 
generating the trajectories using the OME Lagrangian Model and in 
Section 3, the formulae used in calculating the statistics. 

In Section 4 sample results for the years 1978-80 are present- 
ed. Throughout the report emphasis is given to the theory and codes 
developed particularly for this project. The complete set of results 
are stored in data files on the OME's computer system. These files are 
described in the SPECTRUM documentation which have been provided as part 
of the final deliverables for this contract. Based on the examination 
of the sample results as well as the work reported in the ADDENDUM, 
recommendations are given in Section 5. 
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2. TRAJECTORY GENERATION AND MASS BUDGET CALCULATION 

Trajectories have been generated using the OME Lagrangian 
Model, which uses time series of surface pressure data as input. These 
surface pressures are first fitted to IDth degree polynominals. The 
model then calculates wind fields and corresponding trajectories origi- 
nating from selected source points. A more complete description of the 
model is given by Ellenton et al . (1982) and also in the ADDENDUM of 
this report. 

For this project, surface pressure data were sorted and fitted 
to polynominals, and trajectories were generated every three hours for 
the first twelve hours after release, and every six hours until ninety- 
six hours from release (or until loss of the trajectory either due to 
the domain boundary or other computational reasons). For this work, 
computer codes which were developed by the Ministry were used. Trajec- 
tories from sixty sources designated by the Ministry were generated for 
a period of three years (1978-80). These sources are listed in Table 

2.1 and shown in Figure 2.1. Individual trajectories for the months of 
January 1978 and July 1978 for Sudbury, Ontario, and for CONMJO ALLEG 
near Pittsburg, Pennsylvania, are plotted in Figures 2.2 - 2.5. The 
extent of dispersion of trajectories is evident from these figures. In 
Figure 2.6, trajectories for January 1980 for Sudbury are plotted, but 
only until each trajectory reaches the shore. A comparison with Figure 

2.2 indicates how loss of trajectories due to a boundary would affect 
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the mean location of trajectory end points. In Figure 2.6, the earlier 
loss of trajectories in the east would bias their mean location towards 
west. 

Trajectory statistics for only three of the sixty sources have 
been reported here in detail as sample results. The complete set of 
results are stored on the Ministry's computer system. 

Once the trajectories are generated using the kinematic 
trajectory module of the model, they are then weighted according to how 
much pollutant they carry at any given time using the mass budget 
calculation module of the model. This module takes into account the wet 
and dry deposition and chemical transformation of a given pollutant, and 
requires a time series gridded precipitation data. For this project, 
the precipitation data supplied by the Ministry have been sorted and 
gridded using the Ministry's computer codes. 

All the data files generated during the course of the above 
work have been delivered to the Ministry earlier, as each data file was 
created and finalized. 

The trajectory end-points calculated by the Lagrangian Model 
are given in OME modified CMC coordinates. A description of this coor- 
dinate system can be found in Alp and Moran (1983). The weight associ- 
ated with each trajectory end-point given by the mass budget calculation 

3 
module is given as a concentration of the pollutant in yg/m . 
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Table 2.1 The Source Points Used in the Project. 
The coordinates are given in OME modified CMC coordinates, 
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Figure 2.1 The Source Points Used in the Project 
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Figure 2.2: Individual trajectories for January 1978 initiated 

at Sudbury, Ontario. H 
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Figure 2.3: Individual Trajectories for July 1978 initiated 
at Sudbury, Ontario 
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Figure 2.4: Individual trajectories for January 1978 initiated 
at CONMJO ALLEG near Pittsburgh, Pennsylvania 
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Figure 2.5: Individual trajectories for July 1978 initiated 
at CONMJO ALLEG near Pittsburgh, Pennsylvania. 
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SOURCE STATION: SUDBURY ONT 



DATE.- JAN 1980 



3 * S • ? t • JO 11 12 13 1< IS J« ir U It ZO 21 22 23 2« 25 26 27 2« 2S 30 3. 32 33 3« 3S 3e 3- 




S 6 



1 — I — ; — 1 — : — i — 1 — : — i — : — i — : — : — i — i — i — i — r 

10 11 12 13 M IS )6 17 1< IS 20 21 ZZ Z3 24 2S 26 27 21 29 30 31 

CMC X CO-ORD 



32 33 34 35 36 3" 38 3S 



Figure 2.6: 



Individual trajectories for January 1980 initiated at 
Sudbury, Ontario. These trajectories have not been 
continued beyond the shore. 
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%■ ESTIMATION OF STATISTICS 

The main objective of this project was to estimate the first- 
and second-order moments of the x and y coordinates of points in trajec- 
tories generated by the Ministry's Lagrangian Long Range Trajectory 
Model. These statistics were to be estimated at discrete travel-times 
from ensembles of trajectories from the same source that have been 
initiated at different times. 

Before writing the mathematical formulae that were used, which 
may look unfamiliar at first glance, the formulae for the estimates of 
mean x of n observations x . , the standard deviation o , and the covari- 

1 A 

ance o of two random variables (x and y) are written below: 



- i=l ' 



' I (x, - X ) -^ 



= /M_: 




X 



n 

2 (x- - x) (y. - y ) 

and 0. - 



xy n 



Note that is the biased estimator of the variance rather than the 



x 

unbiased estimator 
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^ (x. - x)^ 

s^^M-^ 

n-1 



which tends to be used more frequently. 

For the statistics to be estimated, it was required to compute 
the means, standard deviations and the covariance of the coordinates of 
trajectory points. It was also required to calculate mass-depletion 
weighted statistics (or "decaying puff" statistics) for sulphur dioxide 
and sulphate. This is done by using species concentrations at each 
point in the trajectory as a weighting function. 

To clearly state the formulae for the statistics, let us care- 
fully define the terms that will be represented. At a time-step t from 
a source k, the statistics to be estimated are: 

a) the mean of the x-coordinates of the points, which is repre- 
sented by X (k,t) , 

b) the mean of the y-coordinates of the points, which is repre- 
sented by y (k,t), 

c) the standard deviation of the x-coordinates, which is repre- 
sented by o (k,t), 

d) the standard deviation of the y-coordinate, which is repre- 
sented by (k,t), and 

e) the covariance of x- and y-coordinates, which is represented 
by0j^y(k,t). 
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We suppose that N(k,o) trajectories are initiated from source 
k, and that an integer from 1 to N will be assigned to uniquely repre- 
sent the time each trajectory is initiated. The x-coordinate of the 
point- at time t from the i-th trajectory initiated from source k is 
denoted by x.(k,t), i=l,2....,N(k.t). We define y.(k,t} similarly. To 
denote the weight assigned to the point of the i-th trajectory for 
source k at travel-time t we use w.{k,t). This weighting variable 
w.(k,t) is useful in the unweighted case as well. We define w.(k,t) to 
be a binary weighting function in the unweighted case. That is, either 
w.{k,t) = 1, or w.(k,t) = 0. Since trajectories can be terminated pre- 
maturely and not reach the full anticipated lifetime, in the unweighted 
case where t is the first time step where there is no point, then we 
set w.(k,t) = for all t>t . In the cases where a species concen- 
tration is used as a weighting function, we similarly set w.(k,t) = for 
t>t* (this is reasonable since if there is no point, one cannot calcu- 
late a species concentration). Let W(k,t) denote the accumulated weight 
from all trajectories that are initiated at source k after a travel time 
t. The formulae for estimating the statistics are as follows: 



N 
W{k,t) = I w.(k,t) 
i=l 



N 

x(k,t) = ^— I x.(k,t)w.(k.t) 

W(k,t) 1=1 ^ 



N 
y(k,t) = ^ ^^ y^{k,t)w.(k.t) 



W(k.t) 
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X 



/' \ ^ 2 

^ • ' i=l 



'y^^^^^ " ^ -m^ } [yi(l<.t)-y(k,t)]^w.(lc,t) 



and 



1 ^ 
°xy^^''^^ ' ^nrrn J [xT(k,t)-x(k,t)][y.{k,t)-y(k,t)]w.(k,t) 

These six equations are those used in the unweighted, sulphur dioxide 
weighted, and sulphate weighted cases. 

Five more statistics have been calculated for various initiation 
periods and sources, but these additional statistics are derived from 
the preceeding six. The percentage of initial concentration remaining 

is given by: 

The interpretation of a(k,t) in the unweighted case is the 
percentage of trajectories surviving at travel time t that were 
initiated from source k. The mean concentration is given by: 
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where N(k,t) is the number of trajectories surviving at travel time t 
that were generated from source k. That is, N{k,t) is equal to W(k,t) 
of the unweighted case. There is no point in reporting b{k,t) for the 
unweighted case since it is identically equal to one except when there 
are no trajectories and then there are no statistics at all. The stan- 
dard error of the x-coordinate is 



a (k,t) 
^ / W(k,t) 



and correspondingly, the standard error of the y-coordinate is 



(k.t) 



y / w(k.t) 

It might need to be pointed out that the standard error is an expression 
of the variability of the mean, whereas, the standard deviation is an 
expression of the variability of a single observation. 

If one assumes that the points in the ensemble of trajectories 
N{k,t) are distributed according to the bivariate normal (or Gaussian) 
distribution, then standard statistical analysis leads us to the 
estimate 0(k,t) of the orientation of the major-axis of the probability 
ellipse with the x-axis that is given by 
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(k,t) 

e{k,t) = arctan [-^ ] 

o^ (k.t) 



where.© is reported in degrees and is constrained such that 
-90" < 0(k,t) < 90*. 

A discussion of the estimation of statistics is not complete 
without documentation of the precision of the floating-point representa- 
tion of the values that are calculated by the required computer 
routines. Due to time constraints it was not feasible to perform a 
thorough and detailed analysis of error resulting from floating-point 
representation and floating-point arithmetic. Nevertheless, based on 
established practice and a quick check, we were able to obtain a 
floating-point representation in which we are confident. The computer 
routines were written in the FORTRAN language which allows only three 
different precisions of floating point representation: single precision 
(REAL*4), double precision (REAL*8), and quadruple precision (REAL*16). 
For a number represented in scientific notation (for example, 6.735 x 
10"^), single precision would give 7 significant decimal digits for the 
mantissa, double precision would give 16 significant digits, and quad- 
ruple precision would give 35 significant decimal digits for the 
mantissa. 

Tables 3.1, 3.2 and 3.3 contain statistics estimated for the 
first 10 trajectories in January 1978 that were initiated at Sudbury. 
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Each table was generated using the same program except Table 3.1, which 
was obtained using single precision arithmetic. Table 3.2 was obtained 
using double precision arithmetic, and Table 3.3 was done using quad- 
ruple- precision arithmetic. These tables were produced by a develop- 
mental precursor of the production program and contain a few more 
statistics not calculated in the production runs. These will now be 
described for completeness. The correlation coefficient for two random 
variables is 



X y 



I — 
where S^^ = /^^ I {x. - x) 



n 

I 
i=l 



/' 1 



n 



S = / -^ y (v.- y) 

J T=l 



1 " 



The measure of skewness used is 



Bi = 



^l(Xi-xr 



1 TT 
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The measure of kurtosis used is 



e^ = 



i Kx, - -x)' 



2 ^ 

X 



Some entries in the tables are error flagging values such as negative 
standard deviations (to indicate the variance was calculated to be nega- 
tive - a result of floating-pointing arithmetic error), -11,111 or 
-88,888 to indicate that a division by zero occurred in calculating the 
statistic, or -77,777 to indicate a nonpositive variance was produced in 
the algorithm that estimates the correlation coefficient. 

Comparison of the three tables shows discrepancies between 
single and double precision arithmetic. These discrepancies are signi- 
ficant enough to recommend that single precision not be used. Standard 
statistical computing practice is not to use single precision, so the 
recommendation is not a surprise. The comparison of the eight decimal 
digits in the mantissa for the statistics produced by double precision 
and quadruple precision show consistency for at least the first and 
second order moments. This is much better than the single and double 
precision comparison where at times only one of eight digits are in 
agreement. 
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Tables 3.4 and 3.5 list estimates of statistics for all trajec- 
tories initiated in January of 1978 from Sudbury for double and quad- 
ruple precision, respectively. The values produced by the two preci- 
sions for the first and second order moments are identical. The third 
and four order moments by the two precisions agree to at worst five 
significant digits. 

Therefore, on the basis of this brief analysis, it was decided 
that double precision would be fully satisfactory for our purposes. 
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Table 3.1: Estimated statistics for first 10 trajectories for Sudbury. Single precision, 
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Table 3.2: Estimated statistics for first 10 trajectories for Sudbury. Double precision. 
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Table 3.3: Estimated statistics for first 10 trajectories for Sudbury. Quadruple precision. 
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Table 3.4: Estimates of statistics for trajectories for Sudbury. January 1978. Double precision, 
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Table 3.5: Estimates of statistics for trajectories for Sudbury. January 1978. Quadruple precision. 
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4» RESULTS 

The trajectories are initiated on a tnonth-by-month basis for 
all sixty sources beginning in January 1978 and ending in December 1980. 
So the basic data unit for this project is a set of ensembles of trajec- 
tories initiated in the time-span of one month. 

The data graphically presented in this report are abstracted 

from the data base developed as the principal task of this project. 

This data base stores W(k,t), x(k,t), y(k,t), o (k,t), o„(k,t), and 

X y 

a (k,t) for all pairings of sources k and travel-times t for three 
Ay 

cases: unweighted, sulphur dioxide weighted, and sulphate weighted. 
Results for all three cases using the period comprising 1978 to 1980 are 
tabulated in Tables 4.1, 4.2, and 4,3 for the sources, Sudbury, Baldwin, 
and CONMJO ALLEG near Pittsburgh, Pennsylvania, respectively. Also 
listed in these tables are percentage weight remaining a(k,t), the mean 
concentration b(k,t) {where applicable), and the angle of orientation 
©(k.t). 

Similar tables for the same three sources have been produced 
for the ensemble sets of 1980, 1979, 1978, January, April, July, 
October, January 1978 and July 1978. These are listed in Appendix A, 



h^ r'. ) Concord Scientific Corporafioi 



- 26 - 



The mean trajectory for the month of January 1978 is plotted 
in Figure 4.1 for the source points located at Sudbury, Ontario, at 
Baldwin, Illinois and CONMJO ALLEG near Pittsburgh, Pennsylvania. The 
plotted squares indicate the location of the sources. The other graphic 
symbols indicate the points in the trajectory at each time step with 
line segments joining the centres of the symbols so that the time evolu- 
tion is clear. Figure 4.2 Illustrates the mean trajectories from the 
month of July 1978 for the same three sources. In Figure 4,2 we see the 
eastward movement that had also been expected for January 1978 but did 
not occur. When examining these figures, one should keep in mind that 
the actual dispersion about these means is very large, as has been indi- 
cated by Figures 2.2 to 2.6. Also, the percentage of trajectories 
remaining in the calculation of the statistics (see the corresponding 
tables in Appendix A) should be taken into account in the interpretation 
of the figures. 

A logical period for which to estimate statistics consists of 
the saffTie month in all three years. Figure 4.3 illustrates the mean tra- 
jectory for each of the aforementioned sources using trajectories from 
the month of January in all three years. The mean trajectories for 
January of these three years do not differ greatly from the mean trajec 
tories for January 1978 which are depicted in Figure 4.1. The mean 
trajectories for the three sources estimated from trajectories generated 
in July of 1978, 1979 and 1980 are shown in Figure 4.4. The differences 
between these mean trajectories and their counterparts in Figure 4.2 for 
July 1978 are quite small. 



^ 



Concord Scientific CorDoratior 



T 



27 



The comparison of the January and July mean trajectories for 
the three sources shows stark differences for January and July. The 
July mean trajectories move eastward whereas the January trajectories 
move westward. When the tables in Apendix A corresponding to the 
Figures 4.1-4.4 are examined in conjunction with the figures, it is seen 
that the January patterns start deviating from their original courses 
when the number of trajectories remaining in the calculation of the 
statistics drop below about 95% of their original number. This usually 
happens between 12 and 36 h after initiation of the trajectories. This 
observation suggests the importance of loss of trajectories in estimat- 
ing trajectory statistics. 

A primary goal of this project was to estimate statistics from 
as large an ensemble of statistics as is possible. This is done by 
using all the trajectories generated from January 1978 to December 1980 
for each source. Figures 4.5 through 4.8 depict mean trajectories from 
each source based on this period in groups of fifteen so that all the 
mean trajectories from the sixty sources are illustrated. The most 
striking feature about these mean trajectories is that they appear to 
converge towards the centre of the calculational domain. It can be seen 
from Tables 4.1-4.3 {% remaining) that loss of trajectories due to 
boundary effects may play a significant role in this behaviour. 



^) 



loncord Scientrfic Corporatio' 



28 



Figures 4.9, 4.10 and 4.11 are two dimensional plots of the 
time evolution of the standard deviation of both x-and y-coordinates 
from the full period set of 8,652 trajectories for Sudbury, Baldwin, and 
CONMJO ALLEG near Pittsburgh, Pennsylvania, respectively. The figures 
all look yery similar. For all three sources, the standard deviation of 
the y-coordinate exceeds the standard deviation of the x-coordinate for 
the first hours until sometime between 24 h and 42 h of elapsed travel 
time, when the situation is reversed and remains so until the last hour 
of travel time (96 h). The time evolution of both x and y standard 
deviations is smooth. 

For each source, either standard deviation rises quickly at 

first from zero at initiation time but a virtual plateau is reached at 

about 48 h of travel -time where the value remains between 2.0 and 2.5. 

A cursory examination shows that for t < 36 h, both a and a are appr- 

X y 

oximately proportional to /T. This implies that the constant eddy 
diffusivity assumption used in the derivation of the Gaussian distribu- 
tion which forms the basis of the OME statistical model is physically 
realistic. For larger times, the limiting domain size may be the reason 
for the a's becoming constant. 
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Table 4.1: Statistics for unweighted and sulphur dioxide weighted trajectories . Source - Sudbury 1978-80. 



SOUMCr STATION ISl 
fIJME MCIGHT 



I - ONt SUOrtURY 
« RCHtlNlNG 



f» 


8*52, 


100.0 


3. 


a«sz. 


100. t 


-i. 


SA52. 


loo.o 


«, 


SASZ. 


100.0 


it* 


aiKSz. 


100. • 


I*. 


BA23. 


M.T 


m* 


B«S2. 


98.8 


M. 


8191. 


97.1 


»* 


8705. 


9«.8 


«£• 


r9«8. 


92.3 


«i. 


TTI8. 


•9.S 


S*. 


r«ii. 


n»v 


M. 


T|T5. 


•i«t 


M, 


«9«3. 


?«.■ 


ft. 


6m3B. 


,f».» 


n. 


6183. 


T»*» 


•«. 


61 3S. 


fi*« 


«•• 


58T7. 


if.* 


mm 


S«3S. 


*S.I 



PEDIOOl TSBO 


KCIGHTi src 


UNKEIGHlEn 






1 •*£*»• 


l-ittH 


ISIGHA 


V516M* 


COW»Oi»>«CE 


25.790 


# 


0.0 


15.560 






0.0 


0.0 




25.837 


• 


0.002 


15. S3* 






0.201 


0.221 




25.886 


• 


0.00<i 


15.510 






0.396 


0.*3S 




25.935 


« 


0.006 


15.*85 






n.582 


0.6*1 




25.982 


• 


0.008 


IS.*6I 






0.757 


0.837 




26.0Sfc 


• 


0.011 


I5.kll 






1.0*3 


1.18* 




26.098 


• 


0.0I« 


15.381 


t 




1.3*0 


l.*67 




26.108 


• 


0.017 


15.358 






1.585 


1.677 




26.065 


. 


0.020 


15.331 






1.79* 


1.8** 




26.0«2 


* 


0.022 


15.302 






1.968 


1.966 




2S.9TB 


• 


0.02* 


15.273 






2.106 


2.050 




25.906 


t 


0.026 


15.2*5 






2.225 


2.099 




25.8«0 


r 


0.027 


15.237 






2.327 


2.1*3 




25.781 


• 


0.829 


15.239 






2.396 


2.161 




25.7)1 


s 


0.030 


15.227 






2.*6* 


2.168 




25.692 


• 


0.031 


15.232 






2.509 


2.173 




25.669 


♦ 


0.033 


15.273 






2.SS8 


2.182 




25.631 


♦■ 


0.03* 


15.319 






2.606 


2.190 


o!n6 


25.593 


♦ 


0.035 


15.350 


• 




2.65* 


2.207 


0.»37 



7t«T» 



•O.f 

».? 
«.» 
s.s 
s.« 

«.!• 
f.« 

«.I 

• .T 
f.« 
T.I 
«•■ 
«.« 
3.S 
3.ff 
2.* 
».> 
l.« 



U3 



hk 



•• 


2*eAl.lD2 


»• 


20*57.111 


t» 


17535.56* 


9!» 


15088. 2«l 


la* 


13022.326 


ifr. 


9738.8*7 


2«« 


7266.802 


W* 


5379.035 


3A* 


3969.850 


«t* 


2922.902 


^mM. 


21*2. **6 


i*« 


1569.600 


«9. 


11^1,972 


•ft. 


8*2. 3SS 


T!f. 


616.625 


fi» 


*50.0*9 


M* 


310.827 


«t. 


2*2.603 


ffft. 


178.01* 



ONT SUOOUffT 


PERIOOt 7880 


HEIGHT! S^C 


5LW.WJ9 


01 OH IOC HEIGHTCD 


ftlNlNO 


Me«<4 CONC. 


(MCAN 


r>4e*N 


I5IGHA 


VSIGM* 


COVARU 


180.0 


2.781 


25.790 


♦ 


0.0 


15.560 


I 0.0 


0.0 


0.0 


0.0 


85.0 


2.36* 


25.832 




0.001 


15.527 


• 0.002 


0.209 


0.727 


O.OOS 


72.9 


2.027 


25.879 




0.003 


I5.*9| 


> 0.003 


0.*08 


0.**6 


0.017 


62.7 


1.7** 


25.925 




O.OOS 


15.*52 


t 0.005 


0.598 


0.655 


0.032 


5*.l 


1.505 


25.970 




0.007 


15.*I0 


• 0.007 


0.777 


0.85* 


o.oso 


*0.5 


1.129 


26.038 




o.eii 


15.320 


• 0.012 


1.09* 


1.705 


0.103 


30.2 


0.850 


26.072 




0.016 


15.2*7 


> 0.017 


1.385 


l.*89 


0.208 


22.* 


0.6*1 


26.063 




0.022 


15.185 


> 0.023 


1.6*9 


1.708 


0.3*7 


16. S 


o.*s* 


26.015 




0.030 


15.111 


>. 0.030 


1.878 


1.878 


0.*79 


12.1 


0.366 


25.93* 




0.038 


15.0*3 


t 0.037 


2.06? 


2.00* 


0.556 


8.9 


0.278 


25.825 




0.0*8 


l*.982 


L 0.0*5 


2.211 


2.089 


0.5*5 


6.5 


0.211 


25.713 




0.059 


1*.935 


. 0.05* 


2.3*0 


2.138 


0.500 


*.8 


0.161 


25.598 




O.072 


l*.92l i 


. 0.065 


2.**3 


2.190 


0.363 


3.5 


0,122 


25.*89 




0.086 


1*.926 J 


. 0.076 


2.509 


2.217 


0.181 


2.6 


0.093 


25. *0* 




0.10* 


l*.909 . 


. 0.090 


2.587 


l,f>h 


8.050 


1.9 


0.071 


25.309 




0.123 


I*.9|0 < 


. 0.10* 


?.6U 


2.256 


-8.020 


1.* 


0.05* 


25.2*0 




0.1*5 


l*.9*9 i 


. 0.125 


7.6** 


2.778 


-0.059 


1.0 


0.0*1 


25.166 




0.172 


1*.996 • 


. 0.1*7 


2.687 


2.295 


-0.207 


0.7 


0.032 


25.09* 


• 


0,20* 


15.031 J 


. 0.17* 


2.726 


2.327 


-0.3*5 



TMcra 



•«l 
«.» 
«.l 
♦.f 
%*» 
t.« 
».» 
Uf 

ft.A 
S.t 

i»t 
•t.t 

-0.5 
-1.6 
-2.7 



Table 4.1 (Cont'd): Statistics for sulphate weighted trajectories. Source - Sudbury 1978-80. 



SOUBCE STATION IS I 1 - ONT SUUBJBT 
TIHt KEIdHT « BCXAIHINO ME»»( CONC . 



PEMIOOI rsflu 
I ME AH 



HEiG'^Tt src 

TME4N 



SIJL'»H«TE WEIGMIEO 



0. 

3. 

6. 

9. 
12. 
I*. 
2*. 
30. 
36. 
•2. 
*•. 
S4. 
6«. 
M. 
T2. 
TS. 
8*. 
90. 
96. 



1913. ros 
1991. ao« 

2*43.2?] 
2TTT.bB2 
3023. f9T 

3312.602 

33?6.<ik2 

3201. /B9 

2995. 00% 

2?51.I65 

2*71. B]B 

2IT0.7|2 

10116.367 

1611.669 

I36S.072 

1161.996 

991.701 

012.256 

6«9.532 



lOO.O 

131.6 

161.* 

103. 5 

199. S 

210.0 

219.0 

211.5 

197.9 

181. 

163.3 

1*3.* 

12*. 6 

106.5 

90.2 

76.0 

65.5 

55.0 

*5.6 



0.17S 
0.230 
0.29? 
0.]?l 
0.3*9 
0.30* 
0.309 
0.3B1 
0.365 
0.3** 
0.320 
0.292 
0.263 
0.233 
0.206 
0.102 
0.162 
0.1*2 
0.122 



25.790 
25.H** 
?5.900 
25.951 
26.002 
26.090 
26.1*5 
26.177 
26.178 
26.15* 
26.121 
26.00* 
26.013 
25.952 
25.916 
25.853 
25.796 
25.760 
25.752 



0.0 

0.00* 

0.000 

0.010 

0.013 

0.017 

0.022 

0.027 

0.032 

0.030 

0.0** 

0.050 

0.056 

0.063 

O.OTI 

0.079 

0.0S7 

0.095 

0.105 



15.560 
15.529 
15. ^9* 
15.**8 
15.399 
15.292 
15.175 
15.066 
1*.96T 
I*. 091 
I*. 8*2 
I*. BIS 
1*.797 
I*. 769 
1*.7|9 
I*. 71* 

i*.rt2 

l*.7?l 

i*.7oa 



0.0 

0.005 

0.00* 

0.011 

0.01* 

0.019 

0.023 

0.027 

0.031 

0.03* 

0.03B 

0.0*2 

0.0*7 

0.052 

0.057 

0.063 

0.069 

0.076 

0.0B6 



»SIG^«« 


VSI&M« 


COV«ai«'«CC 


0,0 


0.0 


e.e 


0.19? 


0.209 


0.006 


0.37? 


0.*06 


0.02? 


0.5*3 


0.59t 


0.0** 


0.706 


0.768 


0.068 


0.999 


1.082 


0.133 


1.273 


1.3*2 


0.231 


1.532 


1.53* 


0.383 


1.76* 


1.687 


0.SI6 


1.97? 


I.BOT 


0.616 


2.166 


1.896 


0.655 


2.331 


1.963 


0.65B 


2.**0 


2.028 


O.S** 


2.538 


2.070 


0.421 


2.615 


2. II* 


0.261 


2.6«2 


2.138 


0.19? 


2.727 


2. 159 


0.0*6 


2.75* 


2.18* 


-0.00* 


2.r5i 


2.250 


-0.063 



T*«T« 



9. 

3. 

2. 

I, 
0. 

-0. 

-0. 




o 



Table 4.2: Statistics for unweighted and sulphur dioxide weighted trajectories. Source - Baldwin 1978-80. 



SOO«CF STATION ISl 20 - ILL H«LOWIH 
TIME lieiGKT • »EH«|f4lNG 



Ptaiooi TH10 



HEIGHT J srC 

YME4N 



UHMEISHTEO 



TSIGM* C0VAOU-4CC 



TH€T» 



•i* 


BA52. 




3* 


■452. 




i» 


a6!i2. 




«. 


•A«2. 




tz. 


BAse. 




la. 


iS6«. 




2*. 


a«B6. 




30. 


am. 




36. 


8|T2. 




*;. 


anis. 




*8. 


Tf(30. 




s«. 


TfrlT. 




6». 


7*23. 




66. 


Tl9I. 




T?. 


6966. 




IB. 


6700. 




S*. 


6*28. 




«>e. 


6|90. 




96. 


59*7. 





23.*I0 
23. *3* 
23. 452 
23.*63 
23.*7l 
23.*95 
23.528 
23.576 
23.622 

23.<sao 

23.7*6 
23.788 
23.810 

23.857 
23.889 
23.931 
23.976 
2*. 030 
2*. 76 



0.0 

0.002 

0.00* 

0.007 

0.009 

0.012 

0.015 

O.OIS 

0.020 

0.022 

0.02* 

0.026 

0.027 

0.028 

0.029 

0.030 

0.031 

0.0 32 

0.033 



13.110 
13.126 
11.15* 
13.191 
13.237 
I3.36T 
13.525 
13.702 
13.869 
I*. 02* 
1*.I66 
1*.296 
1*.*0« 
I*. 502 
I*. 595 
I*. 891 
]*.77B 
I*. 858 
1*.908 



0.0 
0.001 

o.oei 

0.007 
0.009 
0.013 
0.015 
0.017 
O.I)|9 
0.021 
0.022 
0.023 
0.02* 
0.025 
0.02r> 
0.027 
0.028 
0.028 
0.029 



0.0 
0.205 

o.*n 

n.6l5 
0.B09 
1.139 
1.*2I 
1.653 
1.8*8 
2.006 
2.128 
2.230 
2.323 
2.390 
2.*** 
2.*78 
2.*96 
2.521 
2.535 



0.0 

0.239 

0.*6« 

0.675 

0.8*7 

1.167 

l.*03 

1.59* 

1.7*8 

1.87* 

1.968 

2.031 

2.088 

7.130 

2.I7S 

2.196 

2.211 

2.229 

2.233 



0.0 

0.009 

0.O36 

0.080 

0.136 

0.233 

0.321 

0.3«3 

0.370 

0.328 

0.253 

0.171 

0.061 

0.039 

•0.073 

-O.IBI 

-0.125 

•0.172 

-0.202 



0.0 
12.1 
11.9 
11.9 
II. • 
18.2 

9.a 

T.« 
6.2 

*.6 
3.2 
2.0 

0.7 
0.* 
-0.7 
-0.9 
-1.2 
•1.8 

-l.a 








souoct 


STITION tSl 


20 


- Ill eiLDMiN 


pcmooi 7eao 


HEJOMTI 5rC 


sutwim 


OlOKlOe WEIOHTfS 






TINC 


KEIOMT 


« 


BEMilNiNG 


MEAN COHC. 


IM€««I 


TMEIM 


«SI6>«« 


T5I0M* 


COVMUNCE 


THCtO 






2*061.102 




100.0 


2.7B1 


23.*I0 






13.110 


* 


0.0 


0.0 


0.0 


0.0 


0,8 






20852.037 




86.7 


2.*10 


23.*32 






13.109 




0.002 


0.212 


0.251 


O.OII 


u.s 






iaii6.**a 




75.3 


2.09* 


23.**7 






13.119 




0.00* 


e.*25 


0.*85 


0,0*3 


13.* 






15753.793 




65.5 


1.021 


23.*S7 






13.139 




0.00a 


0.635 


0.699 


0.094 


13.3 




12! 


13702.78* 




56.9 


I.5B* 


23.*63 






13.167 




0.008 


0.836 


0.895 


0.162 


13.0 




18. 


10265.71* 




♦ 2.7 


1.199 


23.*88 






13.275 




0.012 


1.172 


1.190 


0.26* 


10.9 




2*. 


7679. Bo3 




31.9 


0.908 


23.521 






I3.*I6 




0.016 


l.*56 


l.*2B 


0.351 


9.6 




30. 


5732. I3S 




23.B 


0.690 


23.562 






13.577 




0.022 


1.690 


1.836 


0.302 


T.6 


T" 


36. 


*272.898 




i7.a 


0.523 


23.59B 






13.733 




0.028 


1.883 


1.807 


0.350 


5. a 


»K •» 


*2. 


3178.115 




13.2 


0.396 


23.638 






13.87* 




0.03* 


2.037 


1.9*3 


0.2** 


3.* 


*B. 


235*. a** 




9. a 


0.301 


23.676 






l*.002 




0.0*2 


2.160 


2.033 


0.092 


I.I 


~ — ■ 


5*. 


17*2.65* 




7.2 


0.229 


23.683 






1*.I16 




0.050 


l.tltt^ 


2.09* 


-0,039 


-a. 6 




60. 


129l.a*6 




5.* 


0.17* 


23.682 






1*.227 




0.060 


2.3*9 


2.169 


-0.168 


-I.T 


n 


66. 


950.571 




*.o 


0.132 


23.712 






1*.315 




0.072 


7.*I5 


2.227 


-0.137 


-t.3 


8 


72. 


69*. 9*3 




2.9 


0.100 


23.709 






t*.3B2 




0.087 


2.*5* 


2.290 


-0.201 


-2.7 


TB. 


508.9*5 




2.1 


0.076 


23.717 






1*.*33 




0.102 


2.*5* 


2.310 


-0.375 


-3.8 


o 


a*. 


3A9.B68 




1.5 


0.058 


23.731 






1*.*8I 




0.121 


2.»52 


2.322 


-0.*73 


-*.5 


Q. 


90. 


271.260 




I.I 


0.0** 


23.75* 




0.150 


1*.521 




0.1*1 


2.*65 


2.330 


-0.585 


-5.5 


1 

n 

n 
o 

E 
a 

o 


96. 


199.9*8 




O.B 


0.0 3* 


23.793 




0.17* 


1*.536 




0.165 


?.*65 


2.328 


-0.610 


-5.7 



Table 4.2 (Cont'd): Statistics for sulphate weighted trajectories. Source - Baldwin 1978-80. 



OUBCE 


<(T*rlON ISl 


20 


• ILL BALOrilM 


PERIOOl 7880 


rttlGIT 


srr 


SUL»H»TE 


KEIOHTEO 






lint 


•EIGHT 


« 


RCMAINlHO 


HtkH CONC. 


IMEAN 


T><E»«< 


HSIGX* 


YSI&H» 


COVAPIAMCE 


TMCT» 


e« 


1513. TflS 




loo.o 


0.175 


23.OI0 


w 


0.0 


3.110 


* 


0.0 


0.0 


0.0 


n.o 




3, 


22*5.201 




loe.3 


0.260 


23.002 


* 


0.000 


3.128 




0.005 


n.i93 


0.725 


0.010 




6> 


2(170. •9' 




189.7 


0.332 


23.070 


* 


0.007 


3.160 




0.00" 


n.387 


0.029 


0.037 




9. 


3350.3)6 




221.3 


0.387 


23.090 


* 


0.010 


3.19* 




0.011 


0.569 


0.617 


0.079 




!?• 


3698.725 




2i>*.3 


0.028 


23.516 


t 


0.012 


3.230 




0.013 


0.7O7 


0.789 


0.133 




|B> 


*0lia.73« 




270.1 


0.O77 


23.570 


« 


0.016 


3.313 




0.017 


l.OOS 


1.057 


0.226 




2*. 


«l<il.5a) 




273.6 


0.090 


23.632 


• 


0.020 


3.O03 




0.020 


1.301 


1.282 


0.310 




}0> 


*0n).28J 




260.3 


0.091 


23.687 


•> 


0.020 


3.097 




0.023 


1.510 


1.077 


0.3$l 




36. 


37 1 |.}$8 




245.2 


0.O5O 


23.732 


• 


0.078 


3.569 




0.027 


1.685 


1.621 


0.331 




*2> 


33'>5.*6* 




220.3 


0.023 


23.781 


1 


0.032 


3.60 3 




0.030 


1.801 


l.t06 


0.280 




*8. 


3056.0117 




201.9 


0.390 


23.826 


• 


0.036 


3. 736 




0.031 


1.969 


1.806 


0.182 




S*. 


27?2.8«2 




179.9 


0.357 


23.855 


* 


O.OOO 


3.810 




0.037 


2.063 


1.921 


0.133 




60 • 


2393.590 




158.1 


0.322 


23.876 


1 


O.OOO 


3.888 




0.001 


2.137 


2.006 


0.031 




66 a 


20«>>7.8e0 




136.6 


0.288 


23.909 


t 


0.0O9 


3.936 




0.0O6 


2.207 


2. 088 


0.011 




T2» 


1776.726 




IIT.O 


0.255 


23.921 


* 


0.053 


3.901 




0.050 


2.238 


2.128 


-0.125 


-l«* 


TB. 


1508. «3T 




99. T 


0.225 


23.965 


♦ 


0.058 


3.950 




0.055 


2.250 


2.153 


-0.263 


-3.8 


•«• 


1272.901 




80.1 


0.198 


20.007 


Jf 


0.06O 


3.980 




0.061 


2.271 


2.179 


-0.J8* 


-0.3 


90 • 


1087. *96 




71.8 


0.176 


20,062 


# 


0.070 


0.010 




0.06A 


2.297 


2.170 


-fl.028 


-0,8 


96. 


936.755 




61.9 


0.158 


20.126 


* 


0.075 


O.006 


t 


0.071 


2.302 


2.ITI 


-8.J66 


-3.7 




•S9 



Table 4.3: Statistics for unweighted and sulphur dioxide weighted trajectories. Source CONMJO ALLEG 1978-80. 



SOURCE sTiTlON IS! 28 - PEN COWJO »LlEr. 
TIte HEIGHT « RtM«tNI>IG 



9* 


«AS2. 


lOD.O 


3, 


B«.S2. 


lon.o 


ft* 


8».S?. 


100.0 


»* 


8^41. 


<i9.<t 


It. . 


■k|8. 


99.6 


|«. 


847?. 


98.0 


M. 


87B8. 


95. B 


Wf 


Sn6S. 


93.2 


«». 


TTSfl. 


89. T 


M« 


T**l. 


H6.0 


*«> 


7t?*. 


HZ. 3 


f*. 


6T63. 


TH.? 


M. 


«4<i%. 


T*.5 


m» 


6138. 


TO. 9 


P-* 


5838. 


6T.5 


ft. 


Slft-i. 


64.3 


«• 


511*. 


*|.9 


9t. 


5055. 


^k* 


m* 


4II34. 


55.9 



P?-jrOD! THOU 


MtlGHT 


t '.rr 


UVMtlGHTEO 




^" •'.' 


XHf AM 


YH£»s 


PISIGW* 


TSIGH* 


COVADUNCE 


?6.160 


> 0.0 


l].6i0 




O.D 


0.0 


0.0 


n.o 


?6.?ll 


> o.oo? 


11.645 




0.00:1 


0.200 


0.743 


0.008 


26.?5>» 


> 0.004 


13.647 




0.005 


n.i9 7 


0.475 


0.0 3S 


J6.795 


f 0.006 


13.658 




0.007 


0.587 


0.691 


0.085 


^6.]^3 


> O.OOB 


13.676 




0.01ft 


n.7ii9 


0.887 


0.158 


?6.3»8 


0.012 


13.732 




ft. oil 


1.093 


1.710 


0.367 


?6.3<>5 


> 0.015 


13.801 




0.016 


I.37S 


1.472 


0.807 


26.331 . 


0.018 


I3.BB6 




0.0)9 


1.625 


1.681 


0.871 


?6.794 . 


. 0.021 


13.975 




0.021 


1.03? 


l.«4| 


1.126 


?6.?32 . 


> 0.023 


14.041 




0.023 


2.005 


1.964 


1.353 


;6.|6H 


, 0.025 


14.105 




0.024 


2.145 


2.054 


1.526 


?6.ftBl < 


» 0.028 


14.155 




ft. 025 


2.263 


2.095 


1.622 


?6.n2l < 


0.030 


14.210 




0.026 


2.380 


2.124 


1.680 


?5.94S ; 


0.03] 


14.273 




0.027 


2.468 


2.153 


1.694 


25.1155 . 


0.013 


14.324 




0.028 


2.542 


2.175 


1.591 


25.787 < 


' 0.035 


14.393 




0.829 


2.612 


2.200 


1.516 


25.730 . 


0.037 


14.442 




0.030 


2.669 


2.219 


1.417 


25.669 ; 


0.038 


14.480 




O.03I 


2.700 


2.228 


1.440 


25.653 ! 


- 0.039 


14.547 


• 


n.03? 


2.737 


2.744 


1.397 



THtr* 



11. « 
tt,* 

!••• 
Ift.fl 
!*.« 
lt.3 
19*5 
!«.« 
It.* 
t%* 
l».9 
t«.S 

tt.« 

lltS 

J«.6 







I 

I 

i 



SOURCE 


STtrlOM ISI 


28 


- PCH CONMJO »LLEG 


PERIOOI 7880 


MEtGHTi src 


SULPMUR 


oiovior HEIGMTCO 




TIHC 


WEtGHT 


% 


REMAINING 


ME*"* COHC. 


■ MEAN 


TME4N 


il51G<4« 


r5lGM« 


COVtRUttCE 


THCT8 


'»•' 


24061.102 




lot.o 


2.781 


76.160 


• 


0.0 


13.650 




0.0 


0.0 


o.n 


0.0 


9.9 


"!•: 


20677.166 




85.9 


2.390 


26,215 


* 


0,001 


13.616 




0.007 


0.209 


n.754 


0.009 


11 »f 


«» 


17887.030 




74.3 


2.067 


26.763 


• 


0.003 


13.626 




0.004 


ft. 412 


0.496 


n.036 


tz.t 


%t 


15495.211 




64.4 


1.793 


26.100 


♦ 


0.003 


13.622 




0.006 


ft. 607 


0.719 


0.088 


t3.« 


If* 


134)2.538 




55.7 


1.556 


26.328 


* 


o.oor 


13.624 




0.008 


0.793 


ft. 920 


n.166 


!*.• 


lit 


9947.002 




41.3 


1.173 


26.344 


• 


0.011 


13.646 




0.012 


1.125 


1.746 


0.3S0 


i».r 


«%* 


732"'. 006 




30.5 


0.884 


26.316 


I 


0.017 


13.671 




0.018 


1.414 


1.502 


0.627 


1T.6 


H* 


5368.944 




22.3 


0.666 


26.266 


t 


0.023 


11.703 




0.021 


1 . 668 


1.703 


n.B83 


IT.A 


»* 


3876.109 




16.1 


0.500 


26.195 




0.030 


13.738 




0.030 


1.874 


1.837 


l.lil 


tt.« 


mw; 


2787.102 




11.6 


0.375 


26.097 




0.019 


13.741 




0.037 


2.041 


1.946 


1.315 


IT.S 


«•* 


2007,529 




8.3 


0.282 


26.004 




0.049 


13.737 




0.045 


2.189 


2.033 


1.492 


ina 


■s*-» 


1411.890 




6.0 


0.212 


25.880 




0.061 


13.712 




0.054 


2.297 


2.054 


1.526 


\*.t 
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Table 4.3 (Cont'd): Statistics for sulphate weighted trajectories. Source CONMJO ALLEG. 1978-80. 
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CMC X CO-ORD 

Figure 4.1: Mean trajectories for Sudbury, Baldwin, and CONMJO ALLEG. 
January 1978. Note that these trajectories do not neces- 
sarily represent the path followed by any puff of pollutant 
and cannot be used for modelling long range transport 
without a knowledge of the dispersion about them. 
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CMC X CO-ORD 

Mean trajectories for Sudbury, Baldwin, and CONMJO ALLEG. 
July 1978. Note that these trajectories do not necessarily 
represent the path followed by any puff of pollutant and 
cannot be used for modelling long range transport without a 
knowledge of the dispersion about them. 
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CMC X CO-ORD 

Figure 4.3: Mean trajectories for Sudbury, Baldwin, and CONMJO ALLEG. 
January 1978/79/80. Note that these trajectories do not 
necessarily represent the path followed by any puff of pol- 
lutant and cannot be used for modelling long range transport 
without a knowledge of the dispersion about them. 
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CMC X CO-ORD 

Figure 4.4: Mean trajectories for Sudbury, Baldwin, and CONMJO ALLEG. 
July 1978/79/80. Note that these trajectories do not neces- 
sarily represent the path followed by any puff of pollutant 
and cannot be used for modelling long range transport 
without a knowledge of the dispersion about them. 
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CnC X CQ-ORD 

Figure 4.5: Mean trajectories generated form January 1978 to December 
1980 for each source. Sources 1-15. Note that these 
trajectories do not necessarily represent the path followed 
by any puff of pollutant and cannot be used for modelling 
long range transport without a knowledge of the dispersion 
about them. 
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Figure 4.6: Mean trajectories generated from January 1978 to December 
1980 for each source. Sources 16 - 30. Note that these 
trajectories do not necessarily represent the path followed 
by any puff of pollutant and cannot be used for modelling 
long range transport without a knowledge of the dispersion 
about them. 
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CMC X CO-ORD 

Figure 4.7: Mean trajectories generated from January 1978 to December 
1980 for each source. Sources 31 - 45. Note that these 
trajectories do not necessarily represent the path followed 
by any puff of pollutant and cannot be used for modelling 
long range transport without a knowledge of the dispersion 
about them. 
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CMC X CO-ORD 

Figure 4.8: Mean trajectories generated from January 1978 to December 
1980 for each source. Sources 46 - 60. Note that these 
trajectories do not necessarily represent the path followed 
by any puff of pollutant and cannot be used for modelling 
long range transport without a knowledge of the dispersion 
about them. 
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Figure 4.9: Plot of time evolution versus standard deviation of the x- 
and y- coordinates. Source - Sudbury. 
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Figure 4.10: 



Plot of time evolution versus standard deviation of the x- 
and y- coordinates. Source - Baldwin. 
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Figure 4.11: Plot of time evolution versus standard deviation of the x- 
and y- coordinates. Source - CONMJO ALLEG. 
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5, SUMMARY AND RECOMMENDATIONS 

The most significant findings of the present study are those 
given in Figures 4.5 - 4.11, which show the means of the trajectories 
initiated from all sixty sources every three hours for a period of three 
years, and the standard deviations about these means for selected 
sources. The means of the trajectory-end-points do not go in a west to 
east direction which is rather expected for trajectories that are based 
on surface winds. Since the travel of pollutant parcels are expected to 
be also affected somewhat by upper level winds, a comparison of January 
and July 1982 surface and 850mb trajectories were undertaken and 
reported in the ADDENDUM to this report. These comparisons between the 
two different level trajectory statistics pointed to some difficulties 
due to small sample size and possibly due to domain boundary effects to 
a significant degree. It was found however that 850mb trajectories, in 
general, move farther away than the surface trajectories and in a more 
intuitively expected direction. Also, winter time trajectories tend to 
behave in a more unexpected fashion than summer time trajectories. An 
examination of the percentage loss of trajectories as a function of time 
from trajectory initiation strongly suggest the importance of trajectory 
loss in estimating trajectory statistics. 

The results point out to several potential difficulties in the 
Lagrangian and the statistical model formulations. 
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1. The first area of concern is the wind field formulation in the 
Lagrangian model, especially in the calculation of the surface 
winds from surface pressure data, and perhaps more importantly, the 
choice of a single level (e.g., surface versus 850mb) for which the 
calculated winds would be representative of pollutant transport. 
The following suggestions could be made with the aim of improving 
the model assumptions: 

i) Comparison of trajectories with available tetroon-trajectory 
data, (Hoecker, 1977. Pack et al., 1978; Warner, 1981; 
Reisinger and Mueller, 1983), and tracer experiments (e.g., 
CAPTEX), 

ii) Comparison of trajectories with those generated by other avail- 
able trajectory models, 

iii) Continued comparison of Lagrangian Model results (including the 
mass budget calculations) to available air quality and 
deposition data or predictions from Eulerian Models (this 
comparison however is not very useful for evaluating the 
trajectory module only), 

iv) Use of a combination of the trajectories for the two levels for 
pollutant transport modelling. 
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2. Another area of concern is that pointed out by the 850mb-surface 
comparison study and also by an examination of the percentage of 
remaining trajectories, namely boundary effects and early 
termination of trajectories. It is suggested that: 

i) Reasons for trajectory termination be studied quantitatively, 
i.e., determine if a trajectory is terminated due to edge of 
calculation domain, gaps in data, or non-convergence due to 
data sparsity near oceans, and produce statistics on the occur- 
rence of these reasons (this would indicate which source of 
inaccuracy would be tackled first), 

ii) Related to (i) above, if it is found that domain boundary ef- 
fects are significant, trajectories for a subset of the data be 
calculated using a larger calculation domain than presently 
used and their statistics compared to the corresponding statis- 
tics calculated in the present study, 

lii) A larger database be used for the 850mb and surface compari- 
sons, 

iv) The trajectory statistics for "percentage remaining" less that 
- 95 % not be used until (i) above is completed. 
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The aim of the present study was to generate the statistical 
parameters for use in the OME Statistical Model. To further pursue this 
goal, the following recommendations are made: 

3. Examination of the effect of the coordinate system (i.e., the CMC 
grid) on the statistics generated, by comparing them to those cal- 
culated using the "tangent plane" approach described in the interim 
report (Alp and Moran, 1983; note that the latter approach has not 
yet been fully tested), 

4. Test for normality of the distribution of trajectory end-points to 
determine the suitability of the Gaussian distribution used in the 
Statistical Model, 

5. The present results indicate that using x = ut and y = in the 
Statistical Model (see the interim report, Alp and Moran, 1983) 
may be inappropriate. For this reason it is recommended that the 
statistics generated in the present study be implemented in the 
Statistical Model and its predictions be compared to available air 
quality and deposition data. 

6. The present results indicate that at least for certain sources, the 
mean trajectories move very little and the diffusion ellipses are 
close to circular. Thus, in order to quantify the differences in 
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various directions more accurately, a splitting of the statistical 
model calculation domain into sectors may be more appropriate. 
This may be accomplished by using a one-dimensional distribution 
for each sector instead of the presently used two-dimensional dis- 
tribution function. Parameters for such a distribution may be 
determined from the presently used trajectory end-point data, but 
using perhaps f. to indicated mean trajectory movement in each 
sector i from a given source point. It is suggested that consider- 
ation be given to a separate research effort to develop such a 
statistical model, as it appears that this may provide a 
significant enhancement over the present model. 
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APPENDIX A 



Further Results on Trajectory Statistics 
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APPENDIX A 
Further Results on Trajectory Statistics 

Appendix A contains statistics for unweighted, sulphur dioxide 
weighted, and sulphate weighted trajectories for each source (Sudbury, 
Baldwin, and CONMOO ALLEG near Pittsburgh, Pennsylvania) for the 
following periods: 1980, 1979, 1978, January, April, July, October, 
January 1978, July 1978. The table entries are described in Section 4 
of the main text. 
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0.028 


0.8*3 


0.050 


1.165 


0.1*7 


l.**2 


0.290 


1,657 


0.*1T 


I,«35 


0.5*9 


I. "67 


0.680 


2.044 


0.799 


2,071 


0.786 


2,113 


0.659 


2,119 


0.368 


2.103 


0.1*0 


2.130 


-0.12* 


2.188 


-0,2*0 


2.7?0 


-n.isi 


2.268 


-0.*09 



THET» 



0.0 

5.8 

5.0 

4.7 

4.9 

7.4 

9.1 

9.5 

9. A 

10.1 

10.0 

8.9 

T.7 

3.9 

I.* 

-I.? 

-2.3 

-3.3 

-3.6 



SOURCE ST»TI0>< IS! 1 - OM SUOBJOT 

T|HE HEIGHT « I'EMAI'JiNd ME»»I COMC. 



"E'lIODl I"»90 

<MEftN 



HEir.HT! <;rc 
y»<E*') 



0. 

3. 

6. 

9. 
IE. 
IB. 

a*. 

30. 
36. 
»2. 

48. 
54. 
60. 
66. 
TZ. 

Ta, 

8*. 
90. 
96. 



5n6.<i05 

64T.363 

TB5.31Z 

8«9.593 

9AA.I6T 

10ST.J?8 

10TI.«l^5 

10*1. ?9'> 

9T3.42* 

911?, 3T6 

an6.9*e 

69T.7?3 
598. 5?0 
5|l.0i>3 
1(33.968 
36*.<i39 
3n8.076 
2'i8.*9<< 
?l3.9ft? 



loo.o 

127.8 

155.1 

175.7 

191.? 

?08.l« 

311.6 

205.6 

19?.? 

17t).? 

159. i 

137. H 

118.? 

100.9 

85.7 

72.0 

60.8 

51.0 

*?.? 



0.175 
O.??" 
0.271 
0.307 
0.335 
0.367 
0.375 
0.373 
0.35 7 
0.3<t0 
0.316 
0.287 
0.256 
0.228 
0.201 
0.176 
0.1 56 
0.138 
0.J19 



25.790 
25.BI.7 
25.<t99 
25.9fc<. 
25.9fl"» 
?6.053 
26.090 
26.099 
26.068 
26.008 
25.915 
25.769 
25.«.13 
25.533 
25.<.?7 
25.321 
25.19* 
25. |»? 
25.153 



0.0 

0.007 

0.013 

O.OM 

0.023 

0.031 

0.039 

0.0<i7 

0.055 

0.06<> 

0.075 

o.on7 

0.096 
0.107 
0.120 
0.133 
O.I<.& 
0.159 
0.178 



15.560 
15.516 
15. '.68 
15.1.05 
15.339 
15.207 
l5.0Bi> 
Id. 95* 
l«.861 
Ik.SOO 
I*. 7*0 
I*. 656 
1<..63? 
1*-59| 
1».5?? 
1'..537 
l*.5e5 
14.650 
!<..597 



0.0 
0.00« 

0.01* 

0.020 
0.02* 
0.03? 
0.0*n 
0.0*6 
0.05'' 
0.060 
0.067 
0.073 
0.081 
0.087 
0.09* 
0.106 
0.119 
0.13* 
0.150 



«;dL''h»TE wf ir.HTED 



n.o 

0.186 
0.365 
0.5*0 
0.71 I 
1.006 
1.270 
1.507 
1.720 
I .910 
?.l** 
?.?9| 
?.3*l 
2.*1* 
?.*93 
2.53? 
?.555 
2.551 
2.601 



tv.\r,<*t 


COViPIfiCE 


o.n 


0.0 


0.70* 


0.006 


0.*06 


0.0?? 


0.588 


0.0** 


0.759 


0.07* 


1.051 


0.159 


1.320 


O.?N0 


l.*96 


0.*I8 


1.657 


0.551 


1.798 


0.669 


1.^92 


n.T?3 


1.918 


0.610 


1.972 


0.*58 


1.965 


n.itu 


I. '52 


O.IJ* 


2.027 


-0.131 


2.097 


-0.1T7 


?.I*7 


-0.18* 


?.?00 


-n.093 
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0.0 

".? 

9.3 

".5 

8.T 

8.9 

9.9 

10.* 

in. 6 

10.7 

8.9 

6.6 

*.8 

2.9 

1.1 

-1.7 

-1.6 

-1.6 

-0.8 
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SOU«Cf ST*TI0»J IS! I - ONI SUOHU^r 
Tl»tE WEIGMr » OEXtlNiNG 



0« 


738. 


100.0 


3* 


738. 


100.0 


6< 


738. 


loo.o 


9. 


738, 


loo-o 


12< 


738, 


100.0 


IS. 


733. 


99.3 


?*. 


706. 


95.7 


3P. 


^6*. 


90.0 


3<b. 


0.3?. 


85.6 


*2. 


59b. 


80.6 


«8. 


SS8. 


TS,* 


5*. 


•>19. 


14.) 


60. 


aBZ. 


65.3 


66. 


(>62. 


6e.6 


72, 


(■<>6. 


60.* 


78. 


429. 


5B.1 


8fc. 


<.08. 


55.3 


90. 


385. 


5?. 2 


96. 


359. 


48.6 



PEOIOO: J4S 


HEIGH 


! src 


i>Mi«EI'.HTEn 








X««£*«l 


Y-<t;»M 


««;ir.ii* 


Y>^TfiM» 


cnw*ni»NCE 


ThEt* 


?5.790 


0.0 


15,560 


t 0.0 


n.O 


0.0 


0.0 


8,0 


25.817 . 


O.OIO 


15.472 


> o.on 


0.26? 


0.287 


-n.005 


-4.4 


?5.«60 < 


0.018 


15.377 


0.02n 


O.502 


0.553 


-0.028 


-6.? 


25.918 . 


0.027 


15.27? 


I 0.030 


0,729 


0.802 


-n.071 


-7.6 


25.985 < 


► 0.035 


15.157 


L 0.038 


0.943 


1.038 


-0.129 


-8.2 


26.105 < 


0.0«8 


l*.913 


I 0.054 


1.306 


1.465 


-0.206 


-6,9 


26.163 < 


' 0.061 


1'..759 


t 0.067 


1.616 


1.T82 


-0.103 


-2.3 


26.186 < 


> 0.073 


1*.705 


► 0.077 


1.890 


1.978 


-0.051 


-0.8 


26.158 . 


• 0.0811 


10.568 


. 0.0^6 


2.124 


2.156 


0.075 


1.0 


26.033 < 


> O.0Q4 


|(>.b3? 


t 0.097 


2.289 


2.355 


0.37? 


4,1 


25.839 < 


U.I02 


14.309 


I 0.106 


2.411 


2.512 


0.6S4 


«.6 


25.61J < 


' o.ni 


14.203 


. 0.111 


2.535 


2.527 


0.708 


6.3 


25.365 . 


. 0.123 


14.205 


I 0.114 


?.693 


2.511 


0.641 


9.0 


25.155 < 


0.129 


14.259 


I 0.117 


2.78? 


2.518 


8.377 


*.« 


2<>.9ft8 < 


0.135 


14.330 


t 0.1l« 


2.849 


2.484 


0.097 


n.1 


24.733 < 


0.139 


14.458 


. 0.119 


2.883 


2.475 


-0.231 


-1.6 


2*. 600 < 


0.|4I> 


14,616 


L 0.124 


2.>»06 


2.505 


-O.SI" 


-1.5 


2'..'.90 . 


0.I49 


14.775 


L 0.131 


?.923 


2.567 


-1.0*9 


-7.0 


24. (.IS < 


0.153 


14,797 


L 0.135 


?.9no 


2.566 


-1.540 


-10.4 



SOUBCF 5T*TlO»J IS! I - ONT SUOBOBr 

TME KtlGHT t REM/IIMING ME 



C' 



n 
o 

■3 

a 
o 

a 

V 

o 



3> 
«• 

9. 
12. 
18. 

24. 
30. 
36. 
42. 
48. 
54. 
60. 
66. 
72. 
T8. 
84. 
90. 
96. 



3019.960 

2624.983 

2273.880 

1975.25'. 

17?0,599 

1299.393 

953.171 

690,161 

503.455 

361.367 

262.953 

187.857 

133.041 

96.584 

70,732 

51,588 

37,443 

?7,002 

19.369 



100.0 

86.3 

74.8 

65.0 

56.6 

42.7 

31.4 

22.7 

16.6 

12.0 

8.6 

6.2 

4.4 

3.2 

2.3 

1.7 

1.2 

0.9 

0.6 





PERIOD! J«M 


HEIGHT 


I SfC 


»9 COHC. 


X««E*N 


rMttN 


4,119 


25.790 . 


0,0 


15.560 




u.o 


3,557 


25.819 


> 0.005 


15.4f2 




0.006 


3.081 


25.865 


> 0,010 


15.380 




O.on 


2,676 


25.925 


' 0.016 


15.280 




0.017 


2.331 


25.995 


► 0.022 


15.170 




0,024 


1.773 


26.120 


• 0.035 


14.936 




0.039 


1.350 


26,178 


t 0,051 


14,798 




0.055 


1.0 39 


26,186 . 


0.071 


14.733 




0.072 


0.797 


26.150 ! 


> 0.093 


14.589 




0.092 


0.611 


26.023 i 


0.119 


14.457 




0.118 


0.471 


25,866 < 


' 0.149 


14.349 




0.14Q 


0.362 


25.650 < 


► 0.186 


14.236 




0.1 78 


0.2 76 


25,407 < 


0.235 


14.240 




0.211 


0,209 


25,192 < 


' o.2e> 


14,287 




0.249 


0,159 


?4,982 


. 0.34 1 


14,343 




0,287 


0,120 


24,730 


> 0.399 


14,4 35 




0.135 


0,092 


24,596 . 


0.471 


14,537 




0.39'; 


0,070 


24. SOB < 


0.555 


14.656 




o.4n 


0,054 


24,447 


0.553 


14.695 




0.55«> 



SULPHUR OIO«IOE WEIGMTED 
«SIG>«* rSIGM* COV»BI«NCE 



0.0 

0.260 

0.495 

0.717 

0.925 

1.2'9 

1.581 

1.856 

2.0"8 

2.259 

2.419 

2,552 

?,7in 

2,799 

?,865 

2,864 

2.879 

2.885 

2.873 



0.0 


0,0 


0.283 


-0.006 


0.539 


-0.029 


0.774 


-0.068 


0.997 


-0.1 J9 


1.401 


-0.193 


1.684 


-0.12* 


l."85 


-0.014 


2.054 


0.063 


2.?SI 


0.355 


2.418 


0.711 


2.443 


0.731 


2.436 


0.626 


2.4*7 


0.296 


2.416 


-0.029 


2.408 


-0.305 


2.420 


-0.586 


2,455 


-1.169 


2,455 


-1.746 



THET» 



0.0 
-4.9 
-6.7 
-7.6 
-7.9 
-6.7 
-2.8 
-0.6 

o.» 
*.o 
»,* 
&<« 
*.♦ 

2.1 

-0.? 
-2.1 

-4.0 

-8,0 

-11,9 
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SOURCF -STATIOH ISS 1 - ONT SUOflJR* 
TIME WEIGHT * aCM^IMINf, Mt*s CONC. 



«HF*M 



HEir.HTX SrC 



0* 

«* 

«. 

12. 
1«. 

3«. 
M. 
*«. 
•«« 
5«* 

M,. 
72. 

78. 

at*. 

9«. 



176. nb 
2n7.8n5 

^1*.^8l 

25*. 856 
279.989 
2ni.07J 
PA8.259 
;b8.S96 
2?3.r95 

eoi.9o* 

179.150 

155.907 

138.537 

1??.665 

In*. 19* 

90.4*9 

78.ao9 

A9.U0 



loo.o 

I?*. 9 
1*7.* 
166.1 

iao.7 

198.5 

199.3 

190.2 

176.3 

158.7 

1*3.2 

127.0 

110.5 

98.2 

87.0 

73.9 

6*.l 

55.5 

*9,0 



0.191 
0.239 
0.282 
0.J17 
0.3*5 
0.382 
0.398 
0.*0* 
0.393 
0.376 
0.362 
0.3*5 
0.323 
0.300 
0.275 
0,2*3 
0.222 
0.203 
0.193 



25.790 
2'i.fl25 
25.^72 
25.935 
26.01* 
26.168 
26.257 
26. 2H* 
26.26* 
26.139 
25.979 
25.775 
25.579 
25.389 
25.221 
2*. 950 
2*.fl58 
2*. "95 
2*. 908 



0.0 

0.019 

0.033 

0.0** 

0.055 

0.073 

0.091 

0.111 

0.130 

0.1*9 

0.170 

0.19* 

0.?2? 

0.2*2 

0.263 

0.282 

0.302 

0.329 

0.355 



IS. 560 
I5.*73 
15.379 
15.279 
15.171 
I*. 9*2 
1*.795 
1*.69* 

i*.5oe 

1*.367 
14.271 
1*.183 
1*.206 
1*.237 
1*.266 
1*.28* 
I*. 31* 
1*.359 
1*.*21 



0.0 

0.020 

0.03* 

O.OiS 

n.055 

D.0 76 

0.089 

0.105 

0.120 

0.1*0 

0.159 

O.I 72 

0.186 

0.199 

n.2ll 

0.22* 

0.2J8 

0.253 

fl.?7T 



5JLPH»Tr ^EIOHTED 



n.o 

0.250 

n.*T? 

0.680 
0.878 
I .230 
1.532 
I. 812 
?.0*T 
2.235 
2.*19 
2.59IT 
2.769 
7.950 
7.913 
7.876 
2.969 
2.910 
2.9*9 



*5ir,Mft 


COVARIiNCE 


0.0 


0.0 


0.763 


-0.007 


n.*8* 


-0.025 


n.68<i 


-0.0*9 


n.imi 


-0.08* 


1.7*1 


-0.1*1 


l.*89 


-0.0*8 


1.713 


0.120 


1.888 


0.22? 


2. "98 


0.*83 


2.763 


0.733 


2.301 


0.719 


?,32« 


0.»85 


2.3*5 


-0.029 


2.33* 


-0.367 


2,289 


-0.593 


2.761 


-0.862 


2.738 


-1.*I3 


2.267 


-1.705 



T-<f T* 



0.0 

-*•.? 

-6,3 

-6.1) 

-6.2 

-S.3 

-1.2 

2.1 

■3.0 

5.5 

7.1 

6.1 

3.6 

-ft. 2 

-2.5 

-*.l 

-6.0 

-9.5 

-II. I 
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SOUBCf <;T»TtON IS! I - ONT SUDHUBY 

TI"E WtlOHT « PEM4INING 



0> 


720. 


100.0 


3. 


720. 


loo.o 


6. 


720. 


100.0 


9, 


720. 


loo.o 


t?. 


720. 


100.0 


18. 


716. 


")<».* 


?4. 


711. 


<9H.6 


30. 


rs98. 


<S6.9 


36. 


68b. 


<J5.0 


*?. 


673. 


■JS.S 


48. 


66/. 


92.6 


5*. 


652. 


90.6 


60. 


6*2. 


89.2 


66. 


f,^9. 


87.* 


fi. 


A06. 


«<..2 


sk. 


■i'Ji,. 


82.5 


»** 


S76. 


80.0 


96. 


S58. 


77.5 


96. 


■;28. 


73.3 



PE»100: i^"* 


HEIGHT; «;fr 


UMWEISHTET 








Ii<t»M 


TffciV 


»5IGH» 
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cnviRtsNCE 


TH^T» 


25.790 


> 0.0 


15.560 


L 0.0 


n.o 


n.o 


n.o 


O.n 


25.722 


0.007 


15.517 


. 0.010 


0, 194 


0.257 


-n.006 


-9.4 


25.660 f 


0.015 


15.473 


O.OIP" 


0.190 


0.490 


-o.oin 


-11." 


2S.606 < 


0.021 


15.429 . 


. 0.026 


n.576 


0.696 


-n.nfct 


-11.0 


25.556 . 


0.028 


15.382 < 


. 0.031 


0.751 


0.886 


.0.09B 


-9.9 


25.<.'.B < 


0.038 


15.293 


0.046 


1.020 


1.233 


-0.122 


-6.7 


25.332 < 


0.047 


15.218 


. O.O^B 


1.2*5 


1.5*7 


-0.125 


-4.6 


25.191. . 


. 0.054 


15.137 - 


0.05'- 


1.421 


1.7*1 


-0.099 


-2.8 


25.052 . 


0.059 


15.052 < 


. o.ori 


1.553 


1.857 


-0.116 


-1.2 


2*. 936 . 


0.065 


14.999 < 


, 0.0 7S 


1.694 


1.9*2 


-0.168 


-1.4 


2*,B60 < 


. 0.072 


14.971 . 


0.0 79 


1.858 


2.037 


-0.226 


-3.T 


2*. 723 . 


0.075 


14.939 . 


. 0.081 


1.923 


2. 127 


-0.191 


-6.0 


2*. 635 . 


0.079 


14.923 ' 


0.08" 


2.003 


2.227 


-0.486 


-6.9 


24.580 < 


. 0.084 


14.907 ( 


0.091 


?.099 


2.287 


-0.579 


-7.5 


2*.*9| . 


0.089 


14.952 


. 0.092 


?.|8n 


2.255 


-0.813 


-9.7 


2<..'.1'> < 


0.091 


14.926 . 


. 0.0 9* 


?.?70 


2.292 


-1.072 


-11.7 


24.289 . 


► 0.096 


15.04* 


. 0.09« 


2.10? 


2,159 


-1.408 


-14.9 


24.187 . 


► 0.100 


15.239 


. 0.101 


2.160 


2.187 


-1.610 


-16.1 


24.079 . 


0.104 


15.320 


0.104 


2.401 


2.179 


-1.847 


-17.8 



SOURCE 5T»TI0»J IS! 1 - ONT SUOSaOY 

TI^E WEIGHT % I>EM»1»JIN6 «£*•< CONC. 
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0. 
3. 
6. 
9. 
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30. 
36. 
*2. 
48. 
5*. 

««. 

n» 

m* 

96. 



1814.524 
1528.337 
1295.309 
Il|0.8h5 

9M.5T9 

725.553 

549.881 

414.961 

315.874 

219. 7*7 

184.266 

140.859 

I0 7.9S4 

83.0*0 

S2.480 

47.845 

16.1S9 

27. I?* 

19,7|4 



100.0 

83.3 

70.6 

60.6 

52.* 

39.5 

30.0 

22.6 

17.2 

13.1 

10.0 

7, 7 

5.9 

4.5 

3.4 

2.6 

2.0 

1.5 

1.1 





PE910D: «PM 


HEIGHT! S^C 


t,H QOHC, 


«ME«N 


rME»N 


2.548 


25,790 




0.0 


15.560 


I 0.0 


2.123 


25.724 




0.005 


15,510 


I 0.007 


1.799 


25.669 




0.010 


15.452 


t 0.013 


1.543 


25.624 




0.016 


15.392 


L 0.02! 


1.336 


25.582 




0.023 


15.337 


'_ 0.028 


1.013 


25.484 




0.035 


15.251 


L 0.0*5 


0. 773 


25.164 




0.050 


15.179 


t 0.064 


0.595 


25.197 




0,066 


15.103 


L 0.081 


0.462 


25.030 




0.085 


15.029 


I n,09" 


0.356 


2*.8i.S 




0.106 


1I..97* 


0.116 


0.276 


2*. 690 




0.132 


14.944 


0.1 J« 


0.216 


2*. 509 




0.159 


14.929 < 


0.167 


0.166 


2*. J 75 




0.190 


14.910 < 


. 0.201 


0.132 


24.272 




0.229 


14. 9J* < 


0.24) 


0.103 


24.137 




0.274 


14.919 < 


. 0.280 


0.081 


24.(128 




0.125 


15.029 . 


0.31* 


0.063 


23.858 




0.375 


15.1*6 < 


> O.*02 


0.0*9 


23. '10 




0.4*2 


15.355 < 


(),*7« 


0.037 


2J.57? 




0.531 


1-.421 ; 


1 0.5S1 



SUL^HJB 0IO«inF WEIGHTED 
ISIGX* YSlGM* COV»RI»MCE 



0.0 
n.|9i 

0.375 
0.5*4 
0,705 
0.955 
1.175 
1.348 
1.51? 
1.6*5 
1 .791 
1.882 
1.97S 
?.085 
2.161. 
2.251 
2.254 
;'.3n2 
2.357 



O.n 


0.0 


0.256 


-0.005 


0.486 


-0.021 


0.689 


-0.050 


0.87H 


-O.0B2 


1.210 


-0.124 


1.490 


-0,|71 


1.6*3 


-0.181 


1.735 


-0.211 


1.790 


-0.27* 


1.871 


-0.401 


1.976 


-0.571 


2.085 


-0.709 


2.192 


-0.87S 


2.215 


-1.111 


2.313 


-1.458 


2.420 


-1.921 


2.487 


-2.209 


2.498 


-'.68* 



T4ET* 



ft.O 

-8.1 

-9.1 

-9.6 

-9.* 

-7.8 

-7.1 

-5.7 

-5.3 

-5.8 

-7.1 

-9.2 

-10.3 

-11.4 

-1 1.6 

-16.0 

-20.7 

-?2.6 

-?5.« 



oo»ct 


St«TlON IS! 


1 - ONT SUDHUBV 


oEuiOOi »3t» 


HEir.H 


! SrC 


SULOMlTE 


^rtGMTEi) 






TFHt 


Kf IGHT 


« REMAINING 


ME»M CUNC. 


IMf «*l 


r*(t»«j 


»5IG'<» 


Y<;|GM* 


COVABU-^CE 


T^ET* 


t« 


li*.658 


lOO.D 


0.159 


25.790 ' 


> 0.0 


15.560 


. 0.0 


0.0 


o.n 


0.0 


n.o 


3* 


1?9.307 


112.8 


0.180 


25.729 


• 0.016 


15.<.9« 


L 0.022 


0,178 


0.249 


-0.002 


-1.* 


«. 


1S7.|<.* 


137.1 


0.218 


25.679 . 


> 0.027 


15.429 


L 0.037 


n.339 


0.469 


-0.006 


-2.8 


9. 


1B2.»07 


159.4 


0.254 


25.635 


► 0.036 


15.365 


. O.Ofc* 


n.486 


0.666 


-O.OIS 


-3.5 


1?. 


^n^.6^o 


176.7 


0.281 


25.590 


> 0.044 


15.313 


L 0.060 


n.5?M 


0.850 


-0.040 


-5.8 


IB. 


27<>.640 


195.9 


0.314 


25.49f, . 


► 0.057 


15.187 


_ n.07T 


0.858 


1.156 


-0.090 


-7.0 


2*. 


213.797 


203.9 


0.329 


25,168 < 


► 0.070 


15.079 


. 0.093 


1.063 


1.418 


-0.1B6 


-9.3 


30. 


2?9.'»03 


200.5 


0.329 


25.196 . 


> 0.082 


14.957 


t n.ioi 


1.2*8 


1.536 


-0.270 


-9.8 


16. 


2?2.9o3 


19(1. <• 


0.326 


25.002 . 


> 0.096 


14.865 


L 0.107 


1.437 


1.604 


-0.396 


♦10.8 


*2. 


2n7.Tl9 


1B1.2 


0.309 


24.749 . 


> O.108 


14.834 


. 0.1l«> 


1.559 


1.676 


-1.533 


-12.4 


«8. 


193.7*5 


169.0 


0.290 


24.547 , 


- 0.124 


14.890 


L n.129 


1.726 


1.797 


-0.769 


-14,5 


5<>. 


177.250 


154.6 


0.272 


24.374 , 


► 0.142 


14.958 


L 0.1<.6 


1.890 


1.945 


-0.939 


-14.7 


60> 


1"59.4»B 


139.1 


0.248 


24.274 1 


r 0.160 


14.928 


L 0.16? 


2.020 


2.048 


-1.202 


-16.4 


66. 


139. 31J 


121. 5 


0.221 


24.230 < 


>■ 0.183 


14.B49 


. O.IBl 


2.160 


2.137 


-1.426 


-17.0 


72. 


115.292 


100.6 


0.190 


24.155 . 


' 0.204 


14.726 


. 0.201 


2.208 


2.157 


-1.798 


->0.2 


78. 


10 1.692 


88.7 


0.171 


24.103 < 


> 0.226 


14.830 


. 0.226 


?.276 


2.?76 


-2.445 


-?5.3 


S<*. 


86.261 


75.3 


0.150 


23.937 < 


► 0.24? 


14.848 


. 0.257 


2.247 


2.188 


-3.351 


-33,6 


90. 


TO-itflB 


61.5 


. 1 26 


23.B34 . 


0.26 7 


15.009 . 


. 0.297 


2.238 


2.490 


-1.786 


-17.1 


96. 


■59.895 


52.2 


0.113 


23.936 . 


0.310 


15.088 


. 0.145 


2,402 


2.667 


-4.270 


-16.5 



€• 



SOU»Cf STATION ISt 1 - ONT SUO'^UPr 
If^t WEIGHT » PEHAI>4ING 



9. 


T**. 


100.0 


'%■ 


TUt, 


lOO.O 


6. 


7«.<.. 


100.0 


9. 


T**. 


lOO.Q 


13. 


T4<.. 


lOO.O 


IB. 


T**. 


loo.o 


?*. 


7*1.. 


loo.o 


30. 


7*4. 


100.0 


36. 


7*1. 


99.6 


*z. 


7 3B, 


99.? 


*«. 


7?9. 


98.0 


5*. 


717. 


96.* 


60. 


707. 


95.0 


66. 


697. 


93.7 


71. 


6T3. 


90.5 


7B. 


fc56. 


ee.? 


B<>. 


636. 


85.5 


90. 


621. 


83.5 


96. 


ft03. 


81.0 



sf.sioo: 1)1 


HEIGHT: srr 


U*4wEI&HTEl 








IH^AN 


ftltS 


«5IGM* 


V5IGM* 


cf>v»oi«vjcr 


TMrT« 


?S.790 


• 0.0 


15.560 


^ 0.0 


1.0 


0.0 


0.0 


0.0 


?5.fl6* 


t 0.005 


15.527 


I n.oo5 


0.1?* 


0.1 ?9 


O.OOT 


?ft.5 


?5.93H 


. 0.009 


15.*9ft 


> n.oo" 


n.?*3 


0.?58 


0.0?8 


?5.1 


?6.0l? 


t 0.013 


l5.*bR 


• 0.01* 


0.3511 


0.385 


0.061 


>5.5 


P6.09* 


? 0.017 


15.*** 


L 0.019 


ri,*67 


0.509 


0.106 


?5.9 


?6.?I6 


' 0.0?* 


I5.*09 


o.n?7 


0.66* 


0.737 


0.??0 


?6.5 


?6.337 - 


t 0.031 


15.390 


. 0.0 35 


O.Bftf, 


0.9*6 


0.160 


?6.7 


?6.*52 . 


! 0.037 


15.386 < 


. 0.04? 


1.0?0 


1.135 


0.5?1 


?6.6 


?6.557 < 


>■ O.O** 


15.38? i 


. 0.0*7 


l.l«5 


1.P77 


n.66? 


?5.? 


?6.670 ' 


• 0.050 


15.39* . 


. 0.051 


1.35* 


1.385 


0.785 


?3.? 


?6.77? < 


' 0.056 


15.*03 ; 


. 0.053 


1.501 


1.4*1 


0.899 


?n.9 


gb.Hbb i 


> 0.060 


I5.*37 ; 


. 0.056 


1.513 


1.507 


0.98? 


70.7 


?6.9T2 . 


► 0.065 


I5.*76 < 


. 0.0 59 


1.716 


1.566 


1.09*, 


?n.* 


?7.077 - 


. 0.068 


15.5?7 . 


. 0.06? 


1.799 


1.635 


1.17? 


19.9 


?T.15e < 


► 0.07? 


15.531 • 


0.06? 


1.863 


1.616 


I.?03 


|9.1 


?7.?50 . 


► 0.075 


15.58* ; 


. 0.063 


1.9? 3 


1.6?5 


l.?T4 


19.0 


37.337 . 


0.078 


15.670 . 


. 0.065 


1.960 


l.6*J 


I.34S 


19.3 


27.*09 . 


. O.OBO 


15.736 . 


. 0.06« 


1.986 


I.».87 


1.355 


19.0 


?7.i.86 ! 


O.DAI 


Ii.793 ; 


. 0.070 


I.99M 


I.7?5 


l.*39 


19.8 



SOU«CE ST»TI0N IS: 1 - ONT SUORUOt 
TIME MEIGHT « REM*|MlMG ME*^ CONC. 




f 
i 
& 

cr 

O 

(I 

c 



3, 
6. 

9. 
12. 

18. 
?*. 
30. 
36. 
*?. 
*B. 
5*. 
60. 
66. 
7?. 
78. 
8*. 
90. 
96. 



1*8*. 9*5 

l?|9.03? 

100 7.9?* 

837.65* 

699.171 

*«8.59? 

3*0.661 

?17.5** 

lf.6.587 

I17.**J 

8?. 13a 

56. 8*9 

39.610 

?7.558 

18.9*? 

1?.99* 

8.733 

5.9o? 

3.877 



100. 

8?. I 

57.9 

56.* 

47.1 

3?. 9 

23.9 

16.0 

11.? 

7.9 

5.S 

3.S 

1,1 

1^9 

1*3 

««9 

6*6 

0*4 

0.3 





PEOlOOl JJL 


HEIGHT 


: sfC 


t<i CONC. 


XME»N 


r<*tikH 


1.996 


?5.790 


r 0.0 


15.560 




0.0 


1.638 


?5.85* 


t 0.00* 


15.526 




o.no* 


1.355 


?5.939 


> 0.008 


|3.*91 




0.008 


1.126 


?6.009 


: 0.012 


I5.*57 




0.013 


0.9*0 


26.077 


r 0.017 


15,425 




0.019 


0.657 


26.201 


> 0.030 


15.36* 




0.033 


0.458 


26.307 


» 0.0*5 


15.301 




0.050 


0.319 


26.405 


» 0.067 


15.2*0 




0.070 


0.225 


?5.510 < 


► 0.09S 


15.190 




0.09? 


0.159 


?6.640 < 


' 0.I3I 


15.165 




0.117 


0.113 


?6.77b 


. 0.179 


J5.157 




0.1*3 


0.079 


26.877 . 


' 0.23? 


IS. 177 




0.179 


0.056 


26.977 . 


> 0.297 


15.209 




n.??? 


0.0*0 


27.n66 < 


> 0.377 


15.2*? 




0.?75 


0.028 


27.138 '< 


. 0.47T 


15.25I 




0.33* 


0.020 


27.208 < 


> 0.599 


15.306 




0.407 


0.014 


27.2?9 . 


■ 0.751 


15.359 




0.500 


0.010 


27.^53 < 


0.9 38 


15.390 




n.62l 


0.006 


?7.?55 < 


1.175 


li.3B7 




0.770 



5ULPMJ0 DIOXIDE t«EIGHTED 
<S1G«« rslOM« COV»RI»HCE 



0.0 

n.l?2 

0.239 

0.351 

0.458 

0.656 

0.84S 

1.031 

I. 2? I 

I.*?* 

1.6? I 

1.753 

1.869 

1.979 

?.0T5 

2.158 

2.219 

?.279 

2.313 



0.0 

O.I?9 

n.?58 

0.38* 

0.506 

O.T?9 

0.9?? 

1.081 

1.19? 

I.?67 

l.?99 

1.3*7 

1.395 

l.**l 

l.*53 

l.*69 

1.478 

1.507 

1.515 



0.0 

O.OOT 

0.027 

0.058 

0.O99 

0.?03 

0.3?3 

0.4*5 

0.550 

0.688 

0.8?7 

n.999 

1.1*6 

1.2 38 

1.305 

1.381 

l.*77 

1.53? 

1.533 



THET* 



0.0 
2*. 5 
?*.9 
25.1 
25.? 
75.? 

7*.* 

22.7 
?0.6 
18.8 
17.5 

18,0 
18.? 
17.6 
16.9 

t6.S 
16.7 
16.* 

16. n 



SOUHCf STATION IS! 1 - ONT SUOHU«r 
TIME I*E1GHT t I»EM«|«*ING ME*N CONC. 



SEPIOD! JJL 



MEir.HTi <;rc 
y»<E»>* 



«« 

3:4. 

6. 

9. 
I?. 
10. 
24. 
30. 
36. 
•Z. 
«8. 
5*. 
60. 
66. 
72. 
70. 
8*. 
90. 
96. 



1*2.176 
215.956 
27*.8«7 
3l9.7ft6 
SSI. 233 
308.90* 
307,37* 
376.652 
3S9.33S 
319.271 
3n8.3|9 
272.2*3 
2*2.379 
210. 7?6 
183.383 
ISS.Sll 
177.062 
lnS.*28 
79.637 



lOO.O 
1S1.9 
193.3 
22*. 9 
2*7.0 
273. S 
272.5 
26*. 9 
252.7 
238.6 
216.9 
191.5 
170.5 
1*8.2 
129.0 
109.* 
89.* 
7*. 2 
56.0 



0.191 
0.290 
0.369 
0.*30 
0.*72 
0.523 
0.521 
0.506 
0.*8S 
0.*60 
0.<>23 
0.380 
0.3*3 
0.302 
0.272 
0.237 
0.200 
0.170 
0.132 



25.790 
25.«67 
25.9*1 
26.010 
26.076 
26.192 
26.276 
26.351 
26.*73 
26.625 
26.817 
26.933 
27.017 
27.052 
27.130 
27.115 
26.9*0 
26.883 
26.8*5 



0.0 

0.008 

. 1 * 

0.019 

0.02* 

0.033 

0.0*3 

O.DS* 

0.066 

0.081 

0. 100 

0.116 

0.132 

0.152 

0.173 

0.200 

0.231 

0.263 

0.306 



15.560 
15.525 
15. ".86 
15.**D 
15.399 
15.300 
15.IB3 
15.069 
1*.996 
1*.9*B 
1*.952 
1».977 
15.025 
15.029 
15,0*0 
15.091 
15.0S8 
15.053 
I*. 998 



0.0 
n.0O9 
0.016 
n.021 

0.027 
0.036 
0.0** 
0.050 
0.055 
0.05" 
0.061 
0.069 
0.0 76 

n.om 

0.091 

n.ioi 
o.un 

0.127 
0.135 



SiiL<»n»TE ^FIGHTEO 
ISIGM* 



n.o 

0.120 

0.23* 
0.3*5 
»>.*5? 

n.65* 

0.8*7 
1.039 
1.257 
1 .49* 
1.75* 
1.915 
2.055 
2.20H 
?,3** 
2.*9o 
2.608 
2.705 
7.731 



YSJOX* 


C0V«9I»*ICE 


0.0 


0.0 


0.129 


0.007 


0.757 


0.02*. 


0.383 


0.0S6 


0.503 


0.09T 


0.711 


0.197 


0.871 


0.779 


0.975 


0.111 


1.0*7 


0.*52 


1.077 


0.616 


1.073 


0.900 


1 , 1 36 


1.125 


1.182 


1.279 


1.198 


1.357 


1.726 


l.*92 


1.257 


1.600 


1,7*1 


1.67* 


1.258 


1.683 


1.705 


1.582 



THET» 



1.0 
7*. 7 
7'5.1 
25 i1 
7S.3 
7*. 2 
71.2 
17.1 
16,0 
1S.» 
16.3 
17.1 
16.9 
15.6 
IS, 2 
I*. 5 

n.8 

13.0 
17.0 



^) 



soo«CE <;r«Tio»i isi 

TI««t «lGHt 



1 - ONT SUOHURY 



B« 


T»4. 


100.0 


3. 


t**. 


100.0 


ft. 


T**. 


100.0 


f. 


T»*i. 


100.0 


t2. 


»♦», 


100. e 


Ifl* 


T3». 


W.l 


e». 


*». 


98.1 


ao. 


TOT. 


M.« 


3»> 


669. 


89.9 


•a.. 


A«l. 


•fr.2 


«8« 


AOl. 


a«»8 


5%. 


S67. 


?».i 


6». 


S*0. 


%«« 


««» 


<;09. 


ili*» 


*Z. 


*es. 


»s-.« 


7i. 


46T. 


6Z.i 


tft* 


i>«S. 


w,§ 


to. 


(.26. 


■ff*» 


«6* 


407. 


««,? 



PEtJIOO! OCT 


HEIO-<t 


: 5rc 


uvwf ir.MTEo 






NMEAN 


Y>«E»'4 


«5ir.«» 


Y5IGH* 


COVtRUVCE 


?5.790 




0.0 


15.560 




0.0 


0.0 


0.0 


0.0 


25.B87 




0.006 


15.557 




0.00» 


0.170 


0.231 


0.003 


?5.'?87 




0.013 


15.563 




0.017 


0.34 7 


0.457 


0.007 


26.085 




0.019 


15.577 




0.025 


n.5?8 


o.*.7a 


n.009 


26.) 76 




0.026 


15.595 




0.032 


0.702 


0.874 


0.009 


?6.3ll 




0.036 


15.635 




0.045 


0.972 


1.230 


0.05* 


a6.63<i 




0.0*6 


15.654 




0.057 


1.236 


1.530 


0.091 


?6.50<) 




0.05» 


15.649 




0.066 


1.4*7 


1.753 


0.089 


26.530 




0.063 


15.620 




n.07l 


1.6*0 


1.888 


. 1 38 


^6.5'i3 




0.07a 


IS.SfcO 




0.080 


1.823 


2.016 


0.166 


26. 471 




0.080 


15.471 




0.085 


1.969 


2.094 


0.294 


26.387 




0.089 


15.430 




0.090 


2.115 


2.137 


0.612 


?6.?72 




0.095 


15.330 




0.093 


2.204 


2.168 


0.721 


26.125 




0.100 


15.177 




0.095 


2.250 


2.1*8 


0.671 


26.037 




0.104 


15.018 




0.097 


2.302 


2.1*9 


0.720 


25.927 




0.107 


14.838 




0.099 


2.308 


2.130 


0.679 


25.899 




O.IIO 


14.765 




0.099 


2.338 


2.103 


0.805 


25.790 




0.115 


14.677 




0.101 


2.365 


2.088 


0.991 


25.680 




0.118 


14.592 


£ 


0,103 


2.378 


2.069 


1 . 1 .1 



TMTT* 



ft.« 
3«S 
1.9 
1*1 
3.3 
3«* 
?.« 
2*9 
?*» 
8>.3 
T.« 
9,« 
T.» 
t.T 
t,3 
«.ft 
tOttl 
18*1 




n 

o 

D 

n 
o 
o 

C) 



SOURCE 


ST»TI0H IS! 


- ONT suoswr 


PESIOOl OCT 


HEIG'^TI %rc 


5Ut«>HUl» 


DIOXIDE WEIGHTED 




tinc 


WEIGHT 4 


9E««*lNtM6 


ME»>( CO«*C. 


IXEJN 


YMEAN 


«5tG<4» 


YSIGMA 


C0V»OUNCE 


THET« 


0* 


1895.6 75 


100.0 


2.548 


25.790 


' 0.0 


15.560 


t 0.0 


0.0 


0.0 


0,0 


8,8 


3. 


1588.457 


83.8 


2. 135 


25.889 


. 0.004 


15.558 


i. 0.006 


O.JTO 


0.229 


0.004 


ni 


6. 


1341.188 


70.7 


1.803 


25.989 


• 0.009 


15.565 


t 0.012 


0.3*7 


0.k5l 


o.oti 


5.1 


9. 


1135.314 


59.9 


1,526 


26.087 


t 0.016 


15.57* 


t 0.02O 


0.529 


0,658 


0.017 


3.5 


12. 


9f.2.373 


50.8 


1.294 


26.176 


* 0.023 


15.578 


t 0.027 


0.703 


0.849 


0.026 


3.8 


If. 


691.601 


36.5 


0.938 


26.319 


► 0.037 


15.588 


I 0.P45 


0.973 


1.193 


0.092 


5.5 


e«. 


499.890 


26.4 


0.685 


26.445 


J O.055 


15.581 


I 0.066 


1.228 


1.485 


0.128 


*.9 


10* 


356.020 


18.8 


0.504 


26.513 


► 0.075 


15.561 


L 0.091 


l.*22 


1.711 


0.082 


8*3 


ift. 


249.857 


13.2 


0.373 


26.556 . 


' 0.103 


15.508 


. 0.117 


1,630 


1.853 


0.113 


*.* 


W2* 


IT*. 913 


9.2 


0.273 


26.565 


0.135 


15.400 


. 0.149 


1.785 


1.976 


0.134 


Z.« 


«l« 


122.418 


6.5 


0.204 


26.581 < 


> 0.175 


15.3 30 


. 0.188 


1.9*0 


2.082 


0.289 


«»* 


S*» 


85.524 


4.5 


0.151 


26.SBB < 


' 0.225 


15.270 


. 0.2 34 


2.085 


2.165 


0.595 


T.8 


««. 


59.279 


3.1 


O.IIO 


26.479 ! 


- 0.278 


15.177 . 


. 0.294 


2.139 


2.260 


0.711 


• .A 


«ft» 


40.889 


1. 1 


0.080 


26.358 < 


0.338 


14.976 . 


. 0.358 


2.159 


2.288 


0.628 


r.T 


7B. 


29.040 


1.5 


0.060 


26,322 < 


' 0.408 


14.718 i 


. 0.42f. 


2.198 


2.298 


0.632 


».5 


Ti, 


20.895 


1.1 


0.045 


26.226 ; 


0.472 


14.497 . 


0.491 


2.156 


2.244 


0.569 


T.n 


8*. 


15.594 


0.8 


0.0 35 


26.239 . 


0.54 7 


14.430 ; 


] 0.568 


2.158 


2.243 


0.780 


9*S 


91>* 


11.383 


0.6 


0.027 


26.148 . 


0.650 


14.359 . 


. 0.657 


2.193 


2.217 


0,977 


tl.S 


«br.: 


8.2?l 


0.4 


0.020 


26.001 . 


0.772 


l'..2(*0 . 


0.7,2 


2.214 


2.156 


0.844 


■».« 



P 
o 

o 



SOUBCF ST*TIO^ IS I 1 - ONT SUD8JBY 

TME WEIGHT » OE«*»|MlNG ME»^ CO»(C. 



PEHIOQ! OCT 

<HE«1 



HEIGHT! SFC 



0. 
9. 

6. 

9. 
12. 
16. 
2«. 
30. 
96. 
*?. 
4B. 
5*. 
60. 
66. 
li. 
78. 
B*. 
90. 
96. 



llS.&AO 
n'>.6| J 
ISJ.Hkl 
R3.1?l 
K8.261 
|<i5.3(>3 

l«;3.*l* 

1I3.1TT 
90.371 
75.365 
62.305 
*9.755 
3S.645 
13. 7*,; 
10.0A2 
?7.*T5 
?*.6SI 
?1.938 



loo.o 

113.6 

t?9.e 

137.7 

1*2.0 

139.6 

129.5 

ll<i.7 

95.5 

76.3 

63.6 

52.6 

<i2.0 

32.6 

28.5 

25.* 

23.2 

20.8 

IS. 5 



0.159 
0.181 
0.207 
0.219 
0.226 
0.22* 
0.210 
0.192 
0.169 
0.1*1 
0.125 
0.110 
0.09? 
0.076 
0.069 
0,06* 
0.061 
0.058 
0.05* 



25.790 
25.896 
25.99* 
26.0 76 
26.158 
26.30* 
26.*33 
26.517 
?6.6** 
26.762 
26.955 
27.116 
26.9*6 
26.898 
26.829 
26.652 
26.598 
a6.**6 
26.29<, 



0.0 

0.015 

0.077 

0.0*0 

0.051 

O.075 

0.097 

0.118 

0.151 

0.183 

0.216 

0.263 

0.29? 

0.332 

0.35 7 

0.367 

0. JT5 

0.318 

0.k|7 



15.560 
15.572 
15.599 
15.5*0 
15.559 
15.502 
15.385 
15.265 
15.06* 
1*.7 76 
1*.6*2 
1 * . 559 
1*.382 
13.95* 
13.6** 
13.382 
13.376 
13.311 
13.289 



0.0 

o.n|o 

0,035 
O.O'-.B 
0.061 
0.097 

n.i 12 
n.i39 
n. icib 
0,I8A 
0,221 
0,255 
0,307 
0,3*1 
0.159 
0. 3*7 
0.376 
0,369 
0.374 



5««.PH*TF ifFtOHTCO 
ISIG"* 



n.o 

O.I6'« 
O.IS" 

n.5|7 

n,690 

n.96i 

1.20? 
1.379 
1 .605 
1.737 

2.07? 
2.06? 
2.06* 
?.07S 
2.012 
1.967 
1.926 
1.952 



VSIO"** 


COVtOUMCE 


O.n 


0.0 


0.220 


n.005 


n.*i| 


0,020 


0.616 


0.076 


0.785 


0.030 


l.M* 


0.085 


1.193 


0,101 


1.622 


O.0O1 


1.7*1 


0.B7I 


1,790 


0,109 


1.917 


0.381 


2.009 


0.8*0 


7.169 


0.905 


2.122 


l.i*" 


2.085 


I.I** 


1.901 


0.803 


1.969 


i.on 


1.^32 


1.700 


1.75* 


1 . 385 



TMET« 



0.0 

10. * 

10.1 

S.6 

1.6 

5.7 

*,1 

0.1 

l-». 

2.1 

6.2 

11.1 

12.0 

15.1 

1*.9 

11.2 

1*.7 

17.9 

70.0 




SOUWCF STtTlOW IS! 1 - ONT SUOBKOY 
TIXE KEtGHT * ffMilvlNG 



»» 


?«.6. 


100.0 


3> 


?»5. 


loo.o 


6* 


?*6. 


100.0 


9i 


?«6. 


loo.o 


\i. 


?'.6. 


100.0 


19. 


?**. 


99.? 


?<•. 


7*0. 


97.6 


30. 


?38. 


96.? 


3«.. 


735. 


9S.5 


«?. 


7ae. 


91. r 


«B. 


716. 


97.9 


5*. 


707. 


g(>.i 


60. 


701. 


91.7 


66. 


I^S. 


79.3 


rz. 


198. 


76. « 


7B. 


190. 


73.? 


e<>. 


IT*. 


70.7 


90. 


159. 


6fl.7 


96. 


161. 


65.4 



o^a 


lOUt 7901 


•^EIG^T 


1 5rc 


IJMHEI3HTED 








««F»N 


t<*fltU 


«5ir.>«» 


T<;ir,M« 


C0V»III«»*CE 


r-ifrti 


?=i.790 


♦ 0.0 


5.550 


* 


0.0 


0.0 


0.0 


n.o 


0.0 


?S,795 


. 0.017 


S.*3* 


« 


0.017 


n.369 


n.769 


-fl.017 


-17.9 


?5.B3« 


r 0.03<. 


5.385 


♦ 


0.033 


0.539 


0.516 


-0.073 


-16.6 


25.913 


• o.oso 


5.1*9 


* 


0.0*7 


0.793 


0.7 39 


-0.173 


-15.7 


?6.009 


* 0.066 


5.017 


♦ 


0.060 


1.033 


0.937 


-0.303 


-15.9 


ah.?03 


• 0.095 


I..779 


* 


0.091 


l.*R9 


1.753 


-0.556 


-l*.l 


?6.H.l 


* 0.130 


*.565 


* 


0.099 


1.95* 


1.519 


-0.655 


-10.7 


36. (.ft? 


« 0.1*3 


*.315 


* 


0.11* 


3.196 


1.763 


-0.57* 


-9.0 


?6.i.59 


♦ 0.159 


».0*7 


« 


0.1 JO 


7.** 3 


1.993 


-0.55? 


-5.3 


26. 399 


♦ 0.1 79 


3.815 


• 


0. 1*6 


?.69» 


3.709 


-0.199 


-1.6 


36.152 


. 0.193 


3.516 


* 


0.15* 


7.933 


7,755 


-0.*09 


-?.9 


75.966 


♦ 0.313 


3.*06 


« 


0.155 


3.071 


3.735 


-0.195 


-1.? 


35.785 


• 0.33* 


3.333 


♦ 


0.150 


3.330 


?.?73 


-0.1*3 


-0.7 


35.513 


♦ 0.3*9 


3.353 


* 


0.151 


3.*90 


3.3*8 


-0,695 


-1.3 


35.150 


♦ 0.359 


3.371 


* 


0.150 


3.557 


3.190 


-1.396 


-6.3 


3ft. 799 


♦ 0.357 


3.393 


» 


0.151 


3.579 


3.155 


-3.160 


-9.6 


3<..593 


• 0.371 


3.595 


• 


0.159 


3.577 


3.735 


-3.730 


-13.0 


3'.. (.89 


♦ 0.373 


3.H39 


* 


0.185 


3.550 


3.*0I 


-l.*l? 


-15,1 


3'.. 531. 


• 0.37* 


3.995 


* 


0.196 


3.*79 


3o*99 


-3.633 


-16,7 



500BCE 5T«TlOH IS: 1 - ONT 5UDBUBY 
TIME »*£IGMT % 0£M«IMINO ME»-J COMC 



q 
I 

i 



3. 

6. 

9. 
13. 

la. 

2*. 

30. 
36. 
*2. 
*8. 

5*. 
60. 
66. 
72. 
79. 
9*. 
90. 
96. 



1013. 330 

873.93* 

752.218 

650.007 

563.969 

*?3.155 

316.890 

3*0. *70 

1B1.*1' 

13*. 335 

99,*37 

73.656 

5*. 2*7 

39.291 

78.386 

70.*65 

15.199 

11.313 

9.385 



100.0 

86.1 

7*. 2 

6*.l 

55.7 

*1.8 

31.3 

23.7 

17.9 

13.3 

9,9 

7.3 

5.* 

3,9 

3,8 

2,0 

1.5 

1.1 

0.9 





PEBIOOI 7901 


HEIGMTl ^fC 


»-J COMC. 


X«E*N 


rti.t'* 


*.119 


35.790 




0.0 


15.560 


♦_ 


0.0 


3.5*9 


35.797 




0.009 


15,*37 




0.009 


3.059 


35,939 




0,019 


15,301 




0.019 


2.6*2 


35.903 




0.030 


15.191 




0.03'< 


2.293 


35,99* 




0.0*3 


15.053 




O.OS" 


1.73* 


35.1*8 




0.071 


I*. 8*6 




0.058 


1.320 


76,257 




0.103 


1».553 




O.08O 


1.010 


36.3*5 




0.139 


g*.*on 




0.106 


0.772 


36,357 




0.179 


i<i.l35 




0.139 


0.589 


35.397 




0.339 


13.903 




0.177 


0.*50 


36.153 




0.395 


13.679 




0.317 


0.355 


36.030 




0.357 


n.557 




0.3*5 


0.270 


35.893 




0.**9 


13.*53 




0.393 


0.301 


35.6*6 




0.555 


13.3*5 




0.3*0 


0.151 


35. 139 




0.473 


13.327 




0.387 


0.11* 


2*. 953 




0.790 


13.391 




0.4*3 


0.087 


?*,H05 




0.919 


13.556 




0.530 


0.057 


2*. 795 




1.053 


13.761 




n.638 


0.053 


3*.9H5 




1.195 


13.99? 


t 


0.79? 



Sl)LfH'J«» DIOXIDE wEIGMTEn 
»SIG««» rSTGM* COV»RI*NCE 



0.0 

0.367 

0.573 

n.77? 

1,019 

1.46" 

1.93* 

?.155 

7.*10 

3.653 

7.9*3 

3.065 

3.399 

3.*79 

1.58* 

3.575 

3,590 

3.537 

3.4''3 



0.0 


0.0 


0.366 


-O.OIT 


0.500 


-0.07? 


0.706 


-0.163 


0.«9| 


-0.37K 


1.197 


-0.509 


1.43* 


-0.507 


1.6*6 


-0.698 


I.056 


-0.691 


3.1*9 


-0.519 


3.113 


-0.*76 


3.103 


-0.369 


3.158 


-0.?56 


3.131 


-0.B51 


3.06* 


-1.43* 


3.00* 


-1.927 


7.035 


-7.331 


3.1*5 


-3.957 


3.765 


-3.37* 



TMET» 



0.0 

-13.? 

-1*.9 

-15.? 

-15.0 

-13,3 

-10.? 

-9.* 

-6,9 

-*.? 

-3.* 

-1.6 

-1.3 

-4.0 

-6.3 

-9,6 

-10.3 

-13.3 

-15.5 






OUSCE 


<;r«TlON isi 


- om suoHUHT 


period: 7'iini 


HElij'+r: "sfc 


SIILPH4TE 


liElr.HTEn 






T11E 


WEIGHT « 


PEMIININO 


ME»N CO»JC. 


«>«E»M 


T>«E»N 


«SI&'«« 


T5I&«<» 


CO**RItNCE 


THET* 


0. 


4T.0I0 


100. 


0.191 


25.790 < 


. 0.0 


5.560 


I 0.0 


n.O 


o.n 


0.0 


0.0 


3. 


SS.26e 


123.9 


0.237 


25.802 - 


> 0.0 34 


5.439 


. a.P3? 


n.26? 


0.?44 


-0.016 


-13.2 


6. 


K.8.205 


145.1 


0.277 


25.B26 « 


0.052 


5.327 . 


. 0.053 


0.514 


0.438 


-0.051 


-11.0 


9. 


T5.SB0 


160.8 


0.JO7 


25.fl^2 < 


• 0.087 


5.217 


. 0.059 


0.756 


0.600 


-n.llR 


-11.6 


12. 


HO. 598 


171.5 


0.32B 


25.926 . 


0.111 


5.105 


^ 0.094 


0.996 


0.758 


-0.199 


-11.4 


18. 


(16.916 


184.9 


O.JSft 


26.083 • 


' 0.157 


4.898 


. 0.110 


1.461 


1.025 


-0.192 


-10.4 


g". 


«9.353 


190.1 


0.372 


25.191 . 


. 0.195 


4.705 


L 0.130 


1.845 


1.230 


-0.436 


-7.1 


3D. 


n9.6ea 


190. B 


0.377 


26.296 < 


• 0.230 


4.43S 


. n.156 


2.182 


1.474 


-n.519 


-6.2 


36. 


«6.6?3 


184.3 


0.369 


26.315 < 


0.263 


4.122 


L n.ie? 


2.44ft 


1.696 


-0.645 


-6.? 


<•?. 


no.rfci 


171. « 


0.35* 


fb.eiili « 


0.299 


3. 87? 


I n.?09 


?.690 


1.874 


-0.66* 


'S,? 


<>8. 


75.788 


161.2 


0.351 


26.177 4 


0.332 


3.669 


L 0.221 


2.890 


1.925 


-0.671 


-4.6 


5*. 


69.5«2 


147.9 


0.336 


26.1 17 . 


0.374 


3.556 


n.?28 


3.119 


1.905 


-n.584 


-1.4 


60- 


60.022 


127.7 


0.299 


26.045 . 


0.438 


3.404 


. n,251 


3.391 


1.941 


-0.852 


-4.? 


66. 


S3. 550 


113.9 


0.275 


25.950 • 


0.490 


3.269 


. n.?5S 


3.587 


1.853 


-1.506 


-6.7 


TJ. 


*7.779 


101.6 


0.254 


25.873 . 


0.540 


3.280 < 


0.258 


3.731 


1.785 


-1.560 


-6.4 


78. 


*1.I69 


87.6 


0.229 


25.544 . 


0.578 


3.379 


n.253 


3.707 


1.625 


-1.562 


-6.9 


8<>. 


37.196 


79.1 


0.214 


25.583 t 


0.59 7 


3.547 


> 0.236 


3.643 


1.560 


-1.771 


-7.6 


90. 


33.776 


71.8 


0.200 


25.B72 . 


0.508 


3.726 


. 0.27B 


3.535 


1.616 


-2.089 


-9.5 


96. 


31.258 


66.5 


0.194 


26.097 . 


0.612 


4.040 . 


I 0.31? 


3.424 


1.745 


-2.396 


-11.5 







V 



SOU»CF <;T»TlOn ISi I - ONI >;ui)HURr 

TIMt KEIOHT % OEMAl'^ING 



»• 


?*8. 


loo.o 


3. 


7<>8. 


100.0 


6. 


?«8. 


loo.o 


9. 


7«8. 


loo.o 


12. 


?*8. 


100.0 


18. 


?*8. 


100. D 


?*. 


7«8. 


loo.o 


30. 


?*8. 


100.0 


3*. 


?*8, 


loo.o 


*?. 


?*8. 


100.0 


08. 


?«8. 


loo.o 


Sft. 


7«8. 


100.0 


60. 


?4e. 


loo.o 


66. 


?*8. 


loo.o 


72. 


?*l. 


97.? 


78. 


?38. 


96.0 


8<i. 


?3r. 


95.6 


<>o. 


?36. 


95.? 


96. 


?33. 


9<..0 



at" 


[001 THOI 


HEir,H' 


: SFC 


IJMMEI-.HTEn 








X" 


E»N 


rMt*V 


»S15M» 


Y5IQM* 


OVMUSCE 


THEt» 


?5.790 


t 0.0 


15.560 < 


0.0 


0.0 


o.n 


0.0 


0.0 


?5.«71 


; 0.007 


1S.5SH < 


n.no« 


o.nc 


'».)?7 


n.oo* 


18.6 


?S.951 


>■ 0.010 


15.559 


n.oi6 


0.?19 


1.75 3 


n.017 


19.7 


?6.0?9 


' 0.0?1 


15.56? 


o.o?<. 


«).3?5 


n.37fl 


0.01" 


19.8 


?6.10? 


• 0.0?7 


15.567 . 


. n.0 3? 


n.<i?«. 


0.<.97 


0.06" 


»0.5 


?6.?37 


♦ 0.039 


15.5H5 


t n.o*6 


n.611 


0.7|7 


O.l*" 


?!.<. 


?6.366 


• D.0<.9 


15.610 


. o.oso 


n,7T9 


0.913 


0.751 


?7.» 


?6.495 


» O.059 


15. 60'. 


> 0.069 


0.93? 


l.n«2 


0.361. 


77.8 


?6.6?9 


•■ 0.06*) 


15.685 


. 0.078 


1.070 


l.??3 


0.'.77 


77.6 


?6.77I 


♦ 0.076 


15.7 35 


. 0.085 


1.19? 


1.318 


0.579 


??.? 


?6.9?? 


; 0.08? 


15.796 


0.091 


I.?9P 


l.^kO 


n.67| 


71.7 


?7.0 76 


• 0.088 


15.868 


. 0.fl9B 


1.J9P 


1.S39 


0.751 


71.3 


?7.??6 


• 0.09 J 


15.950 


. O.lOfc 


l.'.66 


1.635 


n.sn 


70.7 


?7.369 


♦ 0.097 


16. Ok? 


. 1.110 


I.5?9 


1.7?5 


0.8*5 


19.9 


?7.5?<. 


• 0.103 


16.0?5 


. 0.109 


1.60? 


1.686 


0.97? 


19." 


?7.683 


• 0.107 


16.081 


. 0.110 


1.65) 


1.693 


0.96« 


19.6 


?7.H36 


• 0.109 


16.1 he 


. n.m 


1.680 


1.7 3* 


0.957 


18.7 


?7.>)79 


1 0.110 


16.31* 


- 0.II5 


1.69% 


1.77? 


0.898 


|7.3 


?8.15| 


» 0.10^ 


I6.<>?e 


> n.ii7 


1.651 


l.T8» 


0.973 


19.7 



i 

8 
a 



OUBCE 


STATION ISl 


1 - ONT SUDBJPr 


OE9100I 7907 


MEIG-*?! 5FC 


SOLPH'JH 


OIO«IDE KEIGHTEO 




THE 


WEIGHT 


% REMAINING 


mEA^ CONC. 


XMEAH 


YMEAN 


ISJOIA 


rstOMA 


covAniAtcE 


Ti^ETA 


0> 


*9*.98? 


100.0 


1.996 


?5.790 


> 0.0 


15.560 


I 0.0 


n.o 


0.0 


0.0 


0.0 


3. 


*n6.6?8 


8?. 2 


1.6*0 


75.870 


► 0.005 


15.558 


• 0.006 


0.110 


0.177 


0.00» 


19.4 


6. 


315. 7?7 


67.8 


1.J5* 


?5.9*5 


► 0.01? 


15.555 


. 0.01* 


n.?i8 


0.75* 


n.OI6 


18.4 


9. 


?77.*37 


56.0 


1.119 


?6.0i6 


; 0.019 


15.550 


L 0.0?3 


0.373 


0.379 


0.035 


18.* 


13. 


J79.91* 


1.6. * 


0.9?7 


76.08? 


. 0.078 


15.544 


L 0.033 


n.fc?* 


0.500 


0.060 


18.6 


18. 


|58.<.7I 


3?.0 


0.639 


?5.|98 


; 0.0*9 


15.578 


'_ 0.057 


0.51? 


0.718 


0.177 


18.8 


?*. 


l(ie.9?6 


??.0 


0.*39 


?5.?91 


> 0.076 


15.505 


t 0.087 


0.790 


0.904 


0.701 


17,8 


30. 


75.09* 


15,? 


0.303 


?6.3 77 


• o.no 


15.481 


t 0.1?? 


0.956 


1.056 


0.267 


16.3 


36. 


51.967 


10.5 


0.?I0 


?6.*77 . 


' 0.153 


15.*70 < 


l 0.15? 


1.105 


1.171 


0.338 


15.5 


*?. 


35.965 


7. J 


0.1*5 


76.515 


" 0.?06 


15.473 


► 0.70 7 


1.7 36 


1.?** 


n.*i4 


15.2 


*e. 


7*.61!» 


5.0 


0.099 


76.79? 


> 0.771 


I5.*e6 


•_ n.?59 


1.3*3 


1.787 


0.50? 


15.5 


5*. 


17.019 


3.* 


0.059 


?6.99? 


. 0.3*4 


15.507 


. o.3?n 


1.*I7 


1.319 


0.58* 


16.2 


50. 


1 1 . 955 


i.o 


0.0*8 


?7.|79 


► 0."*?* 


15.551 


t 0.»95 


l.*66 


1.166 


0.630 


16.3 


66. 


e.3?5 


».» 


0.03* 


?7.335 < 


0.;?? 


15.60* . 


n.*96 


1.506 


1.437 


0.6*8 


15.9 


7?. 


5.569 


1.1 


0.0?* 


?7,*96 . 


0.65? 


15.6?? ; 


. 0.618 


1.55? 


1.470 


0.68<. 


15." 


78. 


3.857 


a.e 


0.016 


27.568 « 


0.813 


15.5h8 . 


0.78? 


1.597 


1.5 36 


0.719 


15.7 


8*. 


?.588 


o.s 


0.011 


?7.857 


1.009 


15.770 ; 


0,99* 


1.6?4 


1.599 


n.79i 


16.7 


90. 


\.'?b 


0.3 


0.007 


?8.oja . 


■ l.?3? 


15. 7M? . 


I.?5» 


1.619 


1.65? 


0.865 


|8.3 


96. 


l.l?7 


0.? 


0.005 


?8.|7* 


1.501 


15.856 ; 


1.597 


1.591 


1.695 


0.926 


70,1 



SOUBCF <;T»TIOn IS; 1 - ONT SUOHUfl* 

TIME HEIGHT « OEM»Ii<ING ME»X C0»4C. 



OEPIQUS TRO; 
XHE»i« 






3« 
«» 
■*» 

so* 

«§• 
^» 

m* 

m* 

'ft* 



(>T.3q2 
73.576 
01.977 
105. «1? 
11*. i! 3? 
1?*.02J 
119. SS7 
11*. SO"* 
lnl.97* 

7?.8i.e 

^6.3<i2 

s'>,07r 

*7.<>65 

30.962 
?2.*6e 
15.715 
|0,1B3 



lOO.O 

153.1 

19I..1 

2?3.3 

2<>I.O 

2h|.7 

252.3 

2*1.6 

215.2 

1"»0.2 

153.7 

UD.O 

12*. 7 

101.2 

Bl.R 

65.2 

<>7.* 

33.2 

21.5 



0.191 
0.293 
0.371 

o.*2r 

0.*61 
0.500 

0.*62 
0.(>1 1 
0. J63 
0.29<. 
0.268 
0.23H 
0.193 
0.161 
0.130 
0.095 
0.06' 
0.0** 



25.790 
25.R7I 
25.«*5 
26.01* 
26.079 
26. ISO 
26.222 
26.?3» 
26.? 79 
36.'. 19 
26.668 
?6.96'. 
27.169 
27.3*r 
27.623 
27.812 
28.092 
28.350 
28.532 



0.0 

0.013 

0.02 3 

0.031 

0.0*0 

0.056 

0.0711 

0.090 

O.IOB 

0.1P7 

0.15S 

0.17* 

0.187 

0.20T 

0.225 

0.262 

0.305 

0.J59 

0.465 



15.560 
15.558 
15.55* 
15.5*2 
15.519 
15.ft70 
15.396 
15.339 
15.273 
15.3'.'. 
15.251 
15.368 
15.475 
15.529 
15.557 
15.636 
15.502 
15.517 
15.631 



0.0 

0.015 

O.O?*. 

n.037 

O.OfcT 

0.065 

0.081 

o.m 

0.131 

n.i?9 

0.136 

n.it.5 

0.165 
0.190 
0.2J5 
0.383 
0.377 
0.*3I 



SciLOM^TE JEir.MTED 



0.0 

n.lio 

n.218 

n.321 

0.*?3 

0.63? 

0.806 

n.959 

1.091. 

I. 311) 

1.32* 

l.*20 

l.*38 

l.*36 

I.*0? 

l.*58 

l.**5 

l.*2? 

l.*81 



YS1CH» 


COV»«»I»HCE 


0,0 


0.0 


0.127 


O.OD* 


0.353 


O.OIS 


0.176 


0.03? 


0.I.9R 


0.058 


1.719 


0.)?'' 


0.901 


0.156 


1.0*0 


0.126 


1.137 


0,15? 


1.1*7 


0.356 


1.105 


0,525 


1.108 


0.595 


1.11* 


0.*80 


1.1** 


0.383 


1.181 


0.*93 


1.308 


0.5*9 


1.339 


1.091. 


l.*96 


1.1*1 


1.176 


o.se* 



TWFT« 



t«*9 

9.t 
Ift.T 
16.4 
13.1 
10.S 

t*.l 

?T.r 
?»** 







SOOi^CF STATION IS: ?0 - ILL 9»L0i»I'< 



0. 


2BTa. 


100.0 


3. 


2bT2. 


100.0 


6. 


?R7?. 


lon.o 


9. 


a«TE. 


100.0 


1?. 


?«T2. 


100. 


j». 


?a7g. 


100.0 


?*. 


?l^b?. 


9<i.T 


30. 


20*2. 


99.0 


36. 


2«?*. 


98.3 


*2. 


?792. 


9T.2 


«iS. 


^T<.^. 


95. S 


5*. 


?ft9?. 


93.7 


60. 


S'.SI. 


92.3 


66. 


?S88. 


90.1 


12. 


2S22. 


87.9 


78. 


2u37. 


811.9 


8«>. 


23<>l. 


SI.S 


90. 


??T?. 


79.1 


96. 


?199. 


76. Z 



PEBIOOl 197H 


HEl&Hl: srr 


UNHEf.MTEi) 








X«<t«N 


YHEAN 


««;IGM« 


YSIGM* 


cavM>itsct 


ThET* 


?3.»I0 


> 0.0 


13.110 


> 0.0 


0.0 


0.0 


n.o 


0.0 


?3.'.?<. 


•■ 0.00<> 


13.116 < 


• 0.005 


n.206 


0.?*? 


0.003 


*.3 


?1.<.3fl 


> 0.009 


11.129 


. 0.009 


O.*06 


0.*6H 


o.on 


*.6 


23.<.»9 


' 0.011 


13.l<>9 


0.01 1 


0.595 


0.673 


0.0 35 


5.6 


?3.'.56 


► 0.015 


13.175 


t).01». 


0.778 


0.865 


0.070 


6.S 


?3.'.66 


> 0.021 


13.2*3 < 


. 0,022 


1.119 


1.201 


n.i«3 


8.3 


?3.<.78 


• 0.027 


13.350 < 


. 0.027 


l.*29 


I.*65 


"1.317 


8.8 


23.501 


. 0.032 


I3.*9H < 


. O.OJI 


1.70 7 


1.661 


a.*l3 


8.1 


23.52R < 


► 0.037 


13.658 


0.0 J5 


1.9*5 


1.835 


0.*55 


6.9 


23.S75 - 


I- 0.0<>1 


13.827 < 


. 0.037 


2.1*5 


1.967 


0.*39 


5.5 


23.665 


► O.Ofci. 


13.968 


0.039 


2.297 


2.062 


0.»36 


*.7 


23.r31 


> 0.0*7 


1*,096 < 


. 0.0*0 


2.4*2 


2.099 


0.373 


3.6 


23.181 


►■ 0.050 


1*.227 < 


0.0*2 


2.579 


2.1*6 


0,202 


1.7 


23. "60 


. o.os? 


1*.333 


. 0.0*1 


2.666 


2.197 


-0.063 


-O.S 


23.917 


' 0.05<> 


1*.*** 


> 0.0*5 


2.729 


2.2*3 


-0.*06 


-3.1 


23.995 


> 0.056 


1*.55J ' 


0.0*6 


2.7*7 


2.265 


-".SSI 


-*.0 


2<'.101 


. 0.057 


1*.603 


0.0*7 


2.7*1 


2.258 


-0.585 


-».5 


2fc.l99 


i 0.058 


|*.676 


. 0.0*8 


2.75* 


2.291 


-0.701 


-5.3 


2<'.285 


> 0.059 


l*.725 


0.0*9 


2.7*0 


2.297 


-0.7*9 


-5.T 



S0U9Cf S'*TI0N IS1 20 - III e»LD«is 
TIME WEIGHT % PeHtlNlNO H£«N CONC 






o 
o 

o 
o 

Q 



0. 

3. 

6. 

9. 
12. 
18. 
2*. 
30. 
36. 
*2. 
*8. 
5*. 
60. 
66. 
72. 
78. 
8*. 
90. 
96. 



7978, 
6869, 
5915, 
51*1, 
**"S3, 
3378, 
2516, 
19n2, 
1*11, 

1070. 

797, 

595. 

*<i7. 

31*. 

2*8, 

182. 

Ml, 
99, 
7*, 



269 

870 
*67 
6(S* 
0r,0 
062 
676 
0*8 
813 
*b2 
9iS9 
65* 
195 
551 
813 
02* 
-I' 

ruo 

393 



100.0 

86.1 

7*.* 

6*.k 

55.9 

*2.3 

31. B 

23.8 

17,9 

13.* 

10. 

7.5 

5.6 

4.2 

3.1 

2.3 

1.7 

1.3 

0.9 





("EOIOOI 1978 


height: s^r 


«N CONC. 


XME*N 


Y>*E»N 


2.778 


23.410 < 


' 0.0 


3.110 < 


L 0.0 


2.392 


21.*25 


• 0.003 


3.093 


L 0.003 


2.067 


23.**0 < 


. 0.005 


3.082 


L 0.006 


1.790 


23.452 - 


> 0.008 


3.0 76 


I O.OIO 


1.55* 


23.*60 


• 0.012 


3.076 


t 0.011 


1.176 


23.46* 


> 0.020 


13,092 


L 0.021 


0.886 


23.*72 


► 0.029 


3,157 


> 0.030 


0.669 


23.*91 


• 0.0*0 


3,282 


L 0.0 39 


0.507 


23.508 


• 0.052 


l.*26 


> 0.049 


0.393 


23.539 


> 0.067 


3,582 


!! 0.061 


0.291 


23.606 - 


> 0.083 


3.703 


. O.on 


0.221 


23.656 < 


■ 0.102 


3.811 < 


0.085 


0.169 


23.685 < 


•. 0,123 


3.946 - 


. 0.099 


0.129 


23.753 . 


> 0.1*8 I 


4.056 < 


. 0.1 17 


0.099 


23.801 < 


> 0.175 


4.155 . 


, 0.140 


0.0 75 


23.89? < 


1- 0.206 


4.230 < 


. 0.164 


0.05? 


24.003 < 


' 0.238 


4.250 < 


. 0.191 


0.0** 


24.106 < 


0.275 


4. JO* < 


. 0.22». 


D.0 3* 


?<..?2'' i 


- 0.314 1 


4.110 < 


0.262 



SULOMUB OIO«IOE «#E!GMTEO 
KSIGMII V5t&M8 COVARUNCE 



0.0 

0.208 

0.409 

0,600 

n.7«8 

1,139 

l,*59 

1,7*1 

1.984 

2.184 

2.3*2 

2,492 

2.610 

2,71* 

2.783 

2.774 

?.753 

2.7*6 

2.706 



0.« 

0.252 

0.487 

0.70Q 

0,«9«i 

I. '31 

1.608 

1.68* 
I, "69 
2,002 
2,070 
2.063 
2.09 7 
2.1*9 
2.201 
2.212 
2.202 
2.255 
2.261 



0.0 
0.00* 
0.017 
0.0*5 
0.09) 
0.236 
0.398 
0.*9| 
0.472 
0.35* 
0.2*4 
0.1*5 
-O.077 
-0.411 

-0.859 
-0.99* 
-1.120 
-1.2 7? 
-1.16* 



T>^ET« 



n.o 

*.7 

5.7 

7.1 

8.4 

10.3 

10.6 

9.2 

6.8 

».2 

2,5 

1.3 

-0.6 

-3.2 

-6.3 

-7.* 

-8.* 

-9,6 

-9.0 



t) 



t'U) 



SOltSCF STiTlOft ISf ?0 - ILL H«LOi*l»^ 
TIME WEIGHT % OEtilMINO MEas CnnC. 



PERIOD: 1979 

««IE*N 



MEI&HT! ^fr 
yME«N 



0. 

3. 

6. 

9. 
12. 
IB. 
?*. 
30. 
36. 

«a. 

5*. 
60. 
66. 
73. 
78. 
B*. 
90. 
96. 



TT0.6<.9 

9?9.650 

ion*. 811 

l?nq.in3 

13«?.0O3 

1363.113 

1313. •?<• 

12';2.?63 

111.7. 1S<| 

10*9.253 

918. 1T5 

815. *72 

736.331 

6<.2.7T6 

5S5.215 

*«1.05<i 

^76.630 

375.107 



lon.n 

l<i5.<> 

ids.o 

215.9 
238.9 
267.1 
271.3 
265.* 
2*9.2 
228. 3 
208.8 
186.7 
166.3 
1*6.6 
127.9 
110.5 
95.7 
85.0 
7*. 7 



0.175 
0.25* 
0.32* 
0.378 
0.*18 
0.*67 
0.*76 
0.*69 
0.**3 
0.*11 
0.363 
0.3*9 
0.315 
0.285 
0.255 
0.228 
0.205 
0.198 

o.in 



?3.*10 
23.*36 
23. *6* 
23.*90 
23.517 
23.56* 
?3.63* 
23.713 
23.757 
23.80* 
23.883 
23.951 
?*.013 
2*. 131 
2*. 2*6 

?*.ieo 

2*. 510 
2*. 6*5 
2*. 777 



0.0 

0.007 

0.012 

0.017 

0.021 

0.028 

0.016 

0.0*3 

0.051 

0.059 

0.067 

0.075 

0.063 

0.092 

0.101 

0.108 

0.116 

0.123 

0.128 



I 1.110 
13.117 
13.132 
13.1*8 
13.162 
13.200 
13.257 
13.353 
I3.*36 
13.532 
11.6*1 
13.711 
13.800 
13.650 
13.690 
13.90? 
13.693 
13.901 
I 1.900 



0.0 
0.008 
0.01* 
0.019 
0.021 
0.029 
0.036 
0.0*1 
0.0*7 
O.Oi? 
0.057 
0.062 
0.066 
0.071 
0.077 
0.081 
0.090 
0.097 
• 0.105 



Si«.OH»TF WFIGHTEO 
»5KtM« 



n.o 

0.191 
0.373 
0.5*8 
0.719 
1.0*0 
1.323 
1.583 
1.806 
1.998 
2.169 
2.307 
2.396 
2,*9* 
2.551 
7.5** 
2.5*6 
7.531 
2.*76 



I^ICfA 


rov*Ri»itce 


0.0 


0.0 


n.778 


0.006 


n.*32 


0.025 


0.619 


0.059 


0.789 


O.ll? 


1.080 


0.250 


1.31* 


n.*28 


I.SOI 


0.515 


1.650 


0.501 


1.T69 


0.399 


1.860 


0.296 


1.892 


0.293 


1.911 


0.08? 


1.93T 


-0.190 


1.963 


-0.*9* 


1.96* 


-0.531 


1.97* 


-O.60O 


2.009 


-0.5*6 


2.03* 


-0.**1 



t-ie:t» 



0.0 

9.0 

10.3 

11.7 

12.2 

13.5 

13.7 

11.6 

8.7 

5.7 

3.6 

3.1 

O.B 

-I.T 

-6.3 

-*.T 

-S.3 

-5.8 

-*.l 



SOUBCF qI*rlo>i Is: 
TME HEIGHT 



^0 - III HiLOKiN 

« REXirNIHf, 



»• 


mat,. 


100.0 


3. 


?<^H6. 


100. 


h. 


2n86. 


100.0 


«. 


ZHflb. 


100.0 


IZ* 


?R8i>. 


99.9 


18. 


?W|8. 


»7.6 


2». 


?7&r. 


«S.9 


M. 


^Tl^. 


«*»« 


3«< 


2ftS9. 


«t*| 


*2. 


eMo. 


-«tt#* 


4ft. 


p^*/. 


W»3 


5«. 


?466. 


»!«« 


«0- 


PTJ"). 


»1.J 


^. 


?TW. 


«»»» 


fZ. 


??<»3. 


l«*is 


rgk 


?l«7. 


f%«* 


««. 


2n5<.. 


Tl»2 


w». 


1<»T0. 


68.3 


«*. 


IH'J. 


65*« 



BEPIOO: 19 7'y 


height: %FC 


UNiitrOMTEO 






»«E»M 


Y-^EtN 


»S1(}'»* 


YS!GH« 


COVkRIlVCE 


?3.*I0 


> 0.0 


n.iio 


I 0.0 


0.0 


0.0 


0.0 


PJ.*"*. 


>' 0.00<> 


1 J.l?7 


I 0.005 


0.198 


0.?53 


0.012 


?3.<.7<. 


• 0.007 


13.15S 


> 0.009 


n.*oi 


0.*95 


0.0*6 


?3.'.91 


t 0.011 


n.20! 


L 0.013 


0.60* 


O.T20 


0.100 


?3.507 


? 0.015 


13.25* 


• o.niT 


0.T98 


'>.<»21 


o.iro 


?3.S56 


> 0.0?l 


)3.<><i0 


L 0.02? 


1.11? 


1.172 


0.2T6 


?3.f>0<» 


> 0.0?6 


13.630 


►. 0.026 


1.386 


1.368 


0.*** 


?3.673 


t 0.031 


13. BI* 


> 0.0 30 


1.61* 


1.5*0 


0.552 


?3.r33 


t 0.035 


13.985 


L 0.033 


1.T99 


1.678 


0.ST7 


?3.T9« 


0.031 


ltt.141 


. 0.031 


1 .9A9 


1.T99 


n.SQii 


P3.H31 


> 0.0<.| 


lk.299 


. 0.037 


2.0*6 


I.B91 


0.365 


?3.B60 


0.0<.3 


1(>.<>60 < 


. 0.039 


?.117 


1.955 


0.235 


?3.8?5 


. 0.0*5 


I*. 597 . 


. 0.0<.J 


2.190 


2.020 


0.091 


?3.901 


>- 0.0<>& 


l*.*?". < 


0.0<.3 


2.227 


2.059 


0.179 


23.911 . 


O.Oi.8 


14.775 ' 


0.0*5 


2.263 


2.125 


0.13R 


P3.9C.6 


. 0.050 


K..85P < 


0.0*6 


?.296 


2.1S1 


0.0fl5 


?3.965 


. OS I 


1<..93| 


. 0.0<.« 


2.3?9 


2.177 


0,165 


?3.97f, 


«• 0.053 


I*. 99? 


0.0*9 


?.366 


2.183 


0.167 


23.967 i 


0.055 


15.00? : 


1 0.050 


?.*01 


2.192 


0.118 



TMfT* 



ir.e 

15,» 
1 5.5 
1^.0 
19.6 

ii.g 

19.0 
IQ.! 
T,S 
%»9 
3.0 
t*« 
?.l 
I.S 
1.9 
I.T 
1.7 
i.2 





SOURCE 


5T«T10N 1st 


20 


- Ill b*ld«in 


PEHIOOt 1979 


HElGHTi src 


SULPHUR 


OIOKIOE KEIGHTEO 






TME 


HEIGHT 


« 


OEMAINING 


ME«M CONC. 


)<««e«N 


YME«N 


><:|G>4* 


TsrOH* 


Cf)V«ni»NCE 


TMET» 




q. 


80?*. 939 




loo.o 


2.781 


23.*I0 . 


0.0 


13.110 


t 0.0 


0.0 


0.0 


0.0 


•.f 




3. 


69<,8.9flS 




86.6 


2.*08 


23.**3 


. 0.002 


13.110 


► 0.003 


0.208 


0.269 


0.01* 


Ittl 




6« 


6015. 6R6 




75.2 


2.091 


23.*69 < 


. 0.005 


13.12* 


L 0.007 


0.*21 


0.520 


0.055 


ini 




9. 


52*9. S?2 




65.* 


1.819 


23.*97 < 


' 0.009 


13.15! 


t 0.010 


0.635 


0.751 


0.121 


i«.r 




If?. 


*563.7?0 




56.9 


1.582 


23.50? < 


. 0.012 


13.189 


> 0.01* 


fi.836 


0.956 


0.207 


t«.s 




It. 


33'^7.B63 




*l.B 


1.192 


23.561 < 


. 0.020 


13.368 


L O.021 


I.I** 


1.195 


0.332 


l*.2 




■24. 


2*92.989 




31.1 


0.901 


23.613 I 


• 0.028 


13.558 


L 0.028 


l.*Ok 


1 . 38* 


0.50* 


i*.* 




90. 


1860.5*3 




23.2 


0.686 


23.668 i 


0.038 


13.730 


0.036 


1.618 


1.570 


0.609 


11*1 




3«- 


13«(a.0|3 




17.3 


0.522 


23.70* • 


• 0.0*8 


13.891 


O.O** 


1.795 


J. 722 


0.62* 


ii«« 


r^ 


*e. 


1019. 2S1 




13.0 


0.398 


23.726 • 


0.06D 


1*.0*3 


. 0.058 


1.9*0 


1.859 


0.5*1 


•.t 


m 


m. 


770.679 




9.6 


0.303 


23.715 . 


0.073 


1*.203 


. 0.071 


2.026 


1.959 


0.3** 


4.« 


an.. 


56ft. *?' 




7.1 


0.231 


?3.698 . 


0.088 


1*.367 . 


0.085 


2.10* 


2.03* 


0.20* 


*.♦ 


-_. 


«0. 


*?l.b3* 




5.3 


0.176 


?3.675 . 


0.106 


I*.512 - 


0.103 


?.|76 


?.120 


0.016 


0.? 




6«* 


308.300 




3.8 


0.133 


23.682 . 


' 0.125 


l*.60fi 


0.12* 


2.200 


2.169 


0.171 


*.o 


o 


??. 


2?7.7oI 




2.8 


0.102 


23.6*7 4 


0.1*6 


1*.673 


0.1*9 


2.206 


2.?53 


0.127 


1.9 


o 

3 


Ti. 


166. *19 




2.1 


0.078 


23.607 . 


0.169 


1*.721 


0.178 


2.185 


2.301 


0.019 


0.? 


o 


W. 


121. *?9 




1.5 


0.059 


23.586 . 


■ 0.201 


l*.7e6 


. 0.211 


?.211 


2.325 


0.086 


1.0 


o 


-»8. 


<»e.«32 




1.1 


O.Oib 


23. "^61 . 


■ 0.?37 


1*.R06 ■ 


n.?46 


P.237 


2.315 


-O.0O9 


-0.1 


S 

( 1 

o 

O 


c 


9*. 


65.6?* 




0.8 


O.OJS 


2J.S50 . 


0.279 


l*.BOI • 


0.287 


?.263 


2.321 


-0.197 


-2.2 



SOOBCE -iTlTlOM ISt ?0 - ILL H»LOwI<4 PE9I00I 1979 HEIGHft src 

TIHt KEIGMI « PEM«|NING Ht»V CO^IC. XXE*"* THEiM 



SULPHITE KEI5HTE0 



COVtRItNCe 



T«*fT» 



30:* 

ftO. 

tit 
n» 



514,876 

7<;4.ro3 

l|oT.^79 

I31B.17I 
13S9.1?4 

I3II.1S1 
1??0.999 
1132.0St> 

890. dZl 
7fl0.b9* 
6^8.696 
Sn0.63l 
*9l.*B7 

*n9,aio 
343. oir 

Z9?,l)67 



1 00.0 
1*7.5 
188.1 
ai9.3 
241.8 
26S.0 

259.7 

2<t|.8 

22*. 2 

199.2 

176.4 

154.6 

132.4 

IIS.O 

97.3 

81. 1 

67.9 

57.8 



0.175 


21.410 


0.258 


23.449 


0.329 


23.483 


0.3R4 


23.511 


0.423 


23.540 


0.475 


23.607 


0.491 


23.669 


0.493 


23.711 


0.»S9 


23.748 


0.434 


23.780 


0.395 


23.762 


0.361 


23.760 


0.325 


23.754 


0.299 


23.770 


0.259 


23.794 


0.229 


23.775 


0.199 


23,759 


0.174 


23.770 


0.154 


23.762 



0.0 

o.oor 

0.012 
0.017 
0.021 
0.029 
0.034 
0.040 
0.046 
0.052 
0.057 
0.063 
0.069 
0.076 
0.083 
0.099 
0.098 
O.IIQ 
0.122 



13.110 
13.129 
13.166 
13.211 
13.253 



13 
13 
13 
13 



.370 

,499 
.602 
674 



13.733 
13.845 
13.959 
14.063 
14.14? 
14.147 
14.181 
14.293 
14.333 
14.302 



0.0 

0.009 

O.OIS 

0.020 

0.024 

0.029 

0.035 

0.040 

0.045 

0.051 

0.057 

0.062 

0.070 

0.079 

0.096 

0.097 

0.110 

0.120 

0.1 30 



0.0 

0.I9I 

0.390 

q.567 

0.747 

1.019 

1.240 

1.432 

1.599 

1.748 

I. 911 

1.967 

1.922 

1.973 

1.991 

1.957 

1.975 

2.039 

?.0T9 



o.n 

0.236 
0.451 
0.652 
0.933 
1.0 78 
1.273 
1.449 
I. '571 
1.71)7 
1.799 
1.954 
1.954 
7.021 
2.073 
2. 145 
2.230 
2.21T 
2.730 



0.0 

O.OII 

0.041 

0.091 

0.156 

0.2T1 

0.40) 

0.*9t 

q.SOS 

0.497 

0.379 

0.301 

0.162 

0.240 

0.226 

-0.059 

>0.143 

-0.266 

-0.399 



0.0 

16.9 

15.9 

15.7 

15.6 

|4.6 

14.7 

13.5 

11.2 

9,0 

6.6 

4.9 

2.5 

3.5 

3.3 

-0.9 

-2.1 

-1.7 

-5.1 




SOUWCF "JTSTIO*! ISt BP - ILL BALDwrN 

TIME uttOHl « OE>**1^J|NO 



0* 


2«9*. 


100.0 


3. 


2B94. 


100.0 


«• 


?«<>*. 


loo.o 


9. 


?»?<.. 


loo.o 


12. 


gB"".. 


loo.o 


I«. 


PBTft, 


<i<).i 


2*. 


^«^7. 


<*1.1 


30. 


?757. 


95.3 


36. 


?Afl9. 


9a.9 


*2> 


a^i6. 


90. <• 


M. 


?s<.i. 


87. H 


5*. 


?<.59. 


85.0 


Mi 


2TM. 


82.0 


m. 


??86. 


T9.0 


m. 


??0l. 


T6.1 


7i. 


?1I6. 


73.1 


•». 


2033. 


70.2 


to. 


l-JftB. 


67.3 


96. 


IA66. 


f.<..5 



PER 1 001 I WO 


HEIGHT 


: <:rr 


U><WEl5H»Ef) 








HIFAS 


y'^^^H 


KSIGx* 


Y<;f G«i» 


oviBii'^ce 


T-<f T« 


?J.<.IO 


► 0.0 


13.110 




0.0 


n.o 


0.0 


0.0 


0.0 


?3.<.3? 


t 0.00* 


13.136 




0.00* 


0.210 


0.222 


0.012 


15.0 


?3. !.''<• 


0.008 


13.17* 




0.008 


0.*25 


0.*3I 


0.0*8 


1*.8 


23.<.'.9 


- 0.012 


13.223 




0.01? 


0.6*3 


0.626 


0.10* 


l*.l 


?3.<.50 


0.016 


13.281 




0.015 


0.8*7 


o.«oe 


0.169 


13.2 


yS.'.^T 


• 0.022 


13.*20 




0.021 


1.183 


I.II7 


0.23* 


9.5 


23.50* 


0.027 


13.600 




0.025 


l.**3 


1 . 156 


0.190 


5.2 


23.5SS - 


, 0.031 


13.803 




0.030 


1.629 


1.553 


0.099 


2.1 


23. ^.11 


>■ 0.03* 


I3.9T7 




0.033 


1.785 


1.702 


0.0** 


0.8 


23.<>75 < 


0.03T 


16.IIT 




0.0 3*. 


1.900 


1.828 


-0.006 


-0.1 


23.7*8 < 


0.0*0 


1*.2*T 




0.038 


2.013 


1.923 


•0.087 


-1.2 


23.778 


• 0.0*2 


Ifc.S*" 




0.0*1 


2.09J 


2.012 


-0.139 


-1.8 


23.777 . 


0.0** 


l*.*l* 




0.0*3 


2.1** 


2.073 


-0.1 3* 


-1.7 


23.909 


0.0*6 


1*.509 




0.0** 


2.212 


2.IS6 


0.005 


0.1 


23.835 - 


r 0.0*8 


1*.587 




0.0*6 


2.24S 


?.l*l 


0.09k 


1.0 


23.8*1 < 


0.050 


l*.688 




0.0*7 


2.32« 


?.I*S 


0.212 


2.2 


23.9** 


► 0.052 


1*.82* 




0.0*8 


2.353 


2.IT5 


0.1*3 


1.5 


23.888 


. 0.05* 


l*.9*5 




0.050 


2.373 


2.187 


0.15? 


1.6 


23.9*2 . 


0.056 


15.026 


* 


0.051 


2.398 


2.10* 


0,185 


1.8 



Soulier 5T»TI0H ISl 20 - ILL B»LO*IN 
TIME i*E|GHT i BEM»IN!t*0 *«EAM CONC. 



€ 



fy>) 



to 

O 



0. 
3. 
6. 
9. 
12. 

Is. 

I«. 
9«« 
42. 

*»» 

be. 

66. 
T2. 

TB. 
«<t. 
90. 
■96. 



8057.893 

7033.979 

61*5.295 

53'.2.606 

*676.0n5 

3529.7B9 

2650.138 

1969.5** 

1*53.072 

10«.8.*^2 

7n5. 196 

578.572 

*?3.018 

317.720 

2?2.*28 

1^0.*flJ 

115.022 

82.687 

"^9.910 



loo.o 

87,3 

76.3 

66.6 

5B.0 

*3.8 

32.9 

2*.* 

18.0 

13.3 

9.8 

7.2 

5.2 

3.8 

2.8 

2.0 

1.* 

1.0 

0.7 





OERlOOt 1990 


MEIG><TI 5EC 


AM CONC. 


KMEAN 


YHEAN 


2.78* 


23.410 




0.0 


13.110 


L 0.0 


2.*31 


23.*27 




0.003 


13.123 


• 0.003 


2.123 


23.*33 




0.006 


13.1*9 


L 0.006 


1.853 


23.*31 




0.009 


13.IB7 


t 0.009 


1.616 


23.*27 




0.013 


13.23* 


: 0.012 


1.228 


23.**1 




0.021 


13.363 


L 0.019 


0.937 


23.*92 




0.029 


13.528 


L 0.027 


0.71* 


23.531 




0.038 


13.719 


I 0.0 36 


0.5*0 


23.586 




0.0*9 


13.883 


. 0.0*7 


0.*08 


23.653 




0.060 


l*.001 


. 0.059 


0.309 


23.708 




0.075 


I*.lfl9 


> 0.07? 


0.235 


23.697 




0.090 


1*.183 < 


. 0.089 


0.I7H 


23.687 




0.108 


l*.2*l < 


. 0.109 


0.135 


23.698 




0.130 


1*.302 i 


. 0.133 


0.101 


23.669 




0.15* 


l*.338 < 


. 0.150 


0.0 76 


?3.h31 




0.183 


1*.365 < 


. 0.189 


0.057 


23.567 




0.21* 


l*.*27 . 


. 0.225 


0.0*2 


23.53* 




0.251 


1*.*7T < 


. 0.26* 


0.0 32 


^^.'^^? 




0.295 


1*.527 . 


. 0.309 



5ULPMUB OIOXlOE WEIGHTED 
XSIGX* YSIGMA COVAfflAsrCE 



0.0 
0.2IT 

fl.**2 
0.666 
O.STT 
1.22* 

i.*9e 

1.701 
1.B59 
1.9T1 
2.092 
2.16T 
2.220 
2.2T2 
2.29* 
2.315 
2.291 
?.285 
?,2«* 



0.0 


0.0 


0.229 


0.01* 


0.**2 


0.05T 


0.637 


0.120 


0.820 


0.18T 


1.122 


0.220 


1.176 


0.1*1 


1.611 


0.033 


1.785 


-0.0T5 


1.931 


-0.200 


2.032 


-0.3*1 


2.1*5 


-O.tSl 


2.252 


-0.*** 


z.-iyo 


-0.126 


2.392 


-0.01? 


2.397 


-0.015 


2.*17 


-O.20I 


2.*02 


-0.231 


2.388 


-0.186 



theta 



0.0 
16.8 
16.3 
15.2 
13.6 
8.3 
3.6 
0.6 
-1.2 
-2.9 
-*.5 
-5.9 
-5.1 
-1.* 
-0.1 
-0.2 
•I.? 
-2.6 
-2.0 



O 

O 

o 



iOURCF 


•;T*TlOM is: 


20 - III h»lu«im 


°E9|00! 1980 


HFIOHT! %rc 


SIILOMITE 


UFIOHTED 






TIME 


KEIGMI 


» PEM«ININ0 


ME4H CONC. 


XMF*N 


r'<E«M 


XSIGM* 


Y5!GM« 


cov«nt«'4CE 


TMET* 


«. 


Snft.fcOS 


100.0 


0.175 


23.*I0 


> 0.0 


3.110 


t 0.0 


0.0 


0.0 


n.o 


0.0 


3. 


7*.9.B*9 


i'il.O 


0.266 


23.*39 


> 0.007 


1.138 


I 0.008 


'>.l<»8 


0.710 


o.oi? 


15.6 


6* 


9q 1.590 


195.8 


0.3*3 


23. ".62 


r O.012 


3.179 


. 0.011 


1.393 


1.60I 


0.0*3 


15.6 


9. 


nse.2?6 


zzs.r 


O.*00 


23.679 


* 0.017 


3.226 


. 0.017 


0.588 


0.577 


0.087 


l».2 


12* 


l?TT,bT5 


^^^.3 


O.biil 


23.<.92 


>■ 0.022 


3.27? 


. 0.021 


0.773 


0.7*0 


0.131 


»?.» 


IS. 


IfcnS.SeO 


27B.1 


0.<.90 


23.51.0 


► 0.029 


3.356 . 


. 0.027 


1.081 


1.005 


0.1*9 


7,3 


2*. 


i<ii<).j<.r 


eao.3 


0.502 


23.596 


> 0.035 


3.<.52 


0.031 


1.331 


1.7*5 


0.118 


I.* 


3«. 


|]<i6.6<tH 


267.9 


0.'.92 


23.6*0 


► 0.0*1 


3.538 < 


. 0.0*0 


1.508 


1,*70 


0.058 


.|«3 


3«. 


I2-4B.096 


?«.<•. 5 


0.<>60 


23.690 < 


F 0.04' 


3.601 • 


. 0,0*6 


1.6*2 


1.631 


-0.010 


-«'.# 


*«. 


1116. ^55 


220. <• 


0.'»27 


23.759 


► 0.05J 


3.6*5 


. 0,05? 


1.762 


1.75* 


-0.0*9 


iO.9 


«». 


lOnl.JS^ 


197.7 


0.396 


23.831 


r 0.060 


3.726 - 


. 0,059 


1.B97 


1.873 


-0.121 


-1.9 


9«< 


a9<t.i<tt> 


176.6 


0.36* 


23.850 


t 0.066 


3.779 < 


. 0.057 


1.971 


?.no7 


-0.177 


-2.6 


M. 


itt.s?* 


153.5 


0.328 


23.850 


• 0.073 


3. BOB 


. 0,077 


2.039 


2.1*2 


-0.123 


-1.7 


66. 


6AZ.aS3 


130.9 


0.290 


23.803 


t O.OBO 


3.825 . 


O.OB" 


7.068 


2.270 


0.0*3 


0.6 


T2* 


5';3.3|9 


109.3 


0.251 


23.687 


i 0.0B7 


3.785 - 


. 0.100 


7.0*1 


2.3*3 


-0,051 


-0.7 


rt* 


*«>1.735 


91.2 


0.218 


23.669 


► 0.098 


3.759 - 


. 0.109 


7.100 


2.3*9 


-0.1*6 


-1.9 


M> 


3fl2.63' 


75.6 


o.ise 


23.640 


t 0.106 


3.769 . 


. 0.119 


7.079 


2.327 


-0.310 


-».l 


9Q* 


31^.649 


62.7 


0.163 


23.595 


► 0.115 


3.807 


. 0.128 


7.0*6 


2.786 


-0.217 


-3.0 


9*. 


2^9. 5ni 


53.2 


O.lk* 


23.617 


t 0.125 


3.83* 


. 0.137 


2.053 


2.757 


-0,037 


-0,* 



^>, 



n 

o 

n 
o 
ci 

(.r 

n 



n 

n 

o 

u 
u 

u 



SOUOCF STATION IS! en - ill AALOUtN 
T|Mf. KEtGMI « 0EM»IX|I»40 



0. 


738. 


100.0 


3. 


738. 


loo.o 


6. 


738. 


100.0 


9. 


738. 


loo.o 


1?. 


736. 


tJ-J.? 


18. 


fc8T. 


93.1 


?*. 


s^a. 


88.3 


30 • 


<.3*. 


H5.9 


36. 


Ml. 


82.8 


*2. 


S86. 


79.1. 


<>8. 


S70. 


77.2 


5*. 


S52. 


T*.8 


60. 


-^28. 


71.5 


M. 


Ji-JJ. 


66.8 


It. 


ftSa. 


6?.l 


78. 


(.^o. 


S6.9 


e<>. 


*03. 


5*. 6 


90. 


183. 


SI. 9 


<*6. 


170. 


50.1 



ocaiooi J«1 


HEIGHT 


t <;rr 


UMWElfiHTEr) 








«Mt»N 


T>»E»S 


J5IGX* 


y«;!6M» 


COViOIINCE 


TntT« 


23.<.I0 < 


O.D 


ii.no 




0.0 


0.0 


0-0 


n.o 


0.0 


23. ".0? < 


: 0.009 


l?.9?5 




0.012 


0.239 


0.336 


0.010 


9,6 


23,3R3 < 


' o.on 


12.858 




0.02* 


O.iflO 


8.AM 


0.019 


9.6 


PS.IfcS ! 


o.o?s 


12.759 




0.035 


n.7i6 


It.O*! 


0.099 


11.0 


?3.?'*2 < 


0.03<t 


12.687 




0.0*5 


0.932 


1.209 


0.186 


12.1 


?3.t63 < 


0.0<>9 


12.7*9 




0.059 


1.287 


1.55* 


0.265 


9.1 


23.01S < 


0.063 


12.900 




0.070 


1.576 


I.T95 


0.223 


5.1 


aa.«89 . 


. 0.072 


13.107 




0.080 


1.815 


2,019 


O.I 75 


3.0 


P2.005 < 


0.081 


13.3*7 




o.oss 


2.01* 


2.163 


0.0*? 


0.6 


??.75? . 


. 0.089 


13.582 




0.092 


2.159 


2.236 


-0.0*8 


-0.6 


ja.^Bl ! 


0.097 


11.828 




0.096 


2.305 


2.293 


-0.303 


-3.3 


?2.6I0 < 


0.103 


1*.078 




0.098 


2.*25 


2.30T 


-0.528 


-5.1 


??.5I6 . 


■ 0.109 


1*.319 




0.101 


2.*9* 


2.119 


-0.S69 


-5.2 


a?.S38 . 


' 0.116 


1*.*68 




0.105 


2.566 


2.325 


-0.*32 


-3.8 


3?.'iS0 < 


0.123 


1<..578 




0.108 


2.6?* 


2.309 


-0.903 


-7.5 


?2.5I<. . 


■ 0.1 ?8 


1*.775 




0.106 


2.631 


2.171 


-0.691 


-5.7 


??.537 . 


. o.n* 


15.0J7 




0.107 


2.689 


2.1** 


-0.T36 


-5.8 


?a.s7o < 


' 0.116 


15.270 




O.lll 


2.670 


2.1 TO 


-0.713 


-5.7 


??.6<»3 < 


0.1*3 


15.531 




0.113 


2.7*5 


2.182 


-0.816 


-6.2 



SOUBCF 5T*T10N 15! 20 - ILL B4LDWIN 
TIME HElOMT t OEMalNlNG «E»M COSC 




0. 

3. 

6. 

9. 
12. 
18. 
2*. 
30. 
36. 
*2. 
*8. 
5*. 
60. 
66. 
72. 
78. 
8*. 
90. 
96. 



3019. 9f,0 
263*.*69 
2280.95* 
197*. 881 
17(13. «fl* 
1198.370 

e';6,*20 

626.992 

*55.5o5 

329.'*** 

2*3.135 

179.251 

110. *23 

91.953 

6*. 21 2 

*3.B61 

11.937 

23.390 

l7.2** 



100.0 

86.7 

75.0 

65.0 

56.0 

39.* 

28.2 

20.6 

15.0 

10.9 

8.0 

5.9 

*.3 

3.0 

2.1 

1.* 

1.1 

0.8 

0.6 





OEBIODI J«N 


HE I GMT 1 5FC 


»M cose. 


«ME»M 


Y«<E»N 


*.I19 


23.*10 


. 0.0 


13.110 


>. 0.0 


3.570 


23. *0* 


' 0.005 


12.975 


L 0.007 


3.091 


23.187 


t 0.010 


12.836 


L 0,013 


2.676 


23.353 


> 0.016 


12.753 


L 0.021 


2.315 


23.301 


; 0.022 


12.672 


L 0,029 


1.7** 


23.179 


r 0.017 


12.711 


y, 0.04* 


1.31* 


23.052 


• 0.05* 


12.855 


L 0.051 


. 989 


22.955 


> 0.07* 


13.071 


L 0.080 


0.7*6 


22.89* 


i 0.097 


13.32* . 


. 0.102 


0.563 


22.878 


: 0.122 


13.57 7 - 


. 0.126 


0.*27 


22.805 


r 0.152 


11.880 ', 


. 0.150 


0.325 


22.73* 


■ 0.186 


l*.|90 . 


. 0.17* 


0.2*7 


22.673 


' 0.225 


1*.*72 < 


0.201 


0.187 


22.7*6 


. 0.275 


1*.635 . 


. 0.2*0 


0.1*0 


22.790 


; 0.335 


1*.703 < 


0.279 


0.10* 


22.815 


<■ 0.*07 


1*.816 < 


0.108 


0.079 


22.861 


0.*88 


15.036 < 


0.155 


0.061 


22.891 


. 0.568 


15. 2B* < 


. 0.*26 


0.0*7 


22.965 


0.668 


15.^12 < 


0.510 



SULPHUR 01 Oil or WEIr>HTED 
ISIGM* Y5IGM* COVAAIINCE 



0.0 

0.238 

0.*76 

0.709 

0.920 

1.280 

1.588 

1.8*8 

7.063 

2.212 

2.169 

2.*95 

2.568 

2.631 

2.682 

?.69f> 

2 . 756 

7, 7*7 

2.77S 



0.0 


0.0 


0.33* 


0.010 


0.6*2 


0.039 


0.925 


0.095 


1.182 


0.169 


1.522 


0,236 


1.776 


0.203 


2.012 


0.107 


2.1S3 


-0.050 


2.283 


-0.167 


2.119 


-0.391 


2.133 


-0.565 


2.121 


-0.536 


2.30* 


-0.236 


2.233 


-0.520 


2,0*1 


-0.059 


2.005 


-0.118 


2.058 


-0.222 


2.119 


-0.*23 



THCT« 



0.0 

9.7 

9.7 

10.8 

11.3 

8.2 

*.6 

I.B 

-0.7 

-2.0 

-*.0 

-5.2 

-*.6 

-\.t 

-*.l 



-0 
-0 
-1 
-3 



SOOBCE ST»TtOi^ ISJ iO - ILL BALDHIM 
TIME MEIGHT % REMAINIMO MEAM CONC. 



PERIOUI J»V 
KXE*N 



HEIGHTJ S^^C 

rHE«s 



It. 

*i* 
*•« 

Si. 
»9. 
»»• 

tm« 

96* 



1 no. 966 
?|3.5«J 
2TT.910 
253.019 
PS5.690 
?*2.75« 
2?T.559 
205.'?0 

165. 9|5 
151.0*8 

nl.^l' 

108.053 

n6.T8a 

6T.460 
59.763 
53.»65 



too.o 

128.3 
151.5 
168. 6 

181.3 

17?. 1 

|(>1.« 

l'.5.9 

131.0 

117.6 

107.1 

93.0 

76.6 

61.5 

• 7.8 

<>2.* 

37.9 

32.9 



0.191 

0.2<.5 

0.32? 
0.3«>* 
0.372 
0.372 
. 359 
0.337 
0.315 
0.291 
0.27fc 
0.2<>9 
0.219 
0.199 
0.161 
O.I^H 
O.lbO 
0.125 



23.410 
23.<>ll 
23.401 
23.369 
23.317 
23.215 
23.137 
23.10* 
23.076 
23.099 
23.061 
22.997 
23. nl* 
23.176 
23.258 
23.119 
23.160 
23.15* 
23.176 



O.B 

0.017 

0.032 

0.0*5 

0.056 

0.079 

0.102 

0.123 

0.150 

0.173 

0.198 

0.220 

0.2*9 

0.287 

0.331 

0.380 

0.*12 

0.*33 

0.*51 



13.110 
12.975 
12.853 
12.736 
12.625 
12.553 
12.61* 
12.781 
13.029 
13.330 
13.766 
I*. 183 
I*. 612 
I*. 9*1 
I*. 988 
15.0*1 
15.196 
15.365 
15.556 



0.0 

0.0?* 

0.0*2 

0.057 

0.071 

0.091 

D.I 10 

0.131 

0.153 

0.1 73 

0.187 

0.195 

0.207 

0.227 

0.2*2 

(1.2** 

0.255 

0.2'Jl 

0.301 



SilLPM»TE WEIGHTED 
X<;iG'4t 



o.o 

0.23* 
0.*65 
0.689 
0.897 
1.260 
1.593 
1.891 
2.157 
2.3*8 
2.5*9 
2.703 
2.857 
?.98? 
1.081 
1.119 
3.187 
3.16* 
1.085 



rsitt"* 


COVARUSCE 


O.n 


1.0 


0.127 


n.OlO 


0.621 


0.038 


0.886 


0.090 


1.131 


0.150 


l.*57 


0.19T 


l.7|8 


0.158 


1.979 


-0.036 


2.199 


-1.251 


2.3*7 


-0.*86 


2.*12 


-0.791 


2.198 


-0.89* 


2.371 


-0.72* 


2.365 


-0.153 


2.252 


0.000 


2.003 


0.7*1 


1.969 


0.60* 


2.058 


0.*0T 


2.05* 


0.18* 



TKET« 



«*• 

».» 

IttkT 

!«.• 

T.I 

S.ft 

•3.1 

-T.O 
-S.I 

9.0 
*.* 

i.3 
I.I 




SOU«CF "^TtllOK ISl 30 - lUL BALOwlM 
TIME XEIGHT « OEM4IN1H0 



0. 


7?0. 


100.0 


3. 


720. 


100.0 


6. 


7?0. 


100. 


9. 


720. 


100.0 


12. 


720. 


100.0 


18. 


720. 


100. 


2*. 


712. 


'io.t 


30* 


H91. 


96.0 


36. 


ft69. 


9a. 9 


•F. 


fc52. 


^O.fe 


I>8. 


ft27. 


H7.I 


5*. 


^0^. 


B3.6 


Ml. 


■^Sl. 


80.7 


(«. 


S62. 


78.1 


72. 


S*2. 


73.3 


78. 


SOT. 


70.* 


B*. 


*70. 


65.3 


90. 


'•4|. 


M.3 


96. 


<.0*. 


S6.I 



PEOIOO! »'>9 


Htlfi'* 


: 5rc 


ilMWEISHrEO 








XMEHN 


»te»N 


x^ir>«t 


Y«.IGM» 


COV»PI»«»CE 


T^ETi 


23. ".10 




0.0 


13.110 


0,0 


n.o 


0.0 


0.0 


0.0 


23. ".I* 




0.009 


13.1*9 


. 0.011 


1.215 


0.298 


-0.01 1 


-12.8 


23. ".12 




0,016 


13.20* 


0.02? 


il,*?5 


0.580 


-0.0*5 


-n,9 


23.403 




0.023 


n.26S 


. 0.031 


1.625 


0.838 


-0.101 


-J*. 6 


23.386 




0.010 


13.336 < 


0.0*0 


O.BH 


1.078 


-0.179 


-15.0 


23.337 




0.0*3 


13.*60 


. 0.056 


1,167 


l,»90 


-0.361 


-1*.8 


23.271 




0.05* 


13.601 


0.06^ 


l.**3 


1,801 


-o.*e* 


-13.1 


23.211 




0.061 


13.788 


. 0,078 


1,61* 


2.0*1 


-0.59? 


-12.8 


23.130 




0.065 


13,9*6 


. 0.0')7 


1.689 


2,2*7 


-0.781 


-15.3 


?J.073 




0.071 


1*.097 


0.094 


1.816 


?.*12 


-1.021 


-17.? 


23.023 




O.075 


1*.312 


. 0.101 


1.903 


2.523 


-1,33* 


-20.2 


23.026 




0.083 


I*,*e6 


► O.IOS 


2.026 


2.585 


-1,708 


-22.6 


22.99? 




0.090 


l*.67? 


t 0.109 


2.157 


2.638 


-2.105 


-2*.1 


22.951 




0.097 


|*.876 


L O.ll* 


2.296 


2.698 


-2.*?0 


-?*.7 


22.957 




0.103 


15.078 


> 0.119 


2.*02 


2.768 


-2.7*6 


-25.* 


22.967 




0.108 


15.201 


► 0.123 


2.*37 


2.775 


-?.79* 


-25.2 


23.010 




fl.ll 3 


15.227 


• 0.129 


?.*50 


2.798 


-2.83? 


-25.2 


23.069 




0.121 


15.2*3 


t 1,132 


2.53* 


2.77* 


-1.113 


-75.9 


22.910 




0.120 


15. 33* 


> 0.135 


2.*22 


2.715 


-2.825 


-?5.7 



50U9CE 5T»TI0M 1st 20 - ILL B»L0«1N 
TIX| WEIGHT « BEM«|N|NG ME 






0, 
3, 
6, 
9. 

12. 

18. 

2*. 

30. 

36. 

*2. 

*8. 

5*. 

60. 

66. 

72. 

7H. 

8*. 

90. 

96. 



1814,52* 

1574,61 I 

1355.8*3 

1169.353 

1013.403 

766.173 

576.063 

4?a.*55 

315.980 

217.086 

174.951 

129,262 

95.*10 

69.795 

51.442 

16.020 

25.*?3 

l8.J<J2 

12.816 



100,0 

85.8 

73.9 

63.7 

55.2 

41.8 

31.* 

23.* 

17.2 

12,9 

9.5 

7.0 

5.2 

3,8 

2.8 

2.0 

I.* 

1.0 

0.7 





PERIOOi «oR 


HEIGHi 


! Src 


«<J COMC. 


XHE**! 


y<*l^s 


2.5*8 


23.*10 




0.0 


3,1)0 ; 


. 0.0 


2.187 


23.41* 




0,005 


3,1*4 i 


0.007 


1.883 


23.415 




0.011 


3.186 


L 0.015 


1.624 


23.408 




0.018 


1,231 


I 0.021 


1.408 


23.397 




0,025 


3.272 


t 0.032 


1.064 


23.375 




0.0*1 


3.330 


L 0.051 


0.809 


23.331 




0.059 


3.*|7 


►. 0,072 


0.620 


23.269 




0.076 


1.558 


;. 0,096 


0.472 


23.192 




0.091 


3.690 < 


L 0.125 


0,36* 


23.157 




0.114 


3.791 


L 0.157 


0,279 


23.101 




0.138 


3.973 


> 0.195 


0.215 


23.075 




0,170 


*.088 


• 0.232 


0,154 


23.023 




0.212 


».230 < 


0.276 


0,12* 


22.906 




0.262 


*.*b8 < 


. 0.132 


0,095 


22.900 




0.327 


*.587 < 


0.402 


0.071 


22.934 




0.379 


*.6?7 


. 0.4fB 


0,054 


22,946 




0.466 


4.518 


I 0.566 


0,0*2 


23.053 




0.565 


4.416 


> 0.65) 


0.032 


22,875 




0.611 


4,1.69 


. 0.750 



5llLPML(D DIOXIDE WEIGHTFO 
«SI6M« V5lGMt COW»0|«NCE 



0.0 

0.21* 

0.422 

0.617 

0.805 

1.144 

1.407 

1.570 

1.615 

1,751 

1,823 

1,930 

2.068 

2.191 

?.3*3 

2.275 

2.348 

2.424 

2.1«6 



0.0 


o.« 


0.293 


-0.010 


0.562 


-0.038 


0.803 


-0.082 


1.026 


•0.1*1 


1.411 


-0.280 


1.720 


-0.364 


1.989 


-0.431 


2.^20 


-0.686 


2.421 


-0.9*9 


2.573 


-1.237 


2,6*0 


-1.579 


2,698 


-2.072 


2.77* 


-2.25* 


2.882 


-2.710 


2.869 


-2.735 


2.85* 


•2.866 


2.790 


-1.1*0 


2.713 


-2,645 



TiffT* 



0.0 
-11.8 
-11.9 
-12.1 
-12.2 
-12.1 
-10,4 
-9.9 
-14.7 
-17.2 
-20. » 
-21.0 
-25.9 
-25.1 
->6.3 
-?7.9 
-27.5 
-28,2 
-?9,0 



50U»CF «1T/ITI0N IS! ?0 - ILL B»LOHIM 
TIME HEIGHT t 0EW4IHING ME«M COVC. 



0. 

3. 

6. 

9. 
\Z. 
IS. 
!•'. 
30. 
36. 
*2. 
48. 
5<t. 
60. 
«6. 
72. 
78. 
80. 
90. 
9b. 



ll*.6SH 
1S3.179 

ina.us6 

2}S.4]0 
215. TfcB 
2S6.80* 
?SB.a55 
2<.6.7?3 
2?*.S4* 
2n7.2<,5 
lfll.*<.9 
I'.O.SPI 
laa.ogS 
119.208 
99,2';6 
BO. 252 

f.2.*82 

"iO.Tsa 



100.0 

133.6 

I6<>.0 

187.9 

209.6 

22*. 

225.8 

215.2 

H5.8 

180.8 

158.3 

1*0.0 

12*. 2 

I0*.0 

86.6 

70.0 

60.1 

5*.S 

**.3 



1 


PE«1 


30! *>» 


HEir.HT 


t 5FC 


«M cove. 


XMFAN 


YXf av 


o.isq 


23.*I0 




0.0 


3.110 




0.0 


0.213 


23.*15 




O.017 


1.137 




0.02? 


0.261 


23.*ia 




0.030 


3.172 




0.03R 


0.29<» 


23.*20 




0.0*1 


3.I9R 




0.05? 


0.327 


23.*25 




0.051 


3.20* 




0.063 


0.357 


23.*** 




0.0 70 


3.161 




0.081 


0.36* 


23.*JI 




0.08* 


3.198 




0.102 


0.357 


23.*19 




0.096 


3.258 




0.121 


0.336 


23.*02 




0.10* 


3.255 




0.1*? 


0.318 


23.*67 




0.120 


3.217 




0.16? 


0.289 


23.**3 




0.1 3S 


3.23* 




O.I 8* 


0.267 


23. *0* 




0.1*8 


1.181 




0.193 


0.2*5 


23.37* 




0.169 


3.291 




0.209 


0.212 


23.298 




0.196 


3.*37 




0.231 


0.183 


2 3.3*2 




0.225 


1.250 




0.251 


0.158 


23.*32 




0.2*5 


3.060 




0.261 


0.1*7 


23.656 




0.275 


2.856 




0.263 


0.1*2 


23.797 




0.288 


2.8*6 




0.277 


0.126 


23.597 




0.2*9 


2.867 




0.30? 



SULPHITE WflTiHTEO 
iSir,M« v<;lG«« 



0.0 

0.211 

0.*15 

0.605 

'J.7SI 

1.122 

1.1*7 

1.506 

I.550 

1.733 

I. 815 

1.877 

2.011 

2.1*1 

?.?<'> 

2.195 

2.281 

2.280 

I.77B 



0.0 

0.?78 

0.528 

0.760 

0.970 

1.127 

1.6*1 

1.936 

2.132 

2.137 

2.»79 

2.*»3 

2.*92 

2.5*5 

2.505 

2.353 

2.18ft 

2.188 

2.1*9 



covtoitvcr 



n.O 
-O.OOT 
-n.026 
-0.06* 
-0.117 

-n.?6i 

-1.33* 
-0.*05 
-0.666 
-0.808 
-1.123 
-l.*6* 
-7.01! 
-2.095 
-2.330 
-2.292 
-2.3*1 
•?.*8* 
-1.689 



THET« 



0.0 
-9.1 
-8.T 
-9.9 

-in. 6 

-II. 7 
-10.* 
-10. I 
-15.3 
-15.1 
-18.8 
-22.6 
-26.* 
-2*. 6 
-?*.9 
-?5.t 
-?*.2 
-25,5 
-28.1 




SOUWCF '■,T»TI0N IS! ^0 - ILL BlLflKlvi 
II>4E «IEIGH? « >?EM/l|VINO 



0. 


7*fc. 


100.0 


3. 


7*4. 


100. 


6. 


7<«*. 


100.0 


<). 


TlXt. 


100.0 


1^. 


?<•<». 


100.0 


18. 


7*<.. 


loo.o 


?*. 


T**. 


100.0 


30. 


7'.'.. 


loo.o 


36. 


74<(, 


100.0 


*^. 


7**. 


100.0 


oe. 


7**. 


100.0 


5<.. 


yi'O. 


99. S 


60. 


7^^. 


98.(1 


66. 


T?8. 


97.8 


72, 


719. 


96.6 


78. 


713. 


95.8 


Bi.. 


706. 


9l>.9 


90. 


m99. 


9<i.0 


96. 


^94. 


93.3 



PEUIOO: JJL 


•lElG-fr: ^^c 


l)NWE15HfED 








xxt*»< 


f^H^ 


HSIG"* 


Y5IGM* 


COViOI^VCE 


T'<f T» 


?3.<.10 


> 0.0 


13.110 


. 0.0 


0.0 


0.0 


0.0 


0.0 


?3.<.36 


' o.oo; 


13.1«0 


0.00* 


0.131 


0.106 


0.00? 


7.* 


?3.'.61 


; 0.010 


13.^53 


. O.OOfl 


O.?60 


o,?n 


0.008 


7.1 


?3.ii9l 


► 0.0l<. 


13.389 


0.01? 


1.386 


0.315 


0.018 


6.9 


?3.5?2 


r 0.019 


13.*05 < 


. 0.015 


0.5O7 


0.*1* 


0.030 


6.7 


23.598 


' 0.037 


13.551 


. 0.02? 


0.7 30 


0.59* 


0.057 


6.1 


?3.66* 


• 0.03* 


13.689 


. 0.0?8 


0.934 


0.T57 


0.07* 


*.9 


?3.750 


► 0.0*1 


13.816 


. 0.033 


1.117 


0.906 


0.065 


3.0 


a3.fl*3 


« 0.0*7 


13.930 


0.038 


!.?75 


I.O*? 


0.025 


0.9 


a3.'*3'» 


► 0,05? 


l*.03O 


. 0.0*1 


l.ftns 


1.165 


-0.03S 


-1.0 


2'.. 0*1 


> 0.056 


l*.lia 


0.0*7 


1.51* 


l.?90 


-O.IIO 


-2.8 


2'..|3<. 


«• 0.059 


1*.?19 


0.050 


1.599 


1.357 


-0.131 


-2.9 


a*.?2'. 


' 0.061 


l*.33? 


0.05? 


1.655 


l.*05 
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0.890 


?!.? 


%i. 


8.291 


1.7 
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1. INTRODUCTION 

This ADDENDUM to the Final Report documents the work which was 
done as. a result of some of the findings on the trajectory statistics 
for 1978-80. Those trajectories were generated using the OME Trajectory 
Model based on wind fields derived from surface pressure data. The 
means of the trajectory end-points for some of the months in 1978-80 
suggested rather unexpected movement patterns. For this reason, the 
Ministry decided that parts of the code should be examined for possible 
inconsistencies, and also that comparisons should be made between sur- 
face and 850 mb trajectories. Due to data availability, the months of 
January and July of 1982 were chosen for which the trajectories and 
their statistics would be generated. This ADDENDUM documents this work 
in detail. 
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2. EXAMINATION OF THE ONTARIO MINISTRY OF THE ENVIRONMENT 
LAGRAN6IAN LONG RANGE TRAJECTORY MODEL 

The examination of the Long Range Lagrangian Trajectory Model 
developed by the Air Resources Branch of the Ontario Ministry of the 
Environment was restricted to the FORTRAN Source Code responsible for 
calculating the trajectories. This segment of the code is comprised of 
the following subroutines: 

TRAJ, 

NEWXY, 

VELO, 

and CMCLL. 

The purpose of the subroutine TRAJ is to calculate the tra- 
jectories. The program processing logic is straightforward. The sub- 
routine makes extensive use of do-loops to process more than one trajec- 
tory at a time. To simplify the following discussion, the generation of 
only one trajectory will be discussed. 

Although the OME -modified CMC polar stereographic coordinate 
system is used to specify the location of source points and surface 
pressure data, it is not used internally in the Lagrangian Trajectory 
Model. Instead, these locations are transformed to corresponding lo- 
cations on a square of unit length on each side. This unit square is 



I 
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further subdivided into smaller squares. This square grid is used in 
the model because it avoids potential problems with the orthogonal poly- 
nomial scheme that is used to model the surface pressure field. These 
problems, which were discussed in Sykes and Hatton (1976), and Polavarpu 
(1980), are: 

1) domain difficulties because polynomials tend to increase or 
decrease without bound as the coordinates increase, and 

2) locally inaccurate fit resulting from sparse data near the 
boundaries. 

As well, to reduce the effect of the latter problem, the grid 
is increased by 5% on all sides. The region of the polar stereographic 
grid that is mapped into the square grid depends on the availability of 
surface pressure and can change from time to time. The extremities of 
the locations of the data determine the region of the polar stereo- 
graphic grid that is mapped into the square. The largest and smallest 
x-coordi nates and the largest and smallest y-coordi nates of the avail- 
able locations of surface pressure data determine the rectangular region 
in the polar stereographic grid that is mapped into the square grid 
whose coordinates range in value from zero to one. 

The process of generating the next point in the trajectory 
begins with the calculation of the surface geostrophic wind vector u^. 
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from the surface pressure field, at the latest point p^in the trajec- 
tory. This wind is used to calculate a "first-guess" position for the 

new point p, = p + u At, where the superscript indicates the guess 

'^ loo 

number, or iteration number. If the first-guess position falls outside 
of the normalized grid then the trajectory stops at the latest point and 
no further processing for the trajectory will occur. If, however, the 
first-guess position falls within the grid, then a surface geostrophic 
wind ti, is calculated for this new point. A second-guess for the next 
point in the trajectory is calculated assuming constant acceleration 
between the initial and the first-guess points: 



Pi " Pq ^ 2 ^^ 



Once again, the new point p^ is examined as to whether it lies within 
the grid. If it does not, then processing of the trajectory is stopped. 
Otherwise, processing continues and the convergence criterion is 
checked. Let the components of the position pj be written as: 



k+ k -» . k -» 
Pi " ^1^ ^1^ 



If px^ - \\<C 



and 



f >! - ^yi I 



<c 
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then processing stops and x, is the new endpoint. If not then the 
iterative process begins and p,is calculated as follows: 



2* ^ * 
Pi = Po * ( 2 ' ^^ 



In general, for k > 2, 

'Pl - % * ' 2^' " • 

Each time p, is calculated, it is checked to determine whether it lies 

within the grid. If it does then the convergence criterion is tested, 
that is, if 



and 



fxi - ^'-'\ I < C 



,k (k-1) I , r 



k-* 
then the iterative process stops and Pj^ is the new point in the trajec- 
tory. To prevent excessive and futile computation the number of iter- 
ations is restricted so that k < k^^ . If a solution is not found 

max 

within k , iterations (remember p, - the first guess position is not 
max-i i- 

calculated in the iterative procedure since p^^ = p^ + u^ At), then the 
trajectory ends at the previous endpoint and therefore doesn't survive 
for the full expected life-time. It has been determined by B. Ley and 



i)'' 



^ , ^ oncord Scientific CorooratiCi 



- 6 



G. Ellenton of OME that requiring convergence within 8 iterations 

(k = 9) and setting C - 0.001 works quite well. The value of C is 

^ max 

in distance units of the square grid and is approximately equal to 

0.381 kiji on the earth's surface. 

The subroutine NEWXY is written to find the candidate next 
trajectory point. In a straightforward way it accepts the wind speed 
and direction from which it determines the x and y wind velocity compo- 
nents, multiplies these components by time and then adds them to the 
appropriate components of the last trajectory point and thereby creates 
the new point. 

The purpose of the subroutine VELO is to calculate the surface 
geostrophic wind at a point and calculate a mean wind velocity by aver- 
aging the wind velocities of the point and a previous point. The sub- 
routine also contains code to back wind direction and alter wind speed 
but this is no longer used. It is recommended that it be deleted 
entirely from the code. It would also be advisable to remove the 
section of code that calculates the mean wind velocity and set it up in 
another separate subroutine. 

The equations used in VELO to calculate the surface geostro- 
phic wind are developed as follows. Start with the equation for the 
geostrophic wind: 



s; 
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The subscript s serves to indicate that the coordinate system resides on 
the surface of the earth. The variable p is the surface pressure, p the 
density of air at the surface, and f is the Coriolis parameter. The 
unit vector I is taken to point in the upward direction. A thorough 
derivation of the geostrophic wind equation (pressure-gradient form) is 
given in Holton (1979). 

To transform the equation into the polar stereographic coor- 
dinate system note that: 



d = md 



where d is the true earth distance and d is the corresponding distance 
in the polar stereographic coordinate system. The map factor m is given 



1 + sin <t>Q 
^ ^ 1 + sin 4. 



where ♦ is the reference latitude where the projection is "true" (^^ = 
60 • in the CMC Polar Stereographic Projection) and * is the latitude on 
which the distance lies. 
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Now, if &x and by represent distances on the polar stereo- 
graphic projection and dx , and dy , represent the corresponding true 
distances on the sphere of the earth, then from the discussion on map 
factor it is clear that 



dx = m dx 



and 



By = m ay^ 



With these two equations, it is possible to write 



at ~ m at 



at m at 



In vector form these two equations can be written as 



\-l - (2-1) 



Defining i and j to be unit vectors of the polar stereographic 
coordinate system corresponding to the x- and y-axis, respectively, the 
pressure gradient is 



i^J cor. 



;o!a Scientifi: Corporaiior 



9 - 






This can be rewritten as 



^ 5x &x^ ay ay^ -^ * 



so that 



^ m s^ 



> 



or alternatively, 

^^p = m ^p , ... (2.2) 

With the equation relating the wind velocity on the earth's 
surface (equation 2.1) to that on the polar stereographic projection and 
the equation relating the pressure gradient on the earth's surface to 
that on the polar stereographic projection (equation 2.2), one can do 
two simple substitutions in the geostrophic wind equation and obtain its 
counterpart for the polar stereographic projection: 

V = t X ^ ^p . 

The length unit used in this equation is the polar-stereographic unit, 
which equals the conventional metre only at the reference latitude. The 
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grid used internally is the square grid of unit length which is con- 
structed from the polar stereographic grid. The transformation from the 
polar stereographic grid to the square grid is done with scaling factors 

C and C . The dimensions of the polar stereographic grid used varies 
X y 

with data availability. Define w to be the dimension of the grid in 
the x-direction in CMC length units, and w to be the dimension of the 
grid in the y-direction in CMC length units. The scaling factors are 
thus given by 



C = w X 381,000 

X X 



and 



C = w„ X 381,000 



The dimension of C and C is length and the unit is the metre because 

X y 

one CMC grid length is equal to 381,000 polar-stereographic units. 
Choosing a and p as the coordinates of the square grid, the suitable 
transformation from the CMC polar stereographic is described by 



X 

" -T- 

X 



and 



C 

y 



I 
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Noting from this that 



dx = C^da 



mid 



ay = C^dp, 



the geostrophic wind in component form in the OME modified CMC Grid from 
computations in the square grid is 



.EL. i_a£ ... (2.3) 

^ " fp C ba 



_ .mi 1 &£ ... (2.4) 

The units for these wind velocity components are the polar-stereographic 
unit per second. 

The geostrophic wind components for the surface calculated by 
the subroutine VELO are those given by equations (2.3) and (2.4). VELO 
also carries calculations a couple of steps further. With the two com- 
ponents it calculates the wind direction in OME modified CMC grid. 
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= arctan (■^) 



and the wind speed. 



1 T 

s = /u + V 



Furthermore the wind speed is converted into the units of polar-stereo- 
graphic kilometres per hour in the following manner: 



r^ ^ ocnn S km 

s = /u + V X 3600 ^ X j^i^Q^ 



—K ^ 

or s = (/ u + V ) X 3.6 



This conversion to wind speed and direction serves no real purpose for 
when the subroutine NEWXY is used to calculate the next point in the 
trajectory it must use the u and v components of the wind. To do so it 
calculates 



^^ ' AtVr 



where s is in polar-stereographic kilometres per hour, &t is in hours, 
and C is 381 polar-stereographic kilometres so that ^6 is in the unit 
of length of the OME modified CMC grid. The displacements in the a and 
p coordinates are 



Aa = &d cos 9 




I 
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and 



AP = Ad sin 0, 
respectively. 

A noticeable reduction in computer time could be achieved by 
calculating the velocity components alone using units of length per hour 
rather than converting the wind velocity to and from a wind speed and 
direction, and changing the unit of wind speed as is presently done. 

In the subroutine VELO the density at the surface is assumed 
to be a constant 1.29 kg/n? . This is an approximation which affects the 
calculation of the geostrophic wind since density is one of the factors 
in the equation 

t^ = t X fp ^^p 

Now the surface for the purpose of this study has been taken to be mean 
sea level and it must be noted that very few land stations lie at this 
elevation so that the density of the air at mean sea level, if estimated 
by some algorithm, would be a fictitious value. This is the reasoning 
behind which the decision was made to let density be a constant in the 
mean sea level version of OME's Lagrangian Model. Nevertheless, one can 
by a simple calculation estimate the effects of this assumption. Let us 
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examine two meteorological extremes: a cold high and a warm low. 
Suppose the temperature to be -ZO'C and the surface pressure to be 1030 
mb in the cold high. This results in a density of 1.42 kgm" in dry 
air. In this case the constant density assumption reduces the geostro- 
phic wind speed by 9.2%. For the warm low, suppose the temperature is 
30- C and the pressure is 980 mb. The density of dry air is then 1.13 
kgm"-^. In this case the constant density assumptions increases the 
geostrophic wind speed by 14.2%. Clearly the effects of the constant 
density assumption on the calculated geostrophic wind value are not 
severe and the use of a constant density is justified. The value of 
1.26 kgm"^ for the density falls within the range of the two extremes 
given by the preceeding synoptic examples and is close to the density at 
mean sea level of 1.225 kgm'^ given by the ICAO Standard Atmosphere 
(where the pressure is 1013.25 mb and the temperature is 288.15'K). 

The determination of the latitude of a point in the OME- 
modified CMC polar stereographic grid and the subsequent calculation of 
the map factor is performed in the subroutine CMCLL. This subroutine 
also calculates the longitude of the point but this is not used in the 
model. The CMC polar-stereographic grid is constructed so that the 
origin coincides with the North Pole and the positive x axis is aligned 
with 10 -E longitude, so the negative y-axis is aligned with 80 "W. and 
the latitude where the projection is true is 60 °N. Let x^ and y^ 
denote the CMC coordinates. 
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The OME-modified coordinates, denoted x^ and y^, are given by 



Xr = X +26 

E c 



and 



yc^28 



Since the longitude of a point is not used in the model, the calculation 
of latitude will be discussed only. Given the OME-modified coordinates, 
one must transform first to the CMC coordinates. With the aid of Figure 
2.1, the calculation of latitude is straightforward. From the geometry 
of Figure 2.1, it follows that: 



^^" ^ = a(l + sin $^) 



~2 2 

where d = / x + y 



and a is the mean radius of the earth. From the geometry of isosceles 
triangles, 



4, = It - 2 {ti/2 - e) 



'tW* 



(b = 29 
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Noting that ^ = j 4< 



reveals ♦ = | — 26 



or 



/(xg - 26)^ ^ (y - 28)' 



^ = ^ - 2 arctan [/ a(l ^ sin ,) 1 



as the required result. Notice that the mean radius of the earth must 
be expressed in CMC length units, that is, a - 16.7. This is not 
exactly what is coded in CMCLL. 

Instead, what is given is 



/(x - 26)^+ (y - 28)^ 
= 4-2 arctan [/ ^(1 ^ sin I) ] ^ 



0' 



where E = 0.005 radian. The deviation from what the calculation ought 
to be is neglible. But it was interesting to learn that this small 
adjustment was added at some time in the past so that latitudes like 
49.99 ° were rounded to 50 *. This is poor programming style. Cosmetic 
adjustments to printed results should not be incorporated in utility 
subroutines and we recommend that this practice be avoided. 
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The calculation of the map factor by subroutine CMCLL is done 
using 

1 + sin <> 
""•" 1 + sin t 

In summary, examination of the FORTRAN source code for the 
Lagrangian Model revealed no errors in the translation of formulae nor 
in the algorithm responsible for generating the trajectories. Due to 
time limitations, an audit was not done of either the input or output 
segments of the code. We did not examine the code responsible for 
fitting the two dimensional orthogonal polynominals to the mean sea 
level pressure data because we received assurances that this code had 
been well studied in the past and therefore we need not concern our- 
selves with this. 

It may be said that documentation of the model has been rather 
thin in the past. With our analysis we hope that we have been able to 
illuminate some of the facts of the model that have been obscured from 
view, particularly, the assumption in the calculation of the man sea 
level geostrophic wind calculation of constant density. The reasons for 
using this assumption are sound and there is no reason for not docu- 
menting them. 

Aside from a few minor comments that have been made and a 
suggestion or two to make the code more efficient (and thus less expen- 
sive to run), no deficiencies have ben identified in the Ministry's 
Lagrangian Model. 
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Figure 2.1. Polar stereographic projection where 4" = 60"N, 
and a = 6.37 km is the mean radius of the earth. 
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I, MODIFICATION OF MODEL TO CALCULATE 850 mb TRAJECTORIES 

The modification of the Ministry's Long-Range Lagrangian Model 
to calculate 850 mb trajectories instead of surface trajectories is by 
no means complicated. To begin with, we replace the surface geostrophic 
with the isobaric coordinate form of the geostrophic wind equationi 

p f p » 

where the subscript p serves to remind that the gradient must be calcu- 
lated on a constant pressure surface, the wind velocity is for a point 
on that surface, and $ is a function of height z that is defined to be 
the work required to lift a unit mass to a height z above mean sea 
level. An equation for * is 



»(z) = /^ gdz 



The derivation of the isobaric coordinate form of the geostrophic 
equation is well known and will not be repeated here. Nevertheless, a 
discussion of the derivation can be found in Chapter 3 of Holton (1979). 
It is convenient to replace *{z) in the geostrophic wind equation by the 
geopotential height Z that is defined as 



9,. 
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_2 
where g = 9.80665 ms , the average global acceleration of gravity at 

mean sea level. With this substitution, the geostrophic wind equation 

in isobaric coordinates becomes 



v^ = t X 7^ ^ Z 
P T p 

Note that this equation is valid for any constant pressure surface and 
not just 850 mb. Fortunately, geopotential height at selected pressure 
levels are the values reduced from radiosonde ascent and collected for 
dissemination. So no further data preparation is required. 

In passing, we note that one of the touted advantages of the 
isobaric-coordinates form of the geostrophic wind equation over its 
constant altitude counterpart is that density no longer appears as an 
explicit parameter in the equation. The difficulty with the constant 
altitude form of the equation is that if the pressure gradient is the 
same at two points, the geostrophic wind will be different at the two 
points if the density is not the same. 

The only other changes to the model result from the limited 
frequency of observations and the sparse observational network. Radio- 
sondes are released only twice daily, at OOUT and 12UT, and the number 
of stations making upper air observations are far fewer than the niPiber 
of surface observing stations. The 850 mb version uses 12h time steps 
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between trajectory points rather than the 3 h and 6 h interval between 
points for the surface version. At 850 mb it is not feasible to fit an 
orthogonal polynomial surface to the small number of station values, so 
an objective analysis must be done to produce values of geopotential 
height on a grid and these gridded values are then used to fit the 
orthogonal polynomial surface. 

It should be noted that while the modification made to the 
model was done so that 850 mb trajectories could be generated, it is 
certainly possible to use the modified code as is for any of the other 
nine mandatory pressure levels whose geopotential heights are routinely 
reported twice daily around the world. So this version of the 
Lagrangian model that has been developed is not so much the 850 mb 
version of the model as it is the upper air version. 
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4, GENERATION OF SURFACE TRAJECTORIES 

Mean sea level trajectories were generated at three hour 
interv.als in the months of January and July of 1982 for three source 
locations: Sudbury, Ontario; Baldwin, Illinois; and a location in 
southwestern Pennsylvania. The pressure data were obtained from two 
sources. One was mean sea level pressure reported by surface observing 
stations. The other was mean sea level pressure reduced from 1000 mb 
geopotential heights at 28 points from an objectively analysed grid. As 
can be seen in Figure 4.1, these 28 points lie, for the most part, over 
large bodies of water where there are no surface weather stations. The 
limited extent of the gridded 850 mb geopotential heights has resulted 
in a restricted rectangular boundary for trajectory generation. This 
boundary constrains the x-coordinate to lie between 20 and 30, inclu- 
sive, and the y-coordinate to lie between 10 and 20, inclusive. For 
consistency, this same region-of-interest was used for the surface as 
well. To partially compensate for boundary effects on trajectory calcu- 
lation, a number of the additional 28 points of mean sea level pressure 
lie along y == 10 and x = 30. 

For a few time periods there are no pressure data available at 
all (for a variety of reasons). So a decision was made to halt a tra- 
jectory at the latest point if there was no data available at the next 
time step. Moreover, any trajectory that does not survive at least 24 h 
is not output to the data file. 
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The time-step for the first 18 h of travel -time is 3 h, and 
then 6 h thereafter until the trajectory is terminated at 96 h of 
elapsed time. 

Trajectories are intiated 8 times per day at 3 h intervals 
beginning at OOUT on the first day of each month. Since the data for 
each month include two more days worth of the following month, the last 
trajectory is initiated at 21Z of the first day of the following month 
(either 1 February or 1 August). 

The estimated trajectory - statistics for Sudbury are given in 
Table 4.1 for January and in Table 4.2 for July. The number of trajec- 
tories initiated (see the entry under the column marked WEIGHT for time 
Oh) differs between the two months because of lack of pressure data at 
certain times. 

Tables 4.3 and 4.4 report the estimated statistics calculated 
for the set elapsed times of trajectories initiated from Baldwin for the 
monthly periods of January and July, respectively. 

The statistics estimated from the trajectories initiated from 
the site in Pennsylvania are given in Table 4.5 for January and in Table 
4.6 for July. 



Concord Scientific Corporafior 



Table 4.1: Estimated trajectory-statistics at mean sea level for Sudbury, January 1982. 
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Table 4.2: Estimated trajectory-statistics at mean sea level for Sudbury, July 1982. 
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Table 4.3: Estimated trajectory-statistics at mean sea level for Baldwin, January 1983. 
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Table 4.4: Estimated trajectory-statistics at mean sea level for Baldwin, July 1982. 
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Table 4.5: Estimated trajectory-statistics at mean sea level for Pennsylvania, January 1982. 
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Table 4.6: Estimated trajectory-statistics at mean sea level for Pennsylvania, July 1983. 
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Figure 4.1. Points of mean sea level pressure. 
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I, GENERATION OF 850 mb TRAJECTORIES 

The data for the upper air version of the Lagrangian Model 
are t^ken from an objectively analysed grid of 850 mb geopotential 
height. This grid consists of 121 regularly spaced points and extends 
parallel to the x-axis from coordinates 20 to 30, inclusive, and paral- 
lel to the y-axis from coordinates 10 to 20, inclusive. The geopoten- 
tial height observations are available for OOUT and I2UT only. 

For the 850 mb trajectories it was possible to use only a 
time-step of 12 h. If data were unavailable at a time-step then the 
trajectories would end at the previously calculated point where data 
were available. Again, a trajectory would not be reported if it failed 
to attain an elapsed travel-time of at least 24 h. Trajectories were 
generated for 96 h of travel -time. Trajectories are initiated every 
12 h beginning at OOUT on the first day of each month. As with the 
surface data each monthly data set includes two days from the following 
month. So the last time of trajectory initiation is 12 h of the first 
day of either February or August. 

Statistics estimated at each travel -time for trajectories 
initiated from Sudbury are reported in Table 5.1 for January 1982 and in 
Table 5.2 for July of the same year. Similar statistics for trajec- 
tories initiated from Baldwin are listed in Tables 5.3 and 5.4 for 
January and July, respectively. In Tables 5.5 and 5.6 are the 



Concord Sctentitic Cofporatioi 



32 



statistics estimated from trajectories generated from the source in 
Pennsylvania for January and July respectively. 

Some general comments can be made about these results. There 
are significantly fewer trajectories initiated at 850 mb then at the 
mean sea level because at 850 mb observations are available only twice 
a day, whereas the mean sea level trajectory-generation used reported 
data at 3h intervals. As with the surface trajectories, there are fewer 
trajectories initiated in January than in July because of more missing 
data in the month of January than July. It is comforting to note that 
the estimates of the second order moments are of the same order of mag- 
nitude for 850 mb as for the surface. But the angles of the major axis 
of the probability ellipse to the axis of the coordinate system are 
larger and show more variability in time as the trajectory progresses 
for the 850 mb level. 
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Table 5.1: Estimated trajectory-statistics at 850 mb for Sudbury, January 1982. 
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Table 5.2: Estimated trajectory-statistics at 850 mb for Sudbury, July 1982. 






SnuBCr <;T»TIOM 1st I - ONT SlIOHLlQir 

TIME rftfGHI » OFH»|*»lmi 



ir> 


«*• 


lo«.o 


1?. 


M* 


loft.o 


Zo. 


M* 


»5.3 


56. 


*«* 


7b. b 


48. 


Sft. 


Sb.3 


60. 


ZJ. 


35.9 


ri. 


l*. 


25.0 


R*. 


tlf 


17.2 


9ft. 


•» 


U.S 



PE9I001 B70f 1EIr,^Il «5n uimEIOhteo 



ItMFAM 



?S.T90 

?7. «.■>»> 
?T.9JI» 

?M.|03 
af.938 



0.0 

o.on 
o.iin 

0.1*'* 
0.?09 
O.Jfl^ 
0.<K»T 
0.590 
0.5I9 



15.560 
Ift.OM 
Ih.itiO 
16.«|T 
lf).f«3 
16.?6T 
15.912 
I5.*b3 
|'..59* 



0.0 

0.07T 

0.l2ft 

0.1 81 

0.?39 

0.13* 

0.»M 

0.693 

0.H2* 



0.0 

n.58» 

O.Sftl 
1.0*7 
1.255 

i.r«6 

1.95^ 
1.T5I 



0.0 

0.619 

0.98* 

l.?6r 

l.*36 
l.60» 
1.910 
2.?98 
2.332 



0.0 

O.JOT 

• .66* 

l.«6« 

I.*»* 

t.98» 

2.8)6 

3.9JI 

3.288 



«.• 
62.8 

61.8 
66.8 

63.5 
62.3 
61. » 
65. T 
6T,0 



Table 5.3: Estimated trajectory-statistics at 850 mb for Baldwin, January 1982. 
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Table 5.4: 



Estimated trajectory-statistics at 850 mb for Baldwin, July 1982. 
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Table 5.5: Estimated trajectory-statistics at 850 mb for Pennsylvania, January 1982. 



50U«CF >!f»TlON 1st 1 - "It COyxJO »I.LEb 





56. 


loo.o 




<iB. 


B5.7 




2t. 


• «.? 




11. 


19.6 




B. 


1*.3 




*« 


10.' 




Mm 


5.* 




1*. 


$.« 




l'». 


1^. 



PflOUI "^Ol 


HEIGHTJ *5n 


UVKEIGHTfO 








• "tin 


YMEiN 


»S1G«« 


TSIO* 


COVMUWC 


T«*CT« 


?6.I60 * 0.0 


11.650 1 0.0 


0.0 


0.0 




1,1 


^ft.T^l • o.ibf 


i].<>n • 0.1*1 


l.lSi 


0.9M 




33.9 


?6.>«J9 • 0.?9% 


1 J. J"** • o.^6? 


I.52T 


1 . 160 




?••! 


^T.l^^ • o.<»3« 


ll.il? t 0.3M 


l.»51 


1.236 




?).• 


?T,|76 . 0.5?'> 


n.v<>4 • o.soi 


l.kn« 


l.«>|S 




1«.T 


^7.?^<^ • o.i69 


t3.H53 ♦. O.AfS 


1.39 J 


1.65* 




*•• 


26.558 t O.f^O 


13.263 t 0.V5T 


1 . 3SS 


0.T92 


.0 27* 


-■•1 


IT.O'V * 0.19? 


1J.J«9 1 0.»99 


1.5*5 


ft.B6* 


*1 .035 


-23.* 


?ft,5l8 1 0. J<.5 


l».?J5 • O.flB? 


o.*«^ 


0.116 


-0.05T 


-13.* 



■mig 



^) 



Table 5.6: Estimated trajectory-statistics at 850 mb for Pennsylvania, July 1982. 
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i, COMPARISON OF 850 mb AND MEAN SEA LEVEL TRAJECTORIES 

To aid in the comparison of the statistical properties of 850 
mb and mean sea level trajectories, it was decided to construct a 
filtered set (as it came to be called) of trajectories at each level by 
including In each set only those trajectories that have a counterpart in 
the other set that is initiated at the same time and date. Hence, for 
example, the filtered set of sea level trajectories in January 1982 for 
Sudbury has the same number of trajectories at the elapsed travel time 
of Oh as the filtered set of 850 mb trajectories for the same period. 

The estimated statistics for the filtered set of trajectories 
for the two levels, two time periods and three sources are given in 
Tables 6.1 through 6.12. The captions identify each table clearly, so 
the identification of each table will not be repeated here. 

Given the greater density of surface observing stations as 
compared to that of upper air stations, one expects fewer trajectories 
to survive to 96 h of elapsed time at 850 mb than at mean sea level 
because of the increased likelihood of convergence criterion failure 
resulting from poorer representation of the 850 mb geopotential height 
field. Also, one cannot neglect the fact that 850 mb flow has a strong 
westerly component and therefore more trajectories would be expected to 
reach the boundaries of the region of interest of 850 mb than at mean 
sea level. Indeed for Sudbury and the source in Pennsylvania, fewer 
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trajectories survive at 850 mb than at mean sea level. Interestingly 
enough more trajectories survive the 96h at the 850 mb level than at the 
surface for Baldwin, although the numbers are certainly comparable. 
This counter-intuitive result can be due to the fact that Baldwin is the 
westernmost location of the three sources. 

To really appreciate these statistics, it is best to display 
them on a map of North America so that one can see the travel and vari- 
ability of these trajectories in relation to geography. Figures 6.1 and 
6.2 contain the means of trajectory points originating from Sudbury for 
the period January 1982 for the surface and 850 mb levels, respectively. 
The pattern of mean positions for the surface is rather erratic due to 
the small number of trajectories constituting the average. With more 
ensembles of trajectories, the pattern would be more established and 
regular (as was the case discussed in the main text of this report). In 
all the figures of this type, a square is drawn to indicate the source 
location, and triangles are used for the mean location at each time step 
with line segments connecting the points to help illustrate the time 
evolution of the mean surface trajectory. For Sudbury in January of 
1982, the mean surface trajectory is disorderly and is largely confined 
to southern Ontario. The 850 mb mean trajectory starts with an antici- 
pated northeastward motion but soon veers southward and displays a 
disorderly pattern in the states of New York and Pennsylvania. 
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The mean trajectories in January 1982 at the surface and 850 
mb for Baldwin are depicted in Figures 6.3 and 6.4, respectively. These 
mean trajectories are starkly different. The surface mean trajectory 
proceeds toward Manitoba where it ends. The 850 mb mean proceeds north- 
eastward and reaches New York State at the 36 h mark but then it 
doubles-back to complete its journey at 96 h in Ohio. This doubling- 
back is curious and might result from lack of resolution of the geo- 
physical height surface over the ocean or near-boundary effects - either 
reason accounting for trajectories terminating before reaching an 
elapsed time of 96h. 

Figures 6.5 and 6.6 depict mean trajectories originating from 
Sudbury for the averaging period of July 1982 at the surface and 850 mb, 
respectively. Both mean trajectories for the first 36h move in a north- 
eastward direction through southern Quebec. After that, they go their 
own separate ways. The surface mean trajectory shows an erratic pattern 
after 36h, which is thought to be symptomatic of small sample size for 
these travel -times. The 850 mb mean trajectory after 36h proceeds 
nearly parallel to the coastline through southwestern Quebec to 
Maryland. This too most likely results from small sanple size and 
additionally near boundary effects. 

The mean trajectories calculated from Baldwin, using the July 
data, for the surface and 850 mb are illustrated in Figures 6.7 and 6.8, 
respectively. Both proceed in a northeastward direction. The surface 




Concord Scientidc Corporatio 



- 42 - 



mean trajectory reaches past SauU St. Marie to approch Sudbury. The 
850 mb mean trajectory finishes its 96h journey in the tri-state area of 
New York, New Jersey, and Pennsylvania. These mean trajectories to- 
gether .show a remarkably coherent motion. They correspond quite well to 
what one would intuitively expect to see. 

The last four figures of this part of the report comprise an 
illustration of mean trajectories together with the associated standard 
errors in the x and y directions (see main text. Section 4). Unfortu- 
nately, for display purposes it was not possible to give the standard 
errors at each time step, or even at time steps separated by 12 h. An 
uncluttered and clear picture is possible with standard errors depicted 
as line segments at time steps separated by 24 h. As can be appreciated 
from the proceeding figures of mean trajectories for Sudbury, it was not 
possible to produce informative figures with standard errors for the 
surface. The mean trajectories without standard errors for the source 
in Pennsylvania were not presented in this report for they were diffi- 
cult to follow. So the mean trajectories with standard errors every 24 
h, but with each travel point plotted, for Baldwin are illustrated in 
Figures 6.9 through 6.12. In these figures a line segment that is the 
length of two standard errors of the x-coordinate is drawn so that the 
trajectory point is the mid-point of the line segment, and the same sort 
of operation is done to display the standard error of the y-coordinate. 
It is certainly reasonable for a physical scientist to look at these 
figures and think that he is viewing error bars in the x and y direc- 
tions. 
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Figure 6.11, which displays information from the surface 
calculations in January, shows increasing standard errors from 24 h to 
48 h but then the standard errors shrink for 72 h and 96 h. This is a 
small sample effect for at 72 h only 6 of 45 initial trajectories 
remain. Similarly, for the same average period for the 850 mb trajec- 
tories, the standard errors grow from 24 h to 48 h but then the antici- 
pated growth is not seen for 72 h or 96 h. Small sample size is 
probably the cause since only 5 of the 45 trajectories survive at 72 h. 

We are more fortunate with the July trajectories of the 
filtered set for Baldwin because out of 45 initial trajectories, 23 
trajectories at the surface survive at 96 h, and 26 trajectories at 850 
mb survive. The initial sample size for both months is small but at 
least we have the consolation that in July we don't lose all but approx- 
imately 5 % as in the case in January. The standard errors for Baldwin 
in July increase as time progresses at both the surface and 850 mb. On 
comparing Figures 6.11 and 6.12, one notices that the standard errors 
are larger for the surface trajectories at all time-steps except 24 h. 
The standard errors for both levels at 96h are still roughly comparable. 
If one were to superimpose Figure 6.12 over Figure 6.11, one would find 
no overlap of rectangular regions defined by standard errors from the 
two different levels at the same elapsed time. 
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Table 6.1: Estimated statistics for filtered set of trajectories at mean sea level for Sudbury. 
Januarv 1982. 
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Table 6.2: Estimated statistics for filtered set of trajectories at mean sea level for Sudbury. 
July 19R2. 
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Table 6.3: Estimated statistics for filtered set of trajectories at mean sea level for Baldwin. 
January 19R2, 
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6.4: Estimated statistics for filtered set of trajectories at mean sea level for Baldwin. 



July 1982. 
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Table 6.5: Estimated statistics for filtered set of trajectories at mean sea level for Pennsylvania, 
January 1982. 
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Table 6.6: Estimated statistics for filtered trajectories at mean sea level for Pennsylvania. 
July 19R?. 
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Table 6.7: 



Estimated statistics for filtered trajectories at 850 mb for Sudbury. January 1982. 
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Table 6.8: Estimated statistics for filtered trajectories at 850 nib for Sudbury. July 1982. 
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Table 6.9: Estimated statistics for filtered trajectories at 850 mb for Baldwin. January 1982. 
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Table 6.10: Estimated statistics for filtered trajectories at 050 mb for Baldwin. July 1982. 
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Table 6.11: Estimated statistics for filtered trajectories at 850 mb for Pennsylvania. January 1982, 
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Table 6.12: Estimated statistics for filtered trajectories at 850 mb for Pennsylvania. July 1982, 
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Figure 6.1: Means of the trajectory points at mean sea level for 
Sudbury. January 1982. 
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Figure 6.2: Means of the trajectory points at 850 mb for Sudbury, 
January 1982. 
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Figure 6.3: Means of the trajectory points at mean sea level for 
Baldwin. January 1982. 
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Figure 6.4: Means of the trajectory points at 850 mb for Baldwin, 
January 1982. 
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Figure 6.5: Means of the trajectory points at mean sea level for Sudbury, 
July 1982, 
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Figure 6.6: Means of the trajectory points at 850 rab for Sudbury. July 1982. 
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Figure 6.7: Means of the trajectory points at mean sea level for 
Baldwin. July 1982. 
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Figure 6.8: Means of the trajectory points at 850 mb for Baldwin. July 1982. 
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Figure 6.9: Mean trajectories with standard errors for Baldwin at mean 
sea level. January 1982. 
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Figure 6.10: Mean trajectories with standard errors for Baldwin at 
850 mb. January 1982. 
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Figure 6.11: Mean trajectories with standard errors for Baldwin at mean 
sea level. July 1982. 
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Figure 6.12: Mean trajectories with standard errors for Baldwin at 850 mb. 
July 1982. 
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